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of  turbulence  and  intemal  waves  (AIDV)  as  a  function  of  wfA”**  and  deptli. 

The  diurnal  change  in  IDV  was  pronounced  at  both  the  35-m  and  the  60-m  depths,  with  nighttiine  IDV  significaittiy 

gr^  than  daytime  IDV.  This  result  is  consistent  with  TROPIC  HEAT  mictostructute  observations  over  two  periods  rf 
a  few  days  to  weeks,  showing  that  IDV  is  a  useful  proxy  for  tnrbuleace  and  internal  wave  energy  dnriim  periods  when 
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ABSTRACT 


The  diurnal  cycle  of  high-frequency  temperahire  fluctuations  attributable  to 
turbulence  and  internal  waves  in  the  upper  central  equatorial  Pacific  Ocean  is  found  to 
persist  over  most  of  the  period  November  1983  •  October  1987.  Moored  temperature  and 
velocity  data  in  the  upper  300  m  at  0°.  140*^  during  this  four-year  period  were  used  to: 
1)  determine  the  vertical  extent  of  the  diurnal  cycle  of  turbulence  and  internal  waves  and 
evaluate  its  contributicm  to  the  equatorial  zonal  momentum  balance,  and  2)  determine 
whether  die  diurnal  cycle  of  turbulence  and  internal  waves  varied  seasonally  and 
interannually  in  resptMise  to  varying  surface  forcing  and  ambient  stability  conditimis. 
Using  15-minute  spot-sampled  temperature  data  at  35  m,  60  m,  100  m.  140  m,  200  m, 
and  300  m,  isotherm  displacement  variance  (IDV)  was  c<»nputed  for  each  six-hour  time 
period  as  a  proxy  for  vertical  overturning  and  displacement  associated  with  turbulence  and 
internal  waves.  The  time  series  of  monthly  mean  IDV  for  each  of  the  six-hour  time  bins 
showed  the  magnitude  of  the  diurnal  cycle  of  turbulence  and  internal  waves  (AIDV)  as 
a  function  of  season  and  dqith. 

The  diurnal  change  in  IDV  was  pronounced  at  both  die  35-m  and  die  60-m  dqiths, 
with  nighttime  IDV  significandy  greater  than  daytime  IDV.  This  result  is  ctmsistent  with 
TROPIC  HEAT  microstructure  observations  over  two  periods  of  a  few  days  to  wedcs, 
showing  that  IDV  is  a  useful  proxy  for  turbulence  and  internal  wave  energy  during 
periods  when  microstructure  measurements  are  not  available.  The  magnitutte  of  AIDV 
decayed  sharply  with  depth  below  60  m,  indicating  that  little  diumally-modulated  internal 
gravity  wave  energy  had  propagated  down  throu^  the  equatorial  undercurrent  core. 
Hence,  internal  wave  fHxipagation  does  not  appear  to  be  important  for  loss  of  zonal 
momentum  to  the  regions  below  the  core  of  the  undercurrent. 

Monthly  means  of  both  IDV  and  AIDV  at  35  m  and  60  m  suggest  that  turbulence 
and  internal  waves  varied  over  seasonal  and  interannual  time  scales.  Seasonally,  maxima 
of  IDV  and  AIDV  occurred  during  the  boreal  winter,  when  wind  speed  and  mixed-layer 
depth  woe  at  their  armual  maxima  and  shortwave  radiation  was  at  its  armual  minimum. 
Minima  of  IDV  and  AIDV  occurred  dining  the  boreal  spring,  when  wind  speed  and 
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mixed-layer  depth  were  at  their  annual  minima  and  shortwave  radiation  was  at  its  annual 
maximum.  Over  the  48-month  period  of  study.  AIDV  at  35  m  and  60  m  were  significantly 
correlated  with  wind  speed,  zmial  wind,  mixed-layer  depth,  diurnal  range  of  mixed-layer 
depth,  buoyancy  frequency,  and  shortwave  radiation.  These  strong  correlations  are 
consistent  with  the  use  of  IDV  as  a  proxy  for  turbulence-wave  processes  associated  with 
shear  production  and  buoyancy  flux  in  the  marginally  stable  entrainment  zone  between 
the  surface  well-mixed  layer  and  the  undercurrent  ctxe. 
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I.  INTRODUCTION 


A.  INTRODUCTION 

In  recent  years,  we  have  become  increasingly  aware  of  the  variable  nature  of  the 
earth’s  climate  and  of  the  need  to  predict  climatic  variations  arising  from  natural  and 
man-made  causes.  The  large  El  Nino/Southem  Oscillation  (ENSO)  event  of  1982-83 
focused  attention  on  understanding  the  interaction  of  the  ocean  and  atmosphere  in  the 
equatorial  Pacific  Ocean  (Cane.  1983;  Rasmussen  and  Wallace,  1983).  The  Southern 
Oscillation,  an  aperiodic  interannual  fluctuation  between  warm  and  cold  oceanic  states, 
has  its  largest  signature  in  and  over  the  trqjical  Pacific  and  Indian  Oceans,  but  it  is 
known  to  affect  oceanic  and  atmospheric  conditions  globally  (Philander,  1990).  Numerical 
modeling  experiments  have  shown  that  most  of  the  interannual  variability  of  the  tropical 
atmosphere  is  related  to  the  changes  in  the  boundary  conditions  at  the  ocean  surface, 
partictilarly  sea  surface  temperature  (SST)  fluctuations  associated  with  El  Nino  (Lau, 
1981,  Shukla  and  Wallace,  1983).  SST  anmnalies  in  die  equatorial  Pacific  modify  the 
fluxes  of  heat  and  moisture  between  the  ocean  and  atmosphere  and  have  been  linked  to 
displacement  of  large*scale  atmospheric  circulation  patterns  over  the  American  continent 
and  elsewhere  (Bjerknes,  1969;  Wyrtki,  1975;  Hotel  and  Wallace,  1981;  TVenberth  et  al., 
1988). 

In  addition  to  the  direct  physical  role  the  equatorial  upper  ocean  plays  in 
modulating  global  climate,  through  photosynthesis  it  also  plays  an  essential  role  in  the 
cycling  of  oxygen,  carbtm  dioxide  (CO,),  and  other  atmospheric  constituents  through  the 
global  system.  In  this  manner,  the  u{^)er  ocean  regulates  the  levels  of  greenhouse  gases 
in  the  atmosphere,  most  notably  CO,,  and  hence  global  climate  (Takahashi,  1989).  In 
particular,  upwelling  of  CO,-rich  subsurface  waters  make  die  equatorial  Padfic  the  most 
important  oceanic  CO,  source  in  the  oceans  (Takahashi  et  al.,  1992).  Although  the  source 
intensity  of  CO,  in  the  equatorial  Pacific  seasonally  varies  a  relatively  minor  amount 
(about  25%)  during  non-El  Nino  years,  diis  source  is  eliminated  by  the  blanketing  effect 
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of  warm.  low-COj  surface  water  that  spreads  eastward  from  the  western  Pacific  during 
El  Nino  warm  events.  Similarly,  the  upper  equatCHial  PaciEc  constitutes  an  important  link 
in  the  oceanic  food  chain.  These  pervasive  influences  provide  the  motivation  for  studying 
the  upper  equatorial  Pacific  Ocean  (Kraus.  1977). 

The  ocean  circulation  is  driven  at  the  surface  by  wind  stress  and  by  buoyancy 
forcing  due  to  heating,  cooling,  evaporation,  and  precipitation.  This  direct  forcing  by  the 
wind  stress  and  the  buoyancy  flux  is  generally  absorbed  within  a  surface  layer  that  is 
often  nearly  homt^eneous  in  temperature,  salinity,  and  density  due  to  turbulent  mixing, 
and  is  thereftve  commonly  referred  to  as  the  ocean  surface  mixed  layer.  The  thickness  of 
this  mixed  layer  ranges  from  near  zero  during  periods  of  light  winds  and  strong  heating 
to  greater  than  a  hundred  meters  under  conditions  of  strong  cooling  and  high  winds.  The 
direct  efrect  of  surface  forcing  is  not  always  limited  to  the  mixed  layer,  however,  but 
sometimes  extends  into  the  density-stratifred  region  below  the  mixed  layer.  This  transition 
region  of  surface  influence,  in  which  surface  layer  turbulence  provides  die  energy  to 
entrain  the  stratified  fluid  below,  is  often  referred  to  as  the  entrainment  zone. 

TROPIC  HEAT  experiment  (Eriksen,  1985)  miaostructure  measurements  in  the 
central  equat<xial  Pacific  in  November  1984  and  April  1987  showed  a  pronounced  diurnal 
cycle  in  vertical  turbulent  mixing  rate  in  both  the  surface  mixed  layer  and  in  the 
stratified  region  below  the  mixed  layer  (Gregg  et  al.,  1985;  Mourn  and  Caldwell,  1985; 
Mourn  et  al.,  1992a).  Since  the  turbulmit  dissipation  varied  diumally  at  d^tfas 
significantly  greater  than  the  penetration  of  die  surface  well-mixed  layer,  turbulent  mixing 
appears  to  have  been  caused  by  processes  other  dian  cmiventional  mixed-layer  physics 
(Mourn  and  Caldwell,  1985;  Mourn  et  al.,  1989;  Peters  et  al.,  1988).  Since  the  mily 
fordng  having  strong  and  persistent  ntght-to-day  differences  was  the  solar  instdation 
(Mourn  and  Caldwell,  1985),  the  diurnal  cycle  in  dissqMtion  presumably  must  have  been 
modulated  by  the  diumally-varying  surface  buoyancy  flux. 

Although  the  mechanisms  responsible  for  the  diurnal  cycle  of  turbulent  diss4>ati(Hi 
in  the  pycnocline  remain  uncertain,  it  has  been  suggested  that  internal  gravity  waves 
propagating  downward  from  the  base  of  the  mixed  layer  to  die  highly  sheared  pycnodine 
might  be  responsible  for  the  qioradic  de^  mixing  (Gregg  et  al.,  1985;  Toole  et  al.,  1987; 
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Peters  et  al.,  1988;  Dillon  et  al..  1989;  Mourn  et  al.,  1989).  Some  of  these  investigations 
reported  a  diurnal  cycle  of  high-frequency  wavelike  motions  coincident  with  the  diurnal 
cycle  of  dissipation  (Mourn  et  al..  1992b).  According  to  this  hypothesis,  the  internal 
waves  provide  the  additional  energy  necessary  to  initiate  dynamic  instabilities  and 
turbulent  mixing.  Recently.  Wijesekera  and  Dillon  (1991),  McPhaden  and  Peters  (1992), 
Mourn  et  al.  (1992a),  and  Hebert  et  al.  (1992)  have  presented  evidence  from 
microstructure  proftlers,  moored  buoys,  and  towed  thermistor  chains  that  is  consistent 
with  this  hypothesis.  These  studies  present  observations  of  a  distinct,  narrow  frequency- 
band,  and  diumally-varying  signal  in  temperature  records  obtained  in  the  low  Richardson 
number  shear  flow  above  the  Equatorial  Undercurrent  (EUC)  core  (Mourn  et  al.  1992b). 
The  moored  buoy  temperature  records  showed  diumally-varying  high-frequency  wavelike 
motions,  possibly  internal  waves,  where  the  measured  frequency  of  10-30  cph  in  tte 
presence  of  the  strong  mean  flow  of  the  equatorial  undercurrent  (EUC)  implied  an 
intrinsic  frequency  near  the  buoyancy  frequency  of  N  =  5-7  q>h  (McPhaden  and  Peters, 
1992).  The  zonal  wavelength  was  iq)proximately  150-250  m  (Mourn  et  al.,  1992a.b).  It 
is  possible  that  the  observed  fluctuations  represented  a  forced  mode  resulting  from  local 
shear  instability  of  the  mean  flow  or  a  combination  of  internal  waves  and  shear 
instabilities. 

Skyllingstad  and  Denbo  (1994)  used  a  two-dimensional,  non-hydrostatic.  large 
eddy  simulation  model  to  examine  the  hypothesized  diurnal  cycle  of  downward 
propagating  internal  waves  in  the  central  equatorial  Pacific.  Their  model  was  initialized 
and  forced  with  temperature,  velocity,  and  wind  observations  from  TROPIC  HEAT.  Their 
model  results  were  consistent  with  the  TROPIC  HEAT  microstructure  observations  and 
with  the  general  hypothesis  that  internal  waves,  generated  in  or  at  the  base  of  the 
diumally-varying  surface  mixed  layer,  propagate  downward  into  the  stratified,  high-shear 
flow  above  the  EUC.  A  more  thorough  literature  review  of  equatorial  internal  waves  and 
turbulent  mixing  is  presented  in  Chiq)ter  n. 

The  results  from  the  TROPIC  HEAT  observations  and  numerical  modeling  studies 
to  date  have  he^)ed  describe  the  diurnal  cycte  of  turbulent  mixing  in  the  central  equatorial 
Pacific  at  the  times  of  those  observations.  Unfortunately,  these  observations  are 
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temporally  scarce,  being  limited  to  a  few  occasions  of  intense  sampling  over  periods 
ranging  from  a  few  days  to  a  few  weeks.  This  scarcity  of  microstructure  obsoA^ations 
limits  the  ability  to  examine  the  diurnal  cycle  of  internal  waves  and  turbulent  mixing 
under  a  mcne  complete  range  of  surface  forcing  and  dynamic  stability  conditions.  Both 
surface  forcing  and  dynamic  stability  conditions  in  the  central  equatorial  Pacific  are 
known  to  vary  on  many  different  temporal  scales,  from  diurnal  to  seasonal  and 
interannual.  Although  the  available  observations  have  suggested  possible  relationships 
among  the  diurnal  cycle  of  turbulent  mixing  in  the  stradfred  region  and  the  local  wind 
speed,  wind  direction,  and  vertical  velocity  shear  in  the  water  column  (Mourn  et  al., 
1992a;  McPhaden  and  Peters,  1992;  Skyllingstad  and  Denbo,  1994),  the  short  time  span 
of  the  observations  limits  our  ability  to  derive  generalized  relationships  between  forcing 
and  internal  mixing. 

In  addition  to  the  scarcity  of  microstructure  data,  the  TROPIC  HEAT 
microstructure  observations  were  limited  in  their  vertical  extmit  to  the  upper  120  m  or 
so,  roughly  the  dq>th  of  the  EUC  core.  Although  this  dq>th  range  q>pears  adequate  to 
examine  the  directly  forced  mixed-layer  processes  and  the  entrainment  zone  between  the 
mixed  layer  and  the  EUC  core,  it  has  not  allowed  adequate  investigation  of  the  vertical 
extent  of  the  diurnal  cycle  of  the  internal  wavefield.  In  particular,  prior  studies  have  not 
been  able  to  determine  whether  or  not  these  waves  transport  momentum  and  energy  from 
the  surface  generation  region  through  the  EUC  to  the  high  shear  flow  beneadi  the  core. 
The  resolution  of  this  issue  is  important  to  explain  the  equatorial  mtHnentum  balance, 
which  is  discussed  in  further  detail  in  Chapter  n. 

B.  RESEARCH  OBJECTIVES  AND  STRUCTURE  OF  DISSERTATION 

This  dissertation  ctmcems  the  vertical  extent  and  time  variability  of  the  diurnal 
cycle  of  high-frequency  temperature  fluctuations  in  the  central  equatorial  Pacific.  Unlike 
priw  investigaticms  that  examined  the  diurnal  cycle  of  turbulence  and  internal  waves 
using  a  few  short  periods  of  high  temporal  and  vertical  resolution  microstructure  data,  this 
study  examines  the  reqxmse  of  the  mixed  layer  and  internal  wave  field  to  the  time 
varying  background  conditicms  ami  forcing  over  a  four-year  period  using  moored  buoy 
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data.  Specifically,  moored  buoy  data  (McPhaden  and  Taft,  1988;  McPhaden  and  Hayes, 
1990)  from  0",  140^  for  the  period,  November  1983  -  October  1987,  are  utilized  to 
investigate  the  diurnal  cycle  of  the  high-frequency  temperature  variability,  as  described 
by  isotherm  displacement  variance  (IDV),  under  varying  surface  facing,  stratification 
and  velocity  conditions.  These  observations  of  IDV  could  reflect  high-frequency  internal 
gravity  waves,  local  shear  instabilities,  turbulence,  or  some  crunbination  of  these 
processes.  The  magnitude  of  the  IDV  may  be  an  effective  measure  of  all  of  these 
processes,  but  it  generally  is  not  a  means  of  distinguishing  among  them. 

While  the  vertical  range  of  data  sets  used  in  [Mior  studies  extended  only  to  about 
120  meters,  the  mooring  data  span  from  the  surface  to  300  meters,  allowing  a  more 
cmnplete  examination  of  the  vertical  structure  of  the  diumally-varying  high-frequency 
temperature  fluctuations  and  the  relationship  between  these  wavelike  motions  and  the 
EUC.  Statistical  relationships  between  the  vtftical  structure  of  the  diurnal  cycle  of  high- 
frequency  temperature  variability  and  surface  winds,  the  surface  buoyancy  flux,  and  the 
dynamic  stability  are  found.  Furthermore,  the  contribution  of  the  diumally  varying  high- 
frequency  wavefield  to  the  heat  and  momentum  balances  altmg  the  equator  are  discussed. 

In  Chiq)ter  n,  a  more  complete  review  of  the  relevant  literature  concerning 
equatorial  ui^)er  ocean  dynamics  is  presented,  with  particular  regard  to  high-frequency 
internal  waves  and  turbulent  mixing.  Both  analyses  of  observations  and  numerical 
modeling  studies  are  reviewed.  In  Chapter  m,  the  data  sets  are  described  and  data 
processing  methods  are  discussed.  The  procedures  for  computing  IDV  and  the  diurnal 
cycle  of  IDV  are  discussed.  This  is  followed  by  a  description  of  a  sensitivity  study  to 
validate  this  procedure.  In  Cluster  IV,  the  seastmal  and  interannual  variability  of  the 
baclcground  oceanic  and  atmospheric  conditions  during  the  period  of  study  are  {X^nted. 
In  Chapter  V,  the  seasonal  and  interannual  variability  of  high-frequency  IDV  is 
discussed,  and  relationships  between  IDV  and  surface  forcing,  dynamic  stability,  and 
mixed-layo’  dq>th  are  examined.  In  Chapter  VI,  the  focus  is  on  the  diurnal  cycle  of  IDV. 
With  only  a  few  exceptions,  it  will  be  shown  that  the  diurnal  cycle  of  turbulence  and 
internal  waves,  as  indicated  by  IDV,  was  confined  above  tlw  core  of  the  EUC  and  ustially 
above  the  100  m  dq>th.  Relationships  between  the  diurnal  cycle  of  IDV  and  surface 
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forcing  and  dynamic  stability  parameters,  as  well  as  diurnal  cycles  of  these  parameters 
are  discussed.  This  is  followed  by  an  examinatitm  of  the  relationsh4>s  between  the  diurnal 
cycle  of  IDV  and  the  mixed-layer  depth  and  the  diurnal  range  of  the  mixed-layer  depth. 
The  interpretation  of  IDV  in  terms  of  turbulence  and  high-frequency  internal  waves  is 
further  examined  by  investigating  the  diurnal  variability  of  the  wave  momentum  flux  for 
the  month  of  November  1984.  The  ctxnbined  analyses  suggest  that  IDV  is  related  to 
turbulence  w  high-frequency  internal  waves  transporting  hmizontal  momentum  from  the 
surface  layer  to  the  stratified  region,  but  not  through  the  EUC  ewe.  The  duq>ter 
concludes  with  a  discussion  of  two  additional  possible  sources  of  variability,  mooring 
motion  and  gravitational  tides,  which  are  shown  to  be  unlikely  to  generate  the  (rf>served 
diurnal  cycle  of  IDV.  Quq)ter  vn  incluttes  a  discussion  of  the  results,  followed  by  the 
conclusicms  and  recommendations  fw  further  investigation. 
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II.  LITERATURE  REVIEW 


A.  EQUATORIAL  UPPER  OCEAN  DYNAMICS 

The  physical  processes  of  the  upper  equatorial  oceans  are  unique  due  to  the 
presence  of  the  large  vertical  shear  associated  with  the  EUC  and  the  vanishing  of 
vertical  comptment  of  planetary  rotation.  The  westward  trade  winds  almg  the  equatex* 
result  in  upwelling  of  0(3  m  day  ')  driven  by  Ekman  divergence  in  the  surface  waters  and 
geostrophic  convergence  deep&c  in  the  wator  column  (Wyitki,  1981;  Halpon  et  al.,  1989). 
This  westward  wind  stress  and  strong  upwelling  produce  a  tongue  of  relatively  cold 
surface  wato*  and  a  shallow  main  thermodine  in  the  eastern  equatorial  Pasific  and 
relatively  warm  surface  wato-  and  a  deq>  main  thermodine  in  the  western  equatorial 
Pacific.  This  upwelling  provides  a  source  of  cool  water  for  the  surface  layer.  The 
princqMd  cause  of  changes  in  die  depth  of  the  equatorial  diermodine  is  a  horiignntd 
redistribution  of  warm  surface  waters  in  re^xwse  to  the  large-scale  winds  (Wyitki.  1975; 
Merle.  1980;  Duing  and  Leetmaa.  1980;  Stevensen  and  Niiler.  1983;  Molinari  et  al.. 
1985).  Advecdon  strongly  influences  SST.  particularly  in  die  eastern  part  of  die  bnin. 
During  El  Nino,  for  example,  the  eastern  tropical  Pacific  is  warmed  by  die  eastward 
advecdon  of  warm  surface  waters  assodated  with  the  elevadon  of  die  thermodine  in  the 
west  and  its  deqiening  in  the  east.  This  advecdve  warming  competes  with  cooling  by 
iqiwelling.  Since  the  mean  vertical  vdodty  conqxment  vanishes  at  the  ocean  surfiKe. 
mixing  fMxicesses  determine  the  extent  to  which  iqmelling  affects  SST.  Nfixing  processes 
also  affect  the  intensity  of  the  equatorial  currents. 

The  steady-state  ztmal  momentum  budget  at  the  equator  indudes  die  ztmal 
convergence  of  ztmal  momentum,  iqiwelling  of  eastward  momentum,  die  zonal  pressure 
gradient,  the  divergence  of  the  mesoscale  eddy  flux,  and  die  turbulent  verdcal  stress 
(Bryden  and  Brady,  1985;  Dillcm  et  al.,  1989), 
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where  V  and  V  are  mean  zonal  (positive  eastward)  and  meridimud  (positive 
northward)  velocities,  P  is  mean  pressure,  p  is  density,  UV'  .  W'  ,  and  U'W 
are  the  zonal,  meridional,  and  vertical  turbulent  fluxes  of  zonal  mtxnentum,  and  t*  is  the 
zonal  wind  stress.  Subscripts  x,  y,  and  z  denote  partial  differentiation  in  the  zcmal  (positive 
eastward),  meridional  (positive  northward),  and  vertical  (positive  upward)  directions  . 
Analyses  of  observations  (McPhaden  and  Taft,  1988)  and  diagnostic  models  (Bryden  and 
Brady.  1985,1989)  indicate  tiiat  on  annual  time  scales  the  princ4>al  zonal  momentum 
balance  in  the  equatorial  uf^  ocean  is  between  the  vertically  integrated  z<»al  pressure 
gradient  force  and  the  wind  stress,  which  is  vertically  distributed  by  turbulent  mixing 
(Crawford  and  Osborn,  1981). 

Examination  of  the  terms  in  tiie  wtically  integrated  z(»al  rnmnentum  equations 
indicates  that  the  ztxial  pressure  gradient  is  not  balanced  by  the  surface  wind  stress  on 
shorter  time  scales  (Dilltm  et  al.,  1989;  Hebert  et  al.,  1991).  The  strength  of  the  wind 
plays  a  significant  role  in  determining  both  the  magnitucte  and,  likely,  tiie  mechanisms  of 
momentum  tranq)ort  between  the  surface  ami  tiie  core  (rf  tiie  EUC.  During  low  wind 
conditions  the  turbulent  stress  divergence  was  significantly  different  in  magnitude  and 
vertical  structure  from  tiiat  found  during  strong  winds  (Hdiert  et  al..  1991). 

In  November  1984,  the  measured  turbulent  stress  divergence  was  too  large  above 
40  m  to  balance  the  residual  term  in  the  zonal  momentum  budget  of  Bryden  and  Brady 
(1985,1989),  and  decayed  eiqionentially  witii  depth  from  tiie  wind  stress  value  at  the 
surface.  On  the  otha*  hand,  in  ^xil  1987  the  turbulent  stress  divergence  was  smaller  than 
that  required  by  Bryden  and  Brady  (1985,1989).  Using  data  from  tiie  NORPAX  Hawau- 
to-Tahiti  Shuttle  Experiment,  Johnsmi  and  Lutiier  (1994)  presented  evidmioe  for  a  rou^ 
zonal  momentum  balance  at  90*117  m,  even  neglecting  vertical  stresses.  Since  their  more 
directly  measured  ntimate  of  the  zonal  momentum  budget  has  a  vertical  structure  in  good 
agreement  witii  the  available,  dissqiation-derived  estimates  of  turbulent  vertical  stresses. 
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they  concluded  that  the  suggestion  of  an  imbalance  in  Bryden  and  Brady’s  momentum 
budget  below  the  thermocline  was  (a-obably  erroneous.  If  the  Johnson  and  Luther 
momentum  budget  is  correct,  the  proposed  requirement  for  vertically  prr^agating  waves 
to  balance  the  budget  (Dillon  et  al.  1989)  is  not  necessary. 

In  general,  momentum  input  by  wind  stress  is  vertically  distributed  by  turbulent 
mixing.  The  depth  of  mixing  is  dependent  on  the  static  and  dynamic  stability  of  the  water 
column.  This  vertical  distribution  of  momentum  in  the  upper  ocean  dictates  which  modes 
of  equabnial  Kelvin  and  Rossby  waves  are  forced,  and  governs  the  response  of  die  zonal 
pressure  gradient  and  the  EUC  to  changes  in  the  wind  stress  (Schudlich,  1991).  In 
summary,  the  vertical  profile  of  wind  stress  in  the  equatorial  ocean  is  an  imprxtant,  yet 
po<»rly  understood,  aspect  of  the  large-scale  equatorial  circulation. 

In  the  central  equatorial  Pacific,  the  South  Equatorial  Current  (SEC)  flows 
westward  at  speeds  averaging  about  0.3  ms  '.  On  the  equator,  beneath  the  SEC,  the 
eastward  flowing  EUC  attains  speeds  as  high  as  l.S  ms  '.  Resultant  vertical  shear  between 
the  surface  mixed  layer  and  the  EUC  c(»c  often  exceeds  0.02  s  ',  and  is  sufficiently  strrmg 
to  maintain  a  gradient  Richardson  number,  Ri,  of  less  than  unity  above  die  undercurrent 
core  (Gregg  et  al.,  1985;  Chereskin  et  al.,  1986;  Mourn  et  al.,  1989),  where  Ri^!^  (N 
is  the  local  buoyancy  frequency  and  S  is  the  magnitude  of  the  vertical  shear  of  the 
horizontal  current).  Theoretical  and  labm-atory  reaeardi  indicates  that  laminar  flow 
becomes  unstable  to  Kelvin-Helmholtz  wave  formation  and  the  onset  of  turbulence  when 
Ri  is  less  than  a  critical  Richardson  number,  Ri,,  of  0.25.  Once  turbulent,  however,  the 
turbulence  may  cemtinue  until  the  Ridiardson  number  reaches  a  terminal  value,  Ri,  of 
about  1.0  (Stull,  1988).  In  addition  to  a  low  Ri.  indicating  dynamic  instability,  the  tmset 
of  turbulence  often  requires  some  triggering  energy  (Stull,  1988).  Thus,  the  low  mean  Ri 
observed  in  the  high  shear  region  below  the  equatorial  mixed  layer  is  on  the  margin  of 
dynamic  stability  and  subject  to  the  onset  of  turbulence  widi  the  addition  of  arfHitinnai 
kinetic  energy. 

The  equatorial  mixed  layer  system  is  a  smnewhat  qiecial  case  with  r^;ard  to  the 
oceanic  planetary  boundary  layer  (Garwood  et  al.,  1989)  n  goieral,  the  CMiukhov  length 
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scale,  L  *  >  is  the  limiting  depA  for  oceanic  mixed  layers  under  the  influence 

"o 

of  wind  stress,  x,  and  buoyancy  flux  Bg  (Obukhov,  1946).  Complicating  the  equatorial 
mixed-layer  system,  upwelling  may  reduce  die  thickness  of  the  surface  layer  to  a  fraction 
of  the  Obukhov  length  (Muller  et  al.,  1984).  While  L  is  the  limiting  scale  fw  the  surface 
layer,  Garwood  et  al.  (1989)  proposed  that  Rossby  and  Montgmnery’s  (1935)  neutral 
planetary  boundary  layer  scale,  ujf,  where  u.  is  the  friction  velocity  and  /is  the  Coriolis 
parameto-,  may  be  an  iq>prq)riate  limit  for  the  thickness  of  the  entrainment  ztme,  5^.  Off 
the  equauv,  is  often  a  small  fraction  the  mixed-layer  depth,  and  therefore  has  been 
approximated  as  a  discontinuity  in  bulk  models  of  the  oceanic  mixed-layo-.  If  not  fmr 
advection,  the  entrainment  zone  at  the  equator  would  thicken  indefinitely.  With  typical 
mean  upwelling  conditions  along  the  equator,  Garwood  et  al.  (1989)  found  an  entrainment 
zone  thickness  of  about  100  m  under  a  shallow  mixed  layer,  in  good  agreement  with  the 
deep  diurnal  cycle  of  turbulence  observed  during  TROPIC  HEAT. 

The  rate  of  turbulent  kinetic  energy  dissipation  and  the  vertical  turbulent  fluxes 
of  heat,  mass,  and  momentum  are  seldom  measured  directly  but  can  be  inferred  from 
measurements  of  votical  temperature  gradients  and  horizontal  velocity.  Early 
measurements,  by  means  of  fiee-fsUing  instruments  that  were  tracked  very  accurately, 
revealed  that  turbulent  dis84>ation  is  higher  near  thenquator  than  in  extratrq>ical  latitudes. 
Dissipation  is  high  even  below  the  core  of  die  EUC  (Gregg,  1976;  Gregg  and  Sanford, 
1980;  Gaigett  and  Osborn,  1981),  {vesumably  a  consequence  of  the  large  vertical  shear. 
These  early  results  were  subsequently  used  to  devise  panuneterizatums  for  the  turbulent 
mixing  of  mmnentum  and  heat  in  terms  of  Ri  (Mourn  and  Caldwell,  1985;  Peters  et  al., 
1987).  When  used  in  numerical  models,  these  parameterizatioos  lead  to  simulatitms  that 
are  mwe  realistic  dum  those  with  constant  coefficients  of  eddy  viscosity  and  diffusivity 
(Pacanowski  and  Philander,  1981).  Althou^  these  parametorizadons  won  significant 
imivovements,  diey  remain  inadequate  fa-  satisfactory  modeling  of  the  tropical  Pacific 
(Philander,  1990). 
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B.  MICROSTRUCTURE  OBSERVATIONS  OF  THE  EQUATORIAL  UPPER 
OCEAN 

Microstructure  observstiotis  in  November  1984  near  0°,  140^,  as  part  of  the 
TROPIC  HEAT  Experiment,  revealed  a  large  diurnal  cycle  in  vertical  turbulent  mixing 
rates  (quantified  by  die  rate  of  viscous  dissipation  of  turbulent  kinetic  energy)  that  was 
previously  unknown  at  the  equator  (Gregg  et  al.,  1985;  Mourn  and  Caldwell.  1985;  Peters 
et  al..  1988;  Mourn  et  al.,  1989).  During  the  day,  when  strong  solar  heating  stabilized  the 
upper  10-20  m,  little  turbulent  dissipation  was  observed;  what  dissipation  existed  was 
limited  to  depths  less  than  about  10  m.  During  the  night,  when  surface  cooling  and  wind 
mixing  reduced  the  mixed-layer  temperature,  dissipation  increased  by  two  otders  of 
magnitude,  ft'om  10^  to  10^  W  kg  ’,  indicative  of  vigorous  mixing,  down  to  dqiths  of  80- 
90  m,  more  than  twice  the  depth  of  the  convectively  active  surface  mixed  layer,  which 
reached  to  about  35  m  at  night.  The  enhanced  dissipation  at  night  typically  occurred  in 
bursts  of  turbulence  lasting  2-3  hours.  Peak  diss4Mition  rates  in  the  pycnocline  lagged 
those  near  the  surface  by  several  hours  and  were  associated  with  individual  turbulent 
overturns  up  to  10  m  thick  (Mourn  et  al.,  1989). 

Observations  of  increased  mixing  in  tiie  surface  layers  of  the  ocean  at  night,  when 
a  loss  of  heat  to  the  atmosphere  destabilize  the  upper  ocean,  have  been  made  in  many 
locations.  What  is  unique  about  the  measurments  in  the  central  equatorial  Pacific  is  the 
intensity  of  the  turbulence  and  its  presrace  in  the  stratified  region  below  the  mixed  layer. 
In  cmitrast,  observations  from  the  upper  ocean  at  mid  latitudes  show  that  values  of 
dissipation  associated  with  turbulent  mixing  penetrate  only  about  10  m  below  the  mixed 
layer  at  night  even  when  the  mixed  layer  is  very  deep,  0(100  m)  (Shay  and  Gregg.  1986; 
Lombardo  and  Gregg,  1989;  Price  et  al.,  1986).  Although  the  equatorial  pycnocline  was 
well  stratified  (N  ~  5-9  cph),  the  high  vertical  shear  was  sufficient  to  maintain  Ri  between 
0.25  and  1.0,  sug^sting  marginal  dynamic  stability,  and  a  tendency  toward  the  formation 
of  turbulence.  Interestingly,  no  similar  diurnal  cycle  was  ajqwent  in  the  cmnputed  values 
of  Ri  (Chereskin  et  al.,  1986).  This  suggests  tiud  the  relationshq)  between  turbulent 
diss4>ation  and  local  Ri  may  not  be  simple  (Petos  et  al.,  1988).  It  is  possible,  however, 
that  an  important  diurnal  variatitm  in  Ri  mi^t  occur  on  vertical  scales  too  short  to  be 
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resolved  with  the  10-12  m  resolution  of  the  TROPIC  HEAT  data.  Preliminary  analysis 
of  more  recent  data  from  the  Tropical  Instability  Wave  Experiment  (TIWE)  in  December 
1991  suggest  a  possible  diurnal  cycle  of  Ri  in  the  mixed  layer  (Lien  et  al.,  1994); 
however,  these  data  require  further  analysis. 

C.  EQUATORIAL  INTERNAL  WAVES 
1.  Observations 

Although  the  physical  mechanism  for  the  diurnal  cycle  of  turbulent 
dissipation  in  the  pycnocline  remains  uncertain,  several  investigators  have  speculated  that 
internal  gravity  waves  might  play  an  imptxtant  role  (Gregg  et  al.,  1985;  Peters  et  al., 
1988;  Dillon  et  al.,  1989;  Mourn  et  al.,  1989;  Wijesekera  and  Dillon,  1991;  McPhaden 
and  Peters,  1992;  and  Mourn  et  al.,  1992a,b).  The  general  hypothesis  is  that  after  sunset, 
high-frequency  (near  AO  internal  waves  are  generated  just  below  the  mixed  layer  as  a 
result  of  turbulent  convective  overturns.  These  waves  propagate  downward,  become 
unstable,  and  generate  turbulent  eddies  above  the  undercurrent  c(xe  at  locations  where  Ri 
is  reduced  to  subciitical  values  by  superposition  of  wave  and  mean  current  velocities 
(Wijesekera  and  Dillon,  1991).  Dillon  et  al.  (1989)  and  Hebert  et  al.  (1991)  also 
suggested  that  internal  waves  were  likely  to  be  important  in  the  large-scale  zonal 
momentum  balance  of  the  pycnocline.  Both  of  their  turbulent  stress  divergence  estimates 
were  based  on  the  assumption  of  an  i^roximate  local  balance  between  production  and 
dissipation  of  turbulent  kinetic  energy,  neglecting  buoyancy  flux.  They  speculated  that  this 
local  balance  assumption  might  not  be  valid,  because  in  addition  to  vertical  diffusive 
fluxes  caused  by  local  shear  in?  oibility,  internal  waves  might  play  an  imptHlant  role  in 
equatorial  dynamics  by  the  downward  transport  of  zonal  momentum.  It  remains  imclear 
whedier  these  internal  waves,  if  generated,  transfer  momentum  through  the  EUC  into  the 
deep  ocean,  or  t»«ak  and  dissipate  in  the  high  shear  region  above  the  EUC  core.  Jonnson 
and  Luther  (1994)  indicated  that  although  such  mechanisms  may  occur,  they  are  not 
necessary  to  achieve  a  ztmal  momentum  balance  on  the  equatcx. 

In  thetx^cal  studies,  Booker  and  Bretherton  (1967),  Bredierton 
(1966,1969),  and  Muller  (1976)  demcmstrated  the  possibility  of  vertical  transport  of  zonal 
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momentum  by  internal  gravity  waves.  Boundary  layer  perturbations  can  generate  internal 
waves  that  may  transfer  energy  and  momentum  into  the  fluid  interior,  bypassing  the  usual 
turbulent-diffusion  mechanism.  Wijesekera  and  Dillon  (1991)  proposed  two  ways  in  which 
the  eddies  of  a  cooled  boundary  layer  could  result  in  the  excitation  of  internal  gravity 
waves  in  the  underlying  stable  layer,  thermal  forcing  and  shear  forcing  (the  obstacle 
effect).  They  concluded  that  in  the  central  equatorial  Pacific,  shear  facing  is  more 
efticient,  and  therefore  more  likely.  In  a  series  of  observational  and  numerical  modeling 
studies  of  the  atmospheric  planetary  boundary  layer.  Mason  and  Sykes  (1982),  Clark  et 
al.  (1986),  Kuettner  et  al.  (1987),  and  Huff  and  Clark  (1989)  found  that  boundary  layer 
eddies  are  able  to  impinge  upon  the  stratified  layer  above.  In  the  presence  of  vertical 
shear,  these  eddies  represent  a  localized  form  drag,  or  obstacle,  to  the  overlying  flow. 
These  impinging  eddies,  or  obstacles,  then  excite  internal  gravity  waves  (analogous  to  the 
generation  of  mountain  lee  waves).  Similarly  in  the  oceanic  mixed  layer,  Wijesekera  and 
Dillon  (1991)  hypothesize  the  generation  of  downward  prt^agating  internal  gravity  waves 
at  the  base  of  the  mixed  layer  with  a  horizontal  phase  velocity  equal  to  the  average 
advection  velocity  of  the  mixed-layer  eddies.  Based  on  the  average  westward  advection 
of  the  mixed  layer  during  their  November  1984  observation  pmod,  they  concluded 
westward  and  downward  prq)agating  waves  would  be  excited.  McPhaden  and  Peters 
(1992)  pointed  out  that  fa:  several  months  each  boreal  spring  the  near-surface  flow  is 
eastward  (McPhaden  and  Taft,  1988),  implying  downward  and  eastward  wave  prt^agation 
during  those  periods.  In  either  case,  the  mean  vertical  velocity  structure  associated  with 
the  EUC  would  surest  downward  and  westward  wave  momentum  flux.  The  atmospheric 
studies  showed  that  the  horizontal  length  scales  of  the  excited  waves  are  several  times 
larger  than  the  height  of  the  convective  boundary  layer,  and  that  the  waves  tend  to  modify 
the  boundary  layer  from  which  they  were  generated  by  some  sot  of  feedback  mechanism 
which  is  not  yet  understood  (Clark  et  al.,  1986). 

Recently,  several  observational  studies  have  supported  the  hypothesized 
role  of  internal  waves  in  the  turbulent  energy  balance  of  die  pycnocline.  Wijesekera  and 
Dillon  (1991)  found  a  pronounced  diurnal  variatiai  in  both  high-frequency  isqiycnal 
displacement  variance  and  turbulent  kinetic  energy  dissipation  rate.  Their  is(^ycnal 
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displacement  variance  values  were  computed  as  the  root  mean  square  (rms)  vertical 
displacement  from  a  Thorpe-reordered  stable  density  proflle  (additional  discussion  on  this 
computation  is  provided  in  Chapter  m)  (Thorpe,  1977).  Hie  vertical  correlation  of 
isopycnal  displacement  variance  indicated  vertical  length  scales  of  the  wave  field  greater 
than  100  m.  They  found  significant  corelation  at  a  0.1-day  lag  between  internal  wave 
displacement  variance  and  turbulence,  as  inferred  from  calculations  of  die  available 
potential  energy  of  fluctuations  (APEF): 

■=  Chi) 

Po  ^ 

Mj.i 

where  d,  p'.  Po,  and  g  are  the  Thorpe  displacement,  the  density  fluctuation  between  stable 
(reordered)  and  unstable  (observed)  density  profiles,  mean  density  (slOOO  kg  m‘^),  and 

f<* 

the  acceleration  of  gravity  (s9.8  m  s'^),  respectively.  APEP  is  a  measure  of  the  maximum 
potential  energy  that  can  be  released  by  a  hnhuient  overturn  (Dillon  and  Park,  1987). 
From  this  0.1-day  lagged  ctxrelation,  they  surmised  that  the  observed  nighttime  mixing 
in  the  pycnocline  could  be  a  result  of  instability  of  the  high-frequen^  internal  waves, 
which  are  initiated  a  few  hours  pritx  to  the-intense  turbulent  mixing.  Using  TROPIC 
HEAT  2  data.  Mourn  et  al.  (1992a)  found  high  correlations  between  turbulent  diss4>ation 
(from  continuous  microstructure  profrles)  and  high-frequency  internal  gravity  waves  (from 
towed  thennistCM’  chain  data)  along  the  equator  from  140^  to  110”W  in  April  1987. 
During  the  first  4  days  of  their  record,  patdies  of  high-frequency  signal  ^>peared 
intermittently  and  were  concentrated  in  the  period  following  sunset  until  several  hours 
past  sunrise.  The  duration  of  the  patches  was  several  minutes  to  many  hours, 
corresponding  to  several  hundreds  of  meters  to  several  tras  of  kilometers  in  the  zonal 
direction.  The  spectrum  of  isopycnal  displacement  variance  was  d<»ninated  by  a  narrow 
band  of  frequencies  between  18-36  q>h  (corresponding  to  1S0-2S0  m  zonal  wavelength) 
of  internal  waves.  Both  turbulence  and  internal  waves  varied  diumally.  By  the  fifth  day, 
the  high-frequency  ccmiponent  of  the  temperature  variance  was  weak,  and  for  the 
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remaining  10  days  there  was  little  or  no  high-frequency  signal  in  the  record.  Mourn  et  al. 
(1992a)  attributed  the  difference  between  the  first  4  days  and  the  remainder  of  the  period 
to  a  drop  in  the  wind  speed. 

McPhaden  and  Peters  (1992)  analyzed  data  from  a  vertical  array  of  fast 
sampling  (1  minute)  temperature  recorders  near  the  surface  on  a  mooring  at  0°,  140^ 
in  May- June  1987.  Their  analyses  indicated  the  existence  of  a  diurnal  cycle  in  IDV  that 
was  most  pronounced  at  frequencies  of  10-30  q>h.  and  that  was  coherent  over  the  31  m 
vertical  extent  of  the  array  between  30  m  and  61  m.  They  estimated  a  ni^ttime  isotherm 
displacement  amplitude  of  0.5-2  m,  and  that  zonal  wavelengths  were  roughly  100-300  m, 
consistent  with  towed  thermistor  chain  data  of  Mourn  et  al.  (1992a).  They  speculated  that 
the  observed  diurnal  cycle  was  due  to  a  combination  of  internal  waves  remotely  generated 
at  the  base  of  the  nighttime  mixed  layer  and  local  shear  instability  of  the  large-scale  flow. 
They  showed  that  the  q)pearance  of  internal  waves  in  the  pycnocline  at  frequencies 
higher  than  the  local  buoyancy  frequency  (4-8  q>h)  could  be  attributed  to  Doppler  shifting 
by  the  EUC.  In  addition  to  their  analysis  of  fast  sampling  temperature  recorders  for  May- 
Jtme  1987,  McPhaden  and  Peters  (1992)  examined  more  coarsely  resolved  IS-minute 
temperature  time  series  for  November  1984,  April  1987,  and  May-June  1987.  When 
averaged  over  30  days,  these  data  also  indicated  the  i»esence  of  a  diurnal  cycle  in  the 
higher  frequencies.  A  comparison  of  the  diurnal  range  in  IDV  for  these  three  periods  led 
them  to  suggest  that  the  magnitude  of  the  diurnal  cycle  may  be  related  to  the  strengdi  of 
the  wind  fmcing  and  intensity  of  the  vertical  shear. 

There  is  evidence  indicating  that  the  November  1984  and  April  1987 
turbulent  dissipation  data  at  0”,  140"W  are  n<M  rqxnsentative  of  mixing  across  the  entire 
equatorial  Pacific,  but  merely  represent  mixing  resulting  frmn  the  conditions  at  that 
location  at  those  times.  In  results  frcmi  the  WEC88  anise  along  1S0*W  in  February  and 
March  1988,  Carr  et  al.  (1992)  noted  an  absence  of  a  clear  diurnal  signal  of  dissipation. 
With  low  wind  speed,  high  incoming  irradiance,  low  sea  surface  tonperature,  and  a 
shallower  EUC,  environmental  conditions  were  notably  different  fixnn  those  during 
TROPIC  HEAT  1.  Similarly,  the  TROPIC  HEAT  2  results  found  diurnal  variability  of 
high-frequency  temperature  signals  during  the  first  4  or  S  days  (nearer  14()‘W)  and  little 
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or  no  high-frequency  signal  fra*  the  remainder  of  the  experiment’s  14  days  as  the  vessel 
moved  eastward  toward  110^  (Mourn  et  al..  1992a).  In  the  data  collected  to  date, 
pronounced  diurnal  variability  of  turbulent  dissipation  in  the  pycnocline  has  been 
observed  only  at  or  near  0°.  140"W.  and  only  during  periods  of  moderate  or  strong 
westward  winds. 

2.  Linear  Wave  Theory 

This  dissertation  focuses  on  the  diurnal  cycle  of  high-frequency  (near  AO 
isotherm  displacement  variance  at  0°.  140°W.  Although  the  wavelike  phenomena  measured 
by  IDV  may  be  non-linear,  the  detailed  structure  is  not  adequately  resolved  to  determine 
the  degree  of  non-linearity,  and  a  linear  analysis  is  undertaken  to  provide  insight  into 
the  behavior  of  small-amplitude  wavelike  fluctuations  in  this  domain.  Preliminary  to  the 
following  statistical  analysis,  this  section  provides  a  brief  review  of  linear  internal  wave 
theory  in  an  equat(»ial  ocean,  as  applicable  to  the  diurnal  cycle  of  the  internal  waves. 

Following  Boyd  (1989),  the  equations  fw  the  vertical  structure  of  linear 
internal  waves  in  mean  shear  flows  are  derived.  Starting  with  the  linearized  set  of 
Boussinesq  equations  in  non-rotating  coordinates  with  mean  flow  U  -  U(z)  in  the  x- 
direction,  the  govorung  equations  are: 


Dm  +  =  'P. 

(3a) 

D,v  =  -p. 

(3b) 

D,w  +  b^  -pj 

(3c) 

o 

II 

(3d) 

I,  +  Vy  +  Wj  a  0, 

(3e) 

in  which  b  =  gppo'  is  the  buoyancy  perturbatitm  with  perturbation  density  p,  and p  is  the 
perturbation  pressure  divided  by  the  mean  density.  The  perturbation  velocities  are  u  = 
(u.v.w)  in  (Cartesian  coordinates,  with  u  eastward.  £>,  =  d,  +  Ud^  incmpo-ates  advection  by 
the  mean  zonal  flow.  Subscr^ts  x,  y,  z,  and  t  denote  differentiation.  Assuming  negligible 
north-south  variation  of  the  mean  fields  and  a  meridional  wavelengdi  mudi  less  than  500 
km,  the  mean  shear  flow  can  be  consitkred  to  be  horizontally  invariant.  For  plane  waves, 
these  simplifications  lead  to  a  vertical  structure  equation  (Boyd  et  al.,  1993): 
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(4) 


C„(z)  +  m*C(z)  =  0. 

where  m  is  the  local  vertical  wavenumber,  given  by 

o*  o 


(5) 


and  a  is  the  intrinsic  frequency,  i.e.  the  frequency  measured  by  an  observer  moving  at  the 
local  mean  speed  U,  and  K  =  (kj,0)  =  (Arcos6./fsin6,0)  is  the  horizontal  wavenumber. 
When  >  0.  the  solutions  G(z)  are  locally  oscillattM^,  and  when  <  0,  they  are 
exponential.  Equations  (4)  and  (5)  plus  boundary  conditions  completely  specify  the  wave 
solutions  to  (3a-e)  when  both  N(z)  and  U(z)  are  specified.  This  dispersion  relation  is 
shown  to  be  singular  at  the  critical  depths,  at  which  the  intrinsic  frequency  vanl£,hes, 
a  =  (o(l-UcosQ/c)  =  0.  Although  it  might  not  be  applicable  to  the  low  Ri  shear  flow  at 
the  equator,  Booker  and  Bretherton  (1967)  have  shown  that  vertically  prq>agating  waves 
are  strongly  attenuated  in  passing  through  the  critical  layer,  provided  Ri  »  0.25. 

Following  Boyd  (1989),  the  polarization  relations  for  the  remaining 
fluctuation  variables  are  determined; 


(6) 


The  polarization  equations  relate  the  fluctuation  variables  u,  v,  w,  and  p  to  the  vertical 
structure  solution  G{z). 

3.  Spectral  Models 

A  kinematic  deaeration  of  oceanic  internal  waves  is  usually  based  on  the 
hypothesis  that  the  wave  field  is  composed  of  a  random  superposition  of  linear  waves. 
The  Garrett  and  Munk  (GM)  spectral  model  has  beemne  a  standard  descr4>tion  for  deq> 
ocean  internal  waves  (Garrett  and  Munk,  1972,1975;  Munk,  1981).  The  GM  model  was 
empirically  derived  from  observations  of  mie-dimensional  spectra  of  fi:equency  frnnn 
moored  current  meters  and  horiztmtal  wavenumber  fr<»n  towed  thennistors.  The  GM 
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formulation  is  based  on  an  infinite  ocean  having  no  top  or  bottom  boundaries  and  a 
slowly  varying  buoyancy  frequency  profile,  N(z),  compared  with  the  fluctuations  of  the 
waves.  This  permits  the  vertical  variation  of  wave  amplitude  and  vertical  wavelength  to 
be  scaled  by  Niz)  using  the  WKB  ^iproximation.  In  the  upper  ocean,  say  frcmi  the  surface 
to  200  m,  the  GM  model  is  expected  to  have  limited  application  due  to  large  changes  in 
buoyancy  frequency  with  depth  and  the  reflection  of  waves  from  the  surface  boundary 
(Levine,  1987).  Despite  these  limitations,  the  GM  spectrum  has  been  compared  with 
upper  ocean  measurements  made  by  acoustic  doppler  profilers,  mooted  current  meters, 
thermistor  chains,  and  CTD  profiles  (Pinkel,  1975,  1981,  1984;  Kase  and  Siedler,  1980; 
Levine  et  al..  1983a,b;  McPhaden.  1985;  Chereskin  et  al.,  1986).  These  comparisons 
indicate  that  while  some  aspects  of  the  GM  spectrum  are  consistent  with  the  data,  many 
are  not. 

Studies  have  suggested  that  space  and  time  variability  in  the  equatorial 
ocean  may  be  substantially  different  from  variability  at  mid  latitudes  (Wunsch  and  Webb, 
1979;  Eriksen,  1980;  Toole  and  Hayes,  1984;  Mcl%aden,  1985).  In  comparing  the  high- 
frequency  variability  of  equatorial  TROPIC  HEAT  microstructure  data  with  tiiat  observed 
at  other  latitudes  in  terms  of  the  GM  model,  Chereskin  et  al.,  (1986)  found  that 
contemporaneous  isopyoud  displacement  and  velocity  spectra  yielded  q)ectral  eneigy 
levels  which  differed  by  a  factor  of  4.  They  found  tiiat  tills  difference  could  be  accounted 
for  by  vertical  advection  by  the  waves  of  the  background  shear  flow.  They  also  found  that 
the  GM  spectral  energy  level  estimated  from  vertical  displacement  was  similar  to  the 
spectral  energy  level  fa*  many  of  the  observations  frmn  whidi  Levine  et  al.  (1985)  had 
removed  the  scaling  of  the  GM  model. 

Boyd  (1989)  and  Boyd  et  al.  (1993)  used  TROPIC  HEAT  observations  and 
a  simple  spectral  model  to  characterize  the  high-frequency  internal  wavefield  in  the 
strongly  sheared  iq^ier  equatorial  Pacific.  Tlwir  observations  ccmsisted  of  temperature  and 
velocity  data  from  four  moorings  spanning  the  equator  between  1.5*N  to  3*S  along  14(rW 
for  the  period  November  1984  to  June  1985.  In  these  p^)ers.  they  compare  tiieir 
observatimis  and  moctel  ouq>ut  to  the  GM  model  spectrum.  The  GM  model  spectrum  is 
based  on  tiie  following  assumptions  which  are  questitniable  in  tiie  central  equatorial 
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Pacific:  i)  the  no-mean-flow  linear  internal  wave  dispersion  relatitMi  is  valid,  ii)  the 
internal  wavefield  is  vertically  symmetric  (same  amount  of  energy  propagates  up  and 
down)  and  horizontally  isotropic,  iii)  the  internal  wave  energy  spectrum  is  8q>arable  into 
a  functi(»i  of  frequency  times  a  function  of  wavenumber,  where  the  function  of  frequency 
is  independent  of  wavenumber  and  the  function  of  wavenumber  is  independent  of 
frequency  except  for  a  bandwidth  scale  factor. 

Following  Boyd  et  al.  (1993),  several  of  the  GM  model  assumptions  are 
expected  to  be  violated  in  the  strongly  sheared  upper  ocean  zonal  mean  flows  at  the 
equator  along  140^,  including  the  following:  i)  die  linear  dispersion  relation  will  be 
modiHed  by  the  vertical  curvature  of  the  mean  flow  and  by  Doppler  shifting  in  the 
intrinsic  frequency;  ii)  the  mean  zcmal  currents  may  impose  directional  asymmetry  on  the 
energy  spectrum  through  critical  layers  involving  zonally  propagating  waves;  and  iii)  in 
the  neighbm’hood  of  the  ocean  surface,  a  fixed  phase  relationship  between  upward  and 
downward  iM*opagating  components  is  expected.  In  addition  to  these  violatimis  of  the  GM 
assumptions,  the  fact  that  linear  internal  waves  advect  the  mean  flow  readily  leads  to  the 
expectation  that  many  of  the  observed  statUtics  will  differ  significantly  from  the  canmiical 
GM  model.  The  fit  to  coherence  of  velocity  components  at  vertical  separations  is  not 
good  (Boyd,  1989).  Vertical  advection,  which  is  not  included  into  the  GM  model,  has 
been  held  responsible  for  the  small  vertical  crdieience  length  scale  and  die  observed 
decrease  in  that  length  scale  with  increasing  frequency  (Muller  et  al.,  1978). 

D.  REVIEW  OF  EQUATORIAL  MODELING  STUDIES 

Several  modelipg  studies  during  the  past  two  decades  have  illustrated  the  role  of 
the  upper  ocean  mixed  layer  in  the  equatorial  circulation  (Hughes,  1980;  Schopf  and 
Cane,  1983;  Garwood  et  al.  1989;  Schudlich  and  Price,  1992).  In  one  of  the  first  models 
designed  to  eiqilain  tropical  and  equatorial  mixed  layers,  Hughes  (1980)  predicted  a 
sinqile  balance  between  local  wind  and  heating,  modified  by  iq>welling,  resulting  in  a 
mixed-layer  depdi  equivalent  to  the  Obukhov  length  scale.  This  simple  1-1/2  layer  model 
did  not  include  an  undercurrent  and  therefore  no  horizontal  advectirm,  and  was  limited 
by  a  simple  parameterizatirm  of  entrainment  into  a  quiescent  lower  layer.  More  realistic 
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primitive  equatimi  models  of  the  laige-scale  equatorial  circulation,  with  Ri-6epeoduA 
mixing,  showed  that  without  a  mean  surface  heating.  tt»  EUC  is  eliminated  by  a 
westward  wind  stress  in  a  manner  of  days  (Pacanowski  ami  Philander,  1981).  This  result 
was  similarly  found  by  Schopf  and  Cane  (1983). 

The  role  of  diurnal  variations  in  surface  heating  was  explored  by  Kraus  (1987). 
He  used  a  mixed-layer  model  that  parameterizes  local  vertical  fluxes  of  energy  and 
momentum  by  the  transilient  method  of  Stull  and  Kraus  (1987).  This  method  involves  a 
spectrum  of  differently  sized  eddies.  Kraus  (1987)  included  a  representation  of  the  mean 
EUC  and  showed  that  turbulent  mixing  occurs  well  below  the  mixed  layer  at  the  equator. 
Garwood  et  al.  (1989)  and  Schudlich  and  Price  (1992)  used  one-dimensimial  models 
initialized  with  realistic  stratificatitm  and  shear  of  the  EUC  to  simulate  the  diurnal  cycle 
of  diss4>ation  observed  during  TROPIC  HEAT.  Garwood  et  al.  (1989),  using  typical 
central  equatoial  Pacific  forcing,  found  that  instabilities  initiated  just  below  the  surface 
layer  ripple  down  into  the  water  column  and  readi  their  deepest  extent  of  about  90  m 
<mly  after  many  hours.  Their  experiments  suggested  diat  Ri  alone  is  a  poor  predictor  of 
the  rate  of  turbulent  mixing  in  the  pycnocline.  Schudlich  and  Price  (1992)  showed  that 
while  the  daytime  phase  iq>pear8  unaffected  by  die  EUC,  die  nighttime  phase  of  die 
diurnal  cycle  is  strcmgly  affected  by  the  EUC.  The  diurnal  cycle  in  their  model  was 
relatively  insensitive  to  other  equatorial  ocean  features,  such  as  the  zonal  pressure  gradient 
and  upwelling.  Although  neitha’  of  these  models  included  internal  gravity  waves  and  both 
simulated  a  diurnal  cycle  in  dissqMticMi,  Schudlich  and  Price  (1992)  showed  diat  the 
models  do  not  simulate  the  observed  intermitteiioe  of  the  dissqiation  on  hourly  time 
scales.  They  suggested  that  diis  might  be  due  to  die  exclusion  of  internal  waves  diat  may 
trigger  local  shear  instabilities. 

Using  a  two-dunensional,  non-hydrostadc,  large-eddy  simulation  model, 
Slcyllingstad  and  Denbo  (1994)  performed  several  eiqieriments  to  investigate  the  formation 
and  maintenanoe  of  high-fiequency  internal  waves  in  the  central  equatorial  Pacific.  In 
their  numerical  experiments,  internal  waves  were  generated  in  the  surfKe  mixed  layer 
by  a  type  of  Kdvin-Helmholtz  instability  diat  is  d^ndent  on  both  die  flow  Reynolds 
number  and  Richardson  number.  It  would  appear  that  Reynolds  number  in  the  central 
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equatorial  Pacific  is  always  greater  than  their  computed  critkal  Reynolds  number,  and 
therefore  not  important  as  a  predictor  of  instability  in  this  region.  Internal  waves 
amplified  whmi  Ri  was  less  than  O.QS  and  decayed  when  Ri  exceeded  O.OS.  Their  model 
results  showed  that  the  diurnal  beating  cycle  drives  a  diurnal  variation  in  the  mixed-layer 
Ri  value,  which  tends  to  modulate  diumally  the  generation  of  the  internal  waves  with  a 
maximum  intensity  just  before  sunrise.  The  modeled  sensitivity  of  internal  wave  formation 
to  the  surface  wind  stress  direction  could  also  be  explained  in  terms  of  reduced  mixed- 
layer  shear.  In  their  simulations,  easttfly  winds  promoted  strong  mixed-layo*  current 
shear,  low  Ri  and  a  strong  diurnal  cycle  of  internal  waves,  while  westerly  winds  led  to 
weak  shear  in  the  mixed  layer,  higher  Ri,  and  reduced  internal  wave  activity. 
Additionally,  Skyllmgstad  and  Denbo  (1994)  found  diat  the  wave  momentum  flux  was 
comparable  to  bulk  estimates  of  large-scale  fluxes  given  in  Hebert  et  al.  (1991)  and  Dillon 
et  al.  (1989),  indicating  that  internal  waves  can  be  a  significant  {vocess  in  die  zmud 
mmnentum  budget  of  the  equatorial  currents.  In  die  surface  mixed  layer,  dw  internal 
waves  gained  eastward  momentum  at  the  expense  of  the  background  flow.  In  some  cases, 
this  momentum  was  transferred  back  to  die  mean  flow  at  a  critical  level  resulting  in  a 
deceleration  below  die  undercurrent  core.  Otherwise,  the  waves  tended  to  decrease  the 
curroit  velocity  above  die  undercurrent  core. 
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Ill.  DATA  AND  METHODS 


A.  DATA  DESCRIPTION 

1.  EPOCS/TOGA  TAO  Mooring  Data 

The  data  used  in  this '  dissertation  consist  pdnuuily  of  time-series 
measurements  from  taut-wire  surface  moorings  located  at  (f,  140*W  in  4300  m  of  water 
for  the  period  from  November  1983  through  October  1987.  Mooring  dq>loyments  and 
recoveries  were  made  at  roughly  six-month  intervals.  These  moorings  were  dq>loyed  as 
part  of  the  National  Oceanic  and  Atmospheric  Administratimi’s  (NOAA)  Equatorial 
Padfrc  Ocean  Climate  Studies  (EPOCS)  Program  and  later  the  Tropical  Ocean-Global 
Atmosphere  (TOGA)  Pr<)gram’s  Tropical  Atmosphere  Ocean  Array  (TOGA  TAO). 

These  mooriqgs  were  instrumented  with  Yellow  Springs  Instrument  air  ami 
sea  surface  temperature  sensors  at  nominal  hei^ts/depdis  of  3  m  above  die  sea  surface 
and  1  m  below  die  sea  surface,  reflectively,  seven  EGftG  model  610  Vector  Averaging 
Cunent  Meters  (VACMs)  and  six  Sea  Data  Temperature  Recorders  (TRs)  in  the  upper 
300  m,  and  Vector  Averaging  Wind  Recorders  (VAWRs)  at  a  height  of  4  m  above  die 
mean  water  line.  Several  EG&G  Model  6^  Vector  l^asuring  Current  Meters  (VMQ^) 
were  also  used.  The  VACMs  recorded  IS-minutejiveraged  velocities  and  temperatures 
and  the  TRs  recorded  IS-minute  fiot  samples  of  temperature.  The  VAWR  is  an  inverted 
VACM  equipped  with  a  Qimet  diree-cup  anemometer  and  a  pivoted  vane  (Hreitag  et  al., 
1989).  Thoe  have  been  some  changes  in  subsurface  instrument  dqHhs  during  the  4  year 
period  of  this  investigation  due  to  changing  {X’pgram  priorities.  Hgure  3.1  illustrates  a 
typical  mooring  arrangement  Tables  3.1  and  3.2  provide  actual  time  lines  for  instruments 
at  each  of  die  dqiths  used  in  diis  analysis  for  temperature  and  velocity,  respecdvdy. 
There  ve  gfis  in  some  of  the  time  series  due  to  occasional  instrument  failures.  Fuidier 
informadon  <m  the  coUecdon  and  irudal  {vooessing  of  these  data  is  available  in  Freittig 
et  aL  (1987),  Fieitag  and  McPhadra  (1988),  and  Yue  et  al.  (1991).  Moodily  means  and 
seasonal  cycles  of  each  of  these  time  series  measurements  for  the  period  1983-1991  are 
available  in  h^Phaden  and  McCarty  (1992). 


22 


1983  -October  1987 


23 


- CURRENT/WIND 

0“,  140“W  - TEMPERATURE 

0 

50 

£  100 
S  150 

CL 

^  200 
250 
300 

1983  1984  1985  1986  1987 


Table  3.1.  Bar  chart  showing  velocity  and  temperature  data  available  from  current  meters 
and  from  wind  recorders  as  a  function  of  dq>diiU0*,  140n^  during  the  period  1983-1987. 
Arrows  indicate  times  of  mooring  dqjloyments  and  recoveries  and  times  when  CTD  casts 
were  conducted. 
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TaMe  32.  Bar  chart  showing  availability  of  data  from  temperature  recorders  (TRs)  at 
0“,  140W  during  the  period  1983-1987. 
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In  addition  to  these  standard  equatorial  mooring  measurements,  two  of  the 
moorings  were  instrumented  with  a  vertical  array  of  fast-sampling  1 -minute  TRs  designed 
to  resolve  high-frequency  internal  wave  variations.  The  duration  of  these  data  are  typically 
days  to  weeks.  The  first,  in  May-June  1987,  had  fast-sampling  TRs  at  30  m,  34  m,  36  m, 
40  m,  46  m,  59  m,  61  m,  and  101  m.  Processing  of  these  data  is  summarized  in  I^reitag 
and  McPhaden  (1988).  Using  this  data,  McPhaden  and  Peters  (1992)  rqHirted  a  diurnal 
cycle  in  IDV  that  was  most  pronounced  at  frequencies  of  10-30  q>h  that  was  coherent 
between  30  m  and  61  m.  The  second  mooring  was  equipped  with  fast-sampling  TRs  at 
17  m,  24  m,  30  m,  36  m,  47  m.  59  m,  70  m.  82  m,  101  m,  121  m,  160  m,  199  m,  and 
250  m,  an  acoustic  Doppler  current  profiler  (ADCP)  with  8  m  vertical  resolution  between 
10-250  m,  and  a  shtxtwave  radiometer.  This  unique  30  day  data  set  was  obtained  in 
conjunction  with  microstructure  measurements  of  the  Trq>ical  Instability  Wave 
Experiment  (TIWE)  conducted  in  Novembtf/December  1991.  Although  this  dissertation 
is  not  directly  concerned  with  the  analysis  of  diis  sec(»d  fast-sampling  data  set,  some  of 
the  results  of  that  analysis  are  discussed  here  to  assist  in  the  analysis  of  the  15-minute 
time  series  data. 

2.  Climatological  Data 

In  addititm  to  using  the  moored  buoy  data  to  characterize  the  high- 
frequency  temperature  variability,  several  climatological  data  sets  are  utilized  to  provide 
additional  information  about  surface  forcing  parameters.  The  Intonational  Satellite  Cloud 
Climatology  Project  (ISCCP)  data  were  used  to  compute  monthly  mean  surface  irradiance 
values  0°,  140^  for  the  4-year  obsavadon  period.  ISCCP  data  were  also  used  to 
compute  monthly  mean  estimates  of  surface  out-going  longwave  radiation.  The  ISCCP 
data  contains  information  on  cloud  cover,  cloud  type,  and  cloud  heights,  as  well  as 
information  about  some  of  the  atmospheric  composition  and  surface  conditions  (Rossow 
et  al.,  1988;  Rossow  and  Schiffer,  1991).  These  data  are  resolved  globally  at  2.5"  x  2.5" 
resolution  every  three  hours.  The  Comprehensive  Ocean-Atmosphere  Data  Set  (COADS), 
Esbensen  and  Kushnir  (1981),  Oberhuber  (1988),  Dorman  and  Bourke  (1976),  and  Global 
Precipitation  Climatology  Project  (GPCP)  climatologies  were  also  used  to  ctnnpute 
surfaces  fltixes. 
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B.  DATA  PROCESSING 


Preliminary  post-processing  of  the  raw  mooring  data  was  performed  at  NOAA’s 
Pacific  Marine  Environmental  Laboratory  following  Freitag  et  al.  (1991).  This  included 
converting  the  data  to  engineering  units  using  empirically-derived  calibration  coefficients 
and  performing  internal  quality  checks  and  windowing  to  remove  obviously  bad  data 
points.  Since  this  investigation  focuses  on  high-frequency  variations,  the  resulting  15- 
minute  data  were  visually  examined  to  further  verify  the  quality.  Questionable  data  were 
eliminated  and  flagged. 

With  the  objective  of  characterizing  the  diurnal  cycle  of  high-frequency 
temperature  variability  over  seasonal  and  interannual  time  scales,  subsurface  temperature 
and  velocity  data  over  the  4-year  period  from  November  1983  through  October  1987  were 
examined.  To  investigate  the  variability  of  the  diurnal  cycle  over  these  longer  time 
scales,  the  data  were  segmented  by  calendar  month  and  then  into  six-hour  time  bins.  The 
local  time  bins,  2200-0400  (night),  0400-1000  (mwning),  1000-1600  (day),  and  1600- 
2200  (evening),  were  selected  to  correspond  with  observations  of  the  diurnal  cycle  of 
turbulent  dissipation,  and  to  the  obso^ations  by  McPhaden  and  Peters  (1992)  of  the 
diiunal  cycle  of  high-frequency  internal  waves.  The  TROPIC  HEAT  observations  showed 
maximum  rates  of  turbulent  dissipation  and  high-frequency  internal  wave  activity  in  the 
stratified  region  during  the  mtxning  hours  near  sunrise,  conresponding  to  the  0400-1000 
time  bin  in  this  study.  Minimum  values  of  dissipation  typically  occurred  during  the 
evening,  or  the  1600-2200  time  bin.  In  discussing  diurnal  variations,  the  first  two  time 
bins  (night  and  morning)  will  be  collectively  refnred  to  as  nighttime  and  the  runaining 
two  time  bins  (day  and  evening)  will  be  collectively  referred  to  as  daytime. 

In  addition  to  allowing  a  comparison  of  different  six-hour  bins,  and  thus  an 
indication  of  the  mean  diurnal  cycle,  segmenting  the  data  into  six-hour  segments 
effectively  acts  as  a  high-pass  filter  to  reduce  the  effects  of  the  semi-diurnal  and  diurnal 
tides.  Additional  measures  taken  in  computing  IDV  to  filter  the  tides  are  discussed  in 
greater  detail  in  the  next  section. 

In  examining  statistical  relationsh4>s  throughout  this  analysis,  the  standard  error 
of  the  mean  and  95%  confidence  intervals  are  often  used  as  indicattx^  of  the  level  of 
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statistical  significance.  The  standard  error  of  the  mean  was  computed  as  — H— 

where  a  is  the  standard  deviation  of  the  peculation  of  n  observations  (Crowe  et  al.. 

1960).  The  95%  confidence  intervals  were  computed  based  on  the  integral  times  scale 

between  of  the  respective  time  series  using  the  method  of  Davis  (1976).  The  integral  time 
«• 

scale;  =  E  Cff,iAt)C^(,iAt)At  .  determines  the  time  period  required  to  gain  a  new 

"degree  of  freedom"  in  computing  correlations.  The  use  of  the  integral  time  scale  in 
determining  significance  levels  is  particularly  important  due  to  the  strtmg  seasonal  cycles 
of  the  surface  forcing  and  background  conditions  obso^ed  in  the  equatorial  Pacific 
Ocean. 

1.  Isotherm  Displacement  Variance 

As  an  indicator  of  the  seasonal  and  interannual  variability  of  the  high- 
frequency  temperature  variability,  monthly  mean  values  of  IDV  were  computed  for  the 
48  months  of  the  record.  IDV  was  computed  at  each  of  six  depths  instrumented  with  15- 
minute  spot-sampled  temperature  recorders  CTRs);  35  m,  60  m,  100  m,  140  m,  200  m,  and 
300  m.  Although  temperatures  were  recorded  at  15-minute  into^als  at  the  each  of  the 
VACM  depths,  1  m,  10  m,  25  m,  45  m,  80  m,  120  m,  160  m,  250  m,  these  temperature 
values  represented  15-minute  averaged  temperature,  rather  than  instantaneous  temperature. 
The  15-minute  averaging  operator  effectively  removed  the  higher-frequency  waves  of 
interest  here.  It  is  important  that  even  the  15-minute  spot-sampled  TR  data,  with  a 
Nyquist  frequency  of  2  eph,  is  too  coarse  to  resolve  the  high-fiequency  internal  waves 
(10-30  q}h)  observed  by  Mcl%aden  and  Peters  (1992)  using  1-minute  spot-sampled  TR 
data.  Nevertheless,  with  sufficient  averaging,  over  about  30  days,  MePhaden  and  Peters 
(1992)  rqxxled  a  diurnal  cycle  of  IDV  using  15-minute  TR  data. 

Prior  to  computing  variances  fa:  each  of  the  six-hour  rectvds,  each  six-hour 
segment  was  fitted  to  a  cubic  polynomial  to  remove  the  mean,  trend  and  curvature  of  the 
temperature  signal.  The  purpose  of  this  qperatiai  was  to  remove  the  lower  fiequencies, 
particularly  the  tidal  fi'equencies.  Figure  3.2  illustrates  the  simple  method  used  to  remove 
the  denninant  tidal  effects.  The  advantages  of  this  method  are  that  it  is  simple  and  local 
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because  it  was  applied  only  to  six-hour  segments;  hence,  non-stationarity  of  the  data  does 
not  present  a  serious  problem.  Similar  analyses  using  a  32-term  tidal  harmonic  analysis 
routine  and  an  ideal  (cosine  taper)  high-pass  filter  to  the  monthly  data  sets  resulted  in 
undesirable  effects  that  were  attributed  to  the  non-stationarity  of  the  data  over  monthly 
time  scales.  An  examination  of  the  four-year  records  of  raw  temperature  and  velocity  data 
clearly  shows  the  non-stadonary  character  of  die  data. 

After  fitting  with  cubic  polynomials,  variances  of  the  residuals  were  computed  for 

each  six-hour  segment  fm*  each  day  of  each  month  fcx  aU  six  TR  depths.  IDV  was  then 

computed  by  normalizing  these  variances  by  the  mean  temperature  gradient  for  the  same 

six-hour  pmod.  The  mean  temperature  gradient,  dT/dz,  was  computed  from  the  six-hour 

mean  temperatures  at  each  depth  using  centered  differencing,  except  at  the  surface  and 

at  300  m,  where  first  differences  were  used.  A  sensitivity  study  was  perfrumed  to  evaluate 

the  difrerences  between  normalizing  by  the  mean  dT/dz  fm-  the  six-hour  segment  averaged 

over  the  typical  30-day  month  (which  was  more  stable,  but  subject  to  nmi-stationarity 

errors)  or  normalizing  by  the  mean  dT/dz  fw  the  individual  six-hour  segments,  subject 

to  noise  due  to  small  instrumental  errors.  The  results  showed  that  the  possible  error  due 

to  instrument  error  was  much  less  than  that  due  to  non-stationarity  of  the  gradient  over 

30  days  (Figure  3.3).  Thus,  if  a  computed  gradient  was  less  than  the  minimum  resolvable 

gradient  computed  assuming  a  maximum  instrument  errm*  of  0.05  C,  then  that  observation 

was  not  used  in  computing  IDV.  The  minimiun  resolvable  gradient,  graderr  was 

computed  as  graderr  ^  ,  where  Az  is  the  instrument  spacing.  To  greatly 

Az 

improve  the  high-frequency  signal-to-noise  ratio,  the  six-hourly  IDV’s  fcnr  each  of  the  four 
time  bins  w^-e  ensemble-averaged  over  all  days  of  each  month. 

Since  IDV  was  computed  by  dividing  the  six-hour  temperature  variance  by 
the  six-hour  mean  vertical  temperature  gradient,  this  calculation  was  extremely  smsidve 
to  small  temperature  gradients,  such  as  those  occurring  in  the  surface  mixed-layer. 
Tho’efix’e,  it  was  impotant  to  distinguish  between  periods  when  the  instruments  were  in 
the  mixed  layer  and  periods  when  the  instruments  were  below  the  mixed  layo*.  When  the 
35  m  and  60  m  TRs  were  in  the  mixed-layer,  gradients  were  small  and  IDV  would  be 
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Figure  32.  Dlustradon  of  procedure  used  to  hi^-pass  filter  the  temperature  data  to 
reduce  principal  tidal  constituents.  Data  woe  segmented  into  six-hour  time  bins  and  fit 
with  a  cubic  polynomial.  JUustration  shows  that  die  polynomial  removes  the  local  mean, 
six-hour  trend,  and  die  curvature  and  inflecdcm  point  associated  with  the  semi-diurnal 
tide. 
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-  STANDARD  DEVIATION 

- INSTRUMENT  ERROR 


dT/dz  Ettots  (IC*  "Cm'*) 

Temperatuie  Gradient  Error  at  0®,14(W  (Nov  84) 


Figure  3J.  Temperature  gradient  error  analysis  used  to  determine  the  preened 
temperature  gradient  (dT/dz),  which  was  used  to  normalize  temperature  variances  in 
computing  isotherm  di^lacement  variance  (IDV).  The  dashed  line  rqxesents  die 
maximum  error  in  dT/dz  associated  widi  the  instrument  errors.  Maximum  instrument 
errms  were  assumed  to  be  O.QS  C.  The  four  solid  lines  are  standard  deviations  of  dT/dz 
over  the  month  of  November  1984.  These  results  show  that  the  non-stadonarity  of  dT/dz 
over  the  month  is  greater  than  the  possible  error  in  dT/dz  due  to  instrument  eirw. 
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expected  to  become  large  and  meaningless.  Eliminating  six-hourly  segments  having 
vertical  temperature  gradients  less  than  0.004  Cm  '  at  35  m  and  0.002  Cm  '  at  60  m 
served  to  filter  out  time  periods  when  the  TRs  were  in  the  mixed-layer.  In  addition  to  the 
temperature  gradient  tests,  six-hourly  values  of  IDV  were  subjectively  edited  to  ronove 
unrealistically  high  IDV  values  before  computing  monthly  means.  Over  the  four-year 
record,  six-hourly  values  of  IDV  at  35  m  were  eliminated  by  these  two  tests  14.7%. 
i5.1%.  15.1%.  and  9.9%  of  the  time  for  the  night,  morning,  day.  and  evening  time  bins, 
respectively.  Six-hourly  values  of  IDV  at  60  m  were  eliminated  4.1%.  2.9%.  2.6%.  and 
2.6%  of  the  time  for  the  night,  morning,  day.  and  evening  time  bins,  respectively. 

2.  Buoyancy  Frequency 

The  static  and  dynamic  stability  of  the  upper  ocean,  as  well  as  the 
corresponding  generation  and  decay  of  internal  waves  and  turbulent  mixing  are 
fundamentally  related  to  the  density  stratification  of  the  water  column.  Therefore, 
estimates  of  the  buoyancy  frequency. 

p«  as 


were  computed  for  the  48-month  analysis  period.  The  lack  of  a  moored  salinity  time 
series  at  0".14(rw  precluded  the  inclusion  of  high-frequency  salinity  fluctuatims  in  the 
computation  of  N,  Cmiductivity.  temperature,  and  dq)th  casts  (CTDs)  were  conducted 
at  about  six-mmith  intervals  during  each  mooring  dq>loyment  and  recovery  (Table  3.1). 
Examination  of  the  salinity  profiles  from  these  CTTD  casts  suggests  that  although  salinity 
does  change  with  time,  particularly  during  the  migratim  of  a  salinity  fremt  during  the 
passage  of  tropical  instability  waves,  the  vertical  gradient  of  salinity  ai^)ears  to  be 
relatively  constant  compared  with  changes  in  vertical  tenq>erature  gradients.  Using  the 
available  salinity  prt^es.  estimates  of  monthly  salinity  profiles  were  ctmiputed  by  linear 
interpolation.  Actual  monthly  salinities  may  have  varied  substantially  from  tiiese 
estimates.  In  additimi  to  these  monthly  salinity  profrle  estimates,  six-hourly  mean 
temperature  profiles  were  computed  using  both  VACM  and  TR  temperatures  at  14  dq>tiis. 
Using  the  monthly  salinity  values  and  tiie  six-hourly  mean  temperature  values,  six-hourly 
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mean  densities  were  computed  at  each  of  the  14  depths.  These  density  profiles  were  then 
fit  with  a  cubic  spline  (Akima,  1970)  designed  so  that  the  shape  of  the  curve  matches 
the  shape  of  the  data  while  minimizing  oscillations.  The  vertical  density  gradient,  dp^. 
was  computed  from  the  spline  fit  of  the  density  profile.  In  this  manner.  N  was  computed 
at  each  of  the  14  instrument  depths  for  each  six-hour  period. 

3.  Velocity  Shear 

Dynamic  stability  is  determined  in  part  by  the  velocity  structure  of  the 
water  column,  particularly  the  vertical  velocity  shear.  Therefme,  six-hourly  values  of 
vertical  velocity  shear  were  computed  from  six-hourly  mean  velocity  profiles  computed 
from  the  IS-minute  averaged  VACM  velocity  data  at  10  m,  25  m,  45  m,  80  m,  120  m. 
160  m,  and  250  m.  Six-hour  mean  zonal  and  meridional  velocity  components,  u  and  v, 
were  fit  with  a  cubic  spline  (Akima,  1970),  as  described  in  the  previous  section.  Vertical 
shear  of  the  zonal  and  meridional  components,  Birdz  and  dvdz,  respectively,  were 
computed  from  the  spline  fit.  In  a  comparison  study  using  8-meter  vertical  resolution 
acoustic  Doppler  current  profiler  (ADCZP)  data  frmn  1991,  the  Akima  spline  [xrovided  a 
more  realistic  shear  profile  than  that  computed  by  linear  interpolation  from  centered 
difierencing  (Figure  3.4).  Shear  squared  values  in  the  high-shear  region  below  the  EUC 
axe  are  underestimated  due  to  the  greater  vertical  sq;>aration  of  current  meters  below  120 
m.  Although  these  shear  estimates  are  far  from  ideal,  they  are  reastmable  and  likely  die 
best  available  fix  the  time  and  location  of  this  study.  Vector  shear  and  shear  squared  were 
computed  from  the  comptxient  shears. 

4.  Rkhardson  Number 

Gradient  Richardson  number,  Ri,  an  indicator  of  the  dynamic  stability  of 
the  water  column,  is  defined  as  the  ratio  between  the  density  stradficaticm,  N’,  and  the 
velocity  shear  squared,  5’.  Generally,  the  small-scale,  or  local,  instantaneous  Ri  is  desired, 
but  is  often  not  available  due  to  limited  vertical  and  temporal  resolution.  Clearly,  the 
resolution  of  these  time  series  measurements  do  not  satisfy  the  local  and  instantaneous 
conditions.  Nevertheless,  for  the  objectives  defined  here,  i.e.  to  characterize  the  seasonal 
and  interannual  variability  of  the  diurnal  cycle  oS  high-frequency  temperature 
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variability,  or  internal  waves,  mean  Ri  provides  a  useful  indicator  of  mean  processes 
occurring  on  these  longer  time  scales. 

Six-hourly  averaged  Ri  was  computed  from  the  ratio  of  die  six-hourly 
to  six-hourly  S^.  This  ratio  is  likely  to  difrer  from  the  mean  of  the  ratio  of  IS-minute 
to  the  IS-minute  S^.  Monthly  mean  values  of  Ri  were  cmnputed  for  each  of  the  four  six- 
hour  time  bins.  Since  is  sometimes  vanishipgly  small,  some  values  of  Ri,  which  has 
in  the  denominator,  becmne  quite  large.  To  prevent  these  occasional  large  values  from 
significantly  weighting  the  monthly  means,  values  exceeding  ISO  were  truncated  and  the 
remaining  values  were  smoothed  with  a  1-2-1  running  mean  filter.  Due  to  the  averaging 
operators  used  and  the  coarse  votical  resolution  of  the  data,  particularly  velocity,  die 
munerical  values  of  mean  Ri  are  expected  to  be  somewhat  higher  dian  those  measured  on 
smaller  scales  using  microstructure  measurements. 
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EVALUATION  OF  CUBIC  SPLINE  FOR  S*  &  U 


Velocity  Shear  Squared  (■'*) 


F^re  3.4.  Compaiistm  of  vertical  profiles  of  a)  zonal  velocity  and  b)  velocity  shear 
squared  computed  £r(»n  8  m  vertical  resolution  ADCP  data  at  OP,  140^  (sdid)  and 
cubic  q)line  fit  to  7  ADCP  dq)ths  cortespmiding  to  cuirent  meter  dqjths  (dashed). 
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5.  Mix«<l-L4iyer  Depth 

The  thickness  of  the  mixed  layer  determines  both  the  thermal  and  the 
mechanical  inertia  of  die  layer  in  direct  contact  with  the  atmoqihere  and  influences  the 
fluxes  at  the  base  of  the  layer.  The  study  of  mixed-layer  dqith  and  its  variability  is  of 
fundamental  importance  for  the  understanding  and  intenxetation  of  diennal  and  velocity 
fields  of  the  upper  ocean,  for  parsmeterizing  mixed-layer  processes,  and  for  determining 
air-sea  fluxes  (Schneider  and  Muller.  1990).  Numerous  defmitions  for  mixed-layer  depth 
have  been  described  in  the  literature.  Gradients  or  differences  frmn  the  surface  of  either 
temperature,  salinity,  or  density  have  been  specified  that  cannot  be  exceeded  within  the 
mixed  layer.  Wyrtki  (1964)  used  a  temperature  difference  of  0.5  C.  Lukas  and  Lindstrom 
(1987)  employed  critical  gradients  of  temperature  (O.OS  Cm  '.  0.025  Cm*'),  salinity  (0.02 
m  '.  0.01  m'').and  potential  density  Og  (0.01  kgm'^).  Peters  et  al.  (1989)  used  a  difference 
of  Og  of  0.01  kgm*’.  McPhaden  ai^  Peters  (1992)  used  a  temperature  gradient  of 
O.OlCm*'.  Ntme  of  these  definitions  may  be  ideal  for  every  situation,  and  it  is  advisable 
to  examine  individual  data  sets  {xior  to  defining  the  requisite  criteritm.  After  examining 
individual  and  monthly  mean  profiles  of  temperature  and  temperature  gradient,  the 
definition  of  mixed-layer  depth  based  <m  the  depth  when  the  gradient  first  reaches 
O.OlCm*'  was  chosen.  Using  this  criterimi.  mixed-layer  depth  was  computed  for  each  six- 
hour  time  bin  by  linearly  interpolating  the  six-hour  mean  temperature  jxofiles.  Monthly 
mean  mixed-layo-  dq>ths  were  then  calculated  for  each  of  the  four  time  bins. 

As  an  indicator  of  the  diurnal  variability  of  surface  layer  mixing,  daily 
maximum  ranges  woe  computed  fixim  the  four  daily  six-hour  values  of  mixed-layer 
dq)th.  i.e.  the  difference  betwemi  each  daily  maximum  mixed-layer  dqrth  and  each  daily 
minimum  mixed-layer  dq)di.  Mtmthly  mean  diurnal  ranges  of  mixed-layer  dq>th  were 
computed  from  the  daily  values. 

In  this  study,  both  mirnd-layer  depth  and  diurnal  range  of  mixed-layer 
depth  are  used  as  measures  of  the  net  amount  of  surface  layer  mixing.  Since  high- 
frequency  temperature  fluctuations  in  the  thermocline  are  likely  to  be  related  to  surface 
layer  mixing  and  the  processes  which  generate  surface  layer  mixing,  relationships  between 
these  measures  of  surface  layer  mixing  and  IDV  and  die  diurnal  cycle  of  IDV  are 


36 


examined. 


6.  Surface  Heat  and  Buoyancy  Fluxes 

Diurnal  modulations  of  the  upper  ocean  mixing  processes  observed  at 
0",140^  are  driven  fr(»n  the  surface  by  either  mechanical  forcing  of  surface  wind  stress 
or  thermodynamic  fcxcing  of  surface  buoyancy  flux.  Diurnal  variations  in  mechanical 
forcing  are  possible  since  there  are  small  but  significant  diurnal  variaticHis  in  wind  speed 
(Deser,  1994).  These  diurnal  variations  are  likely  to  be  smaller  than  the  pronounced 
diurnal  variation  of  the  surface  buoyancy  flux.  The  surface  buoyancy  flux  determines  the 
amount  and  intensity  of  convective  mixing.  The  ocean  moves  because  of  buoyancy 
contrasts  due  to  differences  of  both  temperature  and  salinity.  These  differences  are  created 
by  the  fluxes  of  heat  and  water  at  the  ocean  surface,  whose  combined  effect  on  buoyancy. 
-gp,  is  called  the  buoyancy  flux,  B,  which  is  given  by: 

where  E  =  ^  evaporation  rate,  P  is  the  precipitation  rate,  is  the 

specific  heat  of  water  (3998  Jkg'^K*'),  5  is  the  surface  salinity  (35.0),  a  »  -p'^dpldT 
is  the  thermal  expansion  coefficient  of  ^water  at  the  surface  (2970  x  10-7  K‘‘). 

p  =  p~'dp/as  is  the  c(»responding  coefficient  for  salinity  (-7.39  x  10^  ),  g  is 
the  acceleration  of  gravity  (9.7976  ms'^),  {/  is  the  surface  wind  ^>eed,  q„  is  the  specific 
hiimidity,  is  the  saturated  specific  humidity,  and  Q,^  is  the  net  surface  heat  flux  (Gill, 
1982).  Buoyancy  flux  is  a  measiire  of  the  input  of  mean  potential  energy  at  the  ocean 
surface  (Oberhuber,  1988). 

Surface  heat  flux  is  governed  by  the  balance  between  shortwave  solar 
irradiance,  (2.,  latent  or  evaporative  flux,  Q^,  sensible  or  conductive  flux,  and 
longwave  back  radiation, 

Q«  =  (2,  +  a  +  C2b  +  eh 

Monthly  mean  values  of  and  B  were  computed  from  monthly  mean  values  of  Q,, 
l2e>  l2h>  methods  and  data  sets  used  in  computing  these  terms  are  givra 
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in  the  following  sev^al  sections. 

Since  it  has  been  hypothesized  that  the  diunial  cycle  of  the  surface 
buoyancy  flux  is  related  to  the  diurnal  cycle  of  turbulent  mixing  in  the  equatorial 
thermodine  (Mourn  and  Caldwell.  198S;  Mourn  et  al..  1989),  surface  buoyancy  and  heat 
fluxes,  as  well  as  the  individual  terms  comprising  them,  are  examined  over  seasonal  and 
interannual  time  scales.  Relationships  between  these  terms  and  surface  layer  mixing,  EDV 
and  the  diurnal  cycle  of  IDV  are  determined. 

Although  there  is  evidence  of  a  diurnal  cycle  of  predpitation  over  the 
tropical  oceans  (Gray  and  Jacobson.  1977;  Dorman  and  Bourke,  1979;  Meisner  aiui  Arkin, 
1987),  observatirMis  of  the  diurnal  cycle  of  turbulent  mixing  in  the  equatorial  Pacific  were 
observed  during  times  when  predpitation  was  not  observed  and  the  surface  buoyancy  flux 
was  controlled  by  the  stirface  heat  flux  (Peters  et  al.,  1988;  Mourn  et  al.,  1989). 
Therefore,  it  is  not  unreasonable  to  assume  that  the  surface  heat  flux,  which  changes 
diumally  with  the  rising  and  setting  of  the  sun,  dominates  the  diurnal  cyde  of  the  surface 
buoyancy  flux  at  0”,  140^. 

7.  Surface  Solar  Irradiance 

Bond  and  McPhaden  (1994)  and  Imawaki  et  al.  (1988)  found  dut  the  day- 
to-night  difference  of  the  sea  surface  cooling,  which  they  define  as  the  sum  of  the  net  Q^, 
(2|„  and  Q^,  is  negligible  con^Mued  with  the  net  solar  irradiance.  As  e]q>ected,  they  found 
large  day-to-night  differences  in  (2,.  Additionally,  Bishop  and  Rossow  (1991)  described 
significant  variability  of  Q,  over  time  scales  varying  from  daily  to  interannual.  They 
attribute  this  variability  to  changes  in  doud  fraction  and  atmospheric  composititm. 

Pa-  the  four-year  period  examined  here,  direct  observations  of  solar 
irradiance  were  not  available  in  the  central  equatorial  Pacific.  However,  mmithly  mean 
estimates  of  surface  solar  irndiance  have  been  computed  for  the  entire  time  series  using 
ISCCP  cloud  data.  The  computatitmally  fast  radiative  transfer  modd  of  Bishop  and 
Rossow  (1991)  for  surface  solar  irradiance  utilizes  the  following  ISCCP  data:  solar  zenitii 
angle  (sj,  atmospheric  water  vapor  profile  and  ozone  column  abundance,  doud  fractitm, 
doud  (^tical  thickness,  visible  surface  reflectance,  and  surface  pressure.  This  algorithm 
uses  the  Frouin  et  al.  (1989)  formula  for  the  solar  irradiance  at  the  surface  under  dear  sky 
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conditions.  The  cloudy  sky  ccnnponent  utilizes  spherical  and  directi<Mud  albedos  derived 
using  optical  thickness  frmn  tables.  These  tables,  which  allow  a  major  improvement  in 
computation  speed,  are  obtained  from  full  Mie  calculations  (Hansen  and  Travis,  1974)  and 
provide  a  physical  link  between  the  observed  cloudy  scene  reflectance  and  the  total  cloud 
transmission  that  is  the  basis  for  all  methods  for  calculating  surface  solar  irradiances 
(Taipley,  1979;  Gautier  et  al.,  1980;  Darnell  et  al.,  1988;  Bish(^  and  Rossow,  1991).  The 
total  incident  solar  irradiance  is  the  sum  of  die  clear  and  cloudy  sky  cmnponents. 

Monthly  mean  surface  solar  irradiance  values  for  the  period  November 
1983  through  October  1987  were  provided  by  Dr.  J.  Beauchamp  at  Goddard  Institute  for 
Space  Studies.  These  monthly  means  were  cttnputed  from  3-hourly  2.5”  x  2.5°  ISCCP 
data.  In  addition  to  the  surface  solar  irradiance  values  computed  using  the  radiative 
transfer  model,  monthly  mean  values  of  surfa^  solar  irradiance  were  computed  from  bulk 
formalae,  QJQt  =  1  -  0.76C  +  0.00190,  where  G,  is  the  insolati<m  under  cloudy 
conditions,  Qo  the  insolation  under  clear  skies  from  the  Smithsonian  formula.  C  cloud 
fraction  from  ISCCP  data,  and  o  is  nocm  solar  altitude  (Reed.  1977;  Weare  et  aL,  1981; 
Seager  et  al.,  1988;  Blumenthal  and  Cane,  1989). 

8.  Longwave  Radiation 

Although  the  diurnal  cycle  of  Itmgwave  radiatimi  has  been  reported  to  be 
small  (Imawaki  et  al.,  1988;  Bond  and  Mcl%aden,  1994)  the  concern  here  is  with  die 
effect  of  seasonal  and  interannual  variations  of  Imigwave  radiation  on  the  diurnal  cycle 
of  upper  ocean  mixing  and  high-frequency  temperature  variability.  During  this  four-year 
study  period,  however,  no  direct  measurements  of  net  longwave  radiation  were  made.  In 
fact,  due  to  the  difficulty  of  making  direct  radiatitm  measurements  at  sea,  there  are  few 
long-term  observations  of  longwave  radiatitm  over  the  oceans  at  any  location.  Data 
acquired  during  the  GARP  Atlantic  Tropical  Eiqieriment  (GATE)  revealed  spatial  and 
tempm-al  variability  of  longwave  radiation  of  about  30  Wm presumably  due  to 
mesoscale  and  small-scale  variations  in  cloudiness  and  other  meteorological  parameters 
(Fung  et  al.,  1984).  Surface  longwave  radiation  depends  on  the  integrated  results  of 
absorption,  emission,  and  scattering  in  the  mitire  inhomogeneous  atmospheric  column. 
Due  to  the  scarcity  of  direct  observations,  empirically  dmved  bulk  formulae  are 


39 


traditionally  used  to  estimate  monthly  mean  net  longwave  radiation  at  the  ocean  surface 
(Clark  et  al.,  1974;  Bunker,  1976;  Wyrtki,  1965;  Hastenrath  and  Lamb,  1978).  Although 
these  methods  do  not  adequately  account  fci'  the  dependence  on  the  concentrations  of 
atmospheric  constituents  and  on  the  vertical  distribution  of  temperature,  humidity,  and 
cloud  prcq)erties,  they  are  generally  the  best  available  estimates.  The  various  bulk 
formulae,  which  use  marine  observations  of  SST,  near-surface  air  temperature,  humidity, 
and  cloudiness,  often  give  dissin.ilar  results  for  given  parameters  (Fung  et  al.,  1984). 

Monthly  mean  net  longwave  radiation  was  computed  following  the 
procedure  given  in  the  Ocean  Observing  System  Development  Panel’s  Background  Report 
on  Surface  Conditions  and  Air-Sea  Fluxes  (Weller  and  Taylor,  1993),  which  uses  the  clear 
sky  ftxtnula  of  Clark  ei. (1974), 

LW^U  =  €o7^(0.39-0.05(ej"®)  ^  4€o7^(r,-r^  . 

The  emissivity,  e,  is  taken  to  be  0.98;  a  -  5.67  x  10^  Wm'^K"*  is  the  Stefan-Boltzmann 
constant,  and  near-surfaoe  air  temperature  (TJ  and  SST  (TJ  are  in  degrees  Kelvin.  Input 
data  were  gathered  from  the  following  sources  in  tvder  of  preference:  mooring  data  for 
T,  and  T„  when  available;  COADS  monthly  means  of  T„  T„  specific  humidity  (q),  and 
sea  level  pressure  (p),  when  available;  and  the  Esbensen  and  Kushnir  (1981)  climatology 
for  T„  T„  q,  and  p.  Vq)or  pressure  (e)  was  computed  as: 

e* - ^ - . 

(q  ♦  €  -  (&J)) 

where  e  =  0.622  (Gill,  1982).  The  cloud  correction  of  the  longwave  radiation  was  made 
using  uic  cloud  correcticm  factor  recommended  by  Clark  et  al.  (1974), 

LWu  -  LW„U{i-bC^, 

where  b  =  0.62  and  C  is  the  cloud  fraction  obtained  from  the  ISCCP  data.  Upward  and 
downward  arrows  indicate  the  net  of  the  upward  and  downward  components  of  longwave 
radiation.  Using  this  procedure  monthly  mean  values  of  net  surface  longwave  radiation 
were  computed  ftx*  the  48-month  observation  period. 

9.  Latent  Heat  Flux 

Of  the  four  terms  in  the  surface  heat  flux  in  the  central  equatorial  Pacific, 
the  latent  heat  flux,  (^,  is  typically  the  largest  source  of  heat  loss  from  the  ocean  surface 
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to  the  atmosphere.  Climatological  values  of  at  0”.  140^  computed  using  bulk 
formulae  range  between  about  -75  and  -110  Wm'^.  Using  available  data  from  several 
sources,  was  computed  as; 


Qt  = 

where  p,  is  the  density  of  air  (1.225  kg  m'^),  /„  is  the  latent  heat  of  valorization.  C,  is  the 
transfer  coefficient  for  latent  heat  or  Dalton  number  (1.15  x  10'^),  q,  is  the  saturation 
specific  humidity,  q^  is  the  specific  humidity  measured  near  the  sea  surface,  and  U  is 
the  surface  wind  speed  (Large  and  Pond.  1982;  Weller  and  Taylor,  1993).  The  latent  heat 
of  vaporization  was  computed  as: 

=  [2.501 -0.00237xrjxl0‘  Jkg  ' 

where  is  the  air  temperature.  The  saturation  specific  humidity  was  computed  using 
Tetens’  formula  (Buck,  1981); 


«  ^  17.502x7^ 

=  (1.007+3.46xlO-*xP)x6.1121xexi)( - ~-=r) 

240.97+7, 

q.  =  0.62197X-7; - 

(P-0.378xP^ 

where  is  the  viq^or  pressure  and  P  is  the  atmospheric  pressure.  Depending  on  data 
availability,  U  and  7,  were  measured  from  the  moored  buoys,  and  q„  and  P  wo’e  taken 
from  the  GOADS  climatology.  When  buoy  observations  of  U  or  7,  were  unavailable,  the 
climatology  of  McPhaden  and  McCarty  (1992)  was  used.  When  CO  ADS  observations  of 
or  P  were  unavailable,  the  Esbensen  and  Kushnir  (1981)  climatology  was  used. 

Since  observations  of  <7,  wnv  available  for  only  19  of  the  48  months  of 
this  study,  the  sensitivity  of  the  calculation  of  Qe  to  q„  was  examined  by  comparing  values 
of  Qe  computed  using  COADS  observations  of  q^  to  values  of  computed  using  the 
Esbensen  and  Kushnir  climatolt^cal  values  of  q^  (Figure  3.5).  The  root-mean  square 
(rms)  dffierence  was  21.31  Wm'*,  or  mtwe  than  25%  of  the  mean  (2,  fw  the  COADS 
data,  -77.4  Wm'^.  Despite  the  large  rms  difference,  the  avrelation  between  these  fields 
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is  0.67. 


10.  SenaUe  Heat  Flux 

Sensible  heat  flux.  Q„.  is  typically  the  smallest  term  of  the  surface  heat 
flux,  with  climatological  values  of  about  -2  to  -5  Wm'^.  Monthly  means  of  sensible  heat 
flux  were  computed  as; 

<?*  “  9/:,CjiTrT^U 

where  is  the  speciRc  heat  of  dry  air  (1004.67  m^s'^K'*),  C,  is  the  transfer  coefficient  for 
sensible  heat  (x  Stanton  number,  and  is  die  sea  surface  temperature.  A  Stanton  number 
of  LISkIO*’  was  used  for  unstable  conditions  (T^T^)  and  0.75x10'^  was  used  for  stable 
conditions  {T^<T^  (Large  and  Pond,  1982;  Weller  and  Taylor,  1993).  Based  on  monthly 
mean  T^-T^  differences,  stable  conditions  occurred  during  48%  of  the  months. 

11.  Precipitation 

Unlike  the  situation  over  land,  direct  measurements  of  precipitation  over 
the  oceans  are  rare,  particularly  over  the  central  equatorial  Pacific.  Three  atlases  of 
precipitation  have  been  used  to  provide  a  range  of  estimates  of  monthly  mean 
precipitation  near  (f,  140^  (Table  3.3).  Dorman  and  Bmirke  (1979)  derived  estimates 
of  precipitation  for  7*  latitude  by  5”  lon^tude  rectangles  over  the  Pacific  using  present 
weather  observations  taken  by  ships.  Using  the  method  of  Tucker  (1961),  measured 
rainfall  amounts  at  land  stations  were  related  to  the  present  weather  reported  in  the 
synoptic  observations.  Corrections  related  to  local  temperature  biases  were  also  applied. 
Over  the  1952-1979  period  of  their  study,  cmly  about  5(X)  observations  were  rqxxled  for 
each  of  the  2°  x  5"  rectangles  surrounding  0”,  14(fW.  They  rqxxted  that  this  is  the 
minimum  number  of  observations  fw  reliable  estimates.  Their  analysis  indicated  annual 
precipitation  of  about  8(X)  mm  (Figure  3.6).  Seasonally,  their  estimates  showed  increased 
precipitation  during  the  boeal  spring  and  decreased  precipitation  during  the  fall  and 
winter. 

Shea  (1986)  produced  a  global  i»rec4>itation  atlas  based  on  objective 
analysis  of  shore  station  data.  Oberhuber  (1988)  used  this  atlas  and  an  analysis  of  die 
surface  salinity  field  of  Levitus  (1982)  to  jx'oduce  a  prec^itadon  climatology.  In  the 
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Figure  3^.  Comparison  between  monthly  mean  values  of  latent  heat  flux  cmnputed  using 
specific  heat  acquired  from  monthly  GOADS  observatiims  and  latent  heat  computed  n«ing 
specific  heat  acquired  from  the  Esbensen  and  Kushnir  (1981)  climatology. 
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Figure  3.6.  Mean  annual  precipitation  over  the  Pacific  Ocean  between  30°S  and  60^ 
from  Dorman  and  Bourke  (1979). 
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LOCATION  or  STATIONS  USCO  FOR  PRCCIPITATION  COMPUTATIONS 
1*30  -  l»7*  (10) 

Figure  3.7.  Locations  of  stations  used  fm*  precipitation  atlases  of  Shea  (1986)  and 
Oberhuber  (1988)  for  the  period  1950-1979. 
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vicinity  of  0°.  140^,  these  analyses  are  based  on  widely-spaced  distant  island  stations 
(Figure  3.7).  The  Oberhuber  climatology  indicates  annual  precipitation  amounts  of  about 
1000  mm.  The  annual  cycle  is  characterized  by  precipitation  rates  of  about  130  mm  m^'‘ 
in  the  boreal  spring  and  about  35  mm  mth'‘  in  the  boreal  fall  (Table  3.3,  Figure  3.8). 
Considering  the  scarcity  of  observations  in  both  climatologies,  the  agreement  is  good. 
Both  atlases  indicate  a  strong  prec4)itation  gradient  north  of  the  equator  associated  with 
the  inter-tropical  convergence  zone  (ITCZ).  This  strong  gradient  suggests  the  possibility 
of  significant  variability  of  precipitation  as  the  ITCZ  varies  seasonally  and  interannually. 


Shea  (1966)/0berhuber  (1988)  Climatology 


SEASONAL  CYCLE  OF  PRECIPITATION  at  0,1 40W 


Figure  3,8.  Seasonal  cycle  of  precipitation  at  0",  140°W  from  Shea  (1986)/Oberhuber 
(1988). 
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Table  3J.  Climatological  values  of  precipitation  at  0°,  140^  frwn  Dorman  and  Bourke 
(1979),  Shea(1986)/Oberhuber  (1988),  arid  Janowiak  and  Arkin  (1991).  Units  are  mm 
mth'. 


Another  indication  of  the  magnitudes  of  precipitatim  and  the  seasonal  and 
interannual  variability  of  precipitation  at  (f,  14(fW  is  provided  by  the  satellite-doived 
precipitation  estimates  from  the  Global  Precipitation  Climatology  Project  (GPCP).  Using 
infrared  (IR)  data  from  geostatitniary  and  polar-orbiting  satellites,  Janowiak  and  Aiidn 
(1991)  derived  pentad  (5-day)  and  mtmdily  rainfall  estimates  in  the  tropics  during  1986- 
1989  using  a  simple  cloud-tq)  temperature-thresholding  algorithm.  Dr.  Janowiak  kindly 
provided  an  updated  time  series  of  predpitatimi  estimates  computed  by  the  same  method 
for  the  period  January  1986  through  December  1993  (Figure  3.9a).  These  data  show  high 
prec4>itation  rates  during  the  1986/1987  and  1991/1992  ENSO  warm  episodes  and 
precipitation  rates  considerably  lower  dian  the  Dorman  and  Bmirke  (1976)  and 
Shea(1986)/Oberhuber  (1988)  climatologies  during  non-ENSO  years.  The  eight-year 
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GPCP  time  series  of  monthly  mean  precipitation  estimates  was  used  to  calculate  the 
GPCP  climatolt^  shown  in  Table  3.3.  This  climatology  indicates  annual  precipitatitm 
at  0”.  14(fW  of  about  257  mm.  considerably  less  than  the  two  tMher  climatol<>gical 
estimates. 

Using  the  GPCP  climatology  for  the  period  Novembtf  1983  -  December 
1985  and  the  GPCP  monthly  mean  obso^raticms  for  the  period  January  1986  -  October 
1987,  a  cmnposite  time  series  of  precipitation  estimates  for  the  study  period  was  derived. 
The  composite  GPCP  time  series  and  the  Sl»a  (1986)/Oberhuber  (1988)  climatological 
time  series  are  compared  in  Hgure  3.9b.  Excq>t  for  a  sharp  peak  in  precipitation  during 
the  spring  of  1987,  the  GPCP  predpitadtm  estimates  are  consistently  much  smaller  than 
the  precipitation  estimates  given  by  the  Shea/Obeibuber  climatolt^es.  Since  both  the 
Shea/Oberhuber  and  Dorman  and  Bourke  climatologies  were  based  on  scarce  greater 
confidence  is  given  to  the  GPCP  data.  Unfortunately,  these  data  have  not  been  validated 
with  reliable  surface  measurements  in  the  central  equatorial  Pacific  ,  so  uncertainty 
remains.  Therefore,  buoyancy  flux  calculaticms  were  performed  using  precipitatitm 
estimates  firtxn  both  the  GPCP  time  series  and  the  Shea/Oberfauber  climatology. 
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1986  1987  1988  1989  1990  1991  1992  1993 
6PCP  Monthly  Precipitation  at  0,140W  (Jan86  ~  Dec93) 


Monthly  Precipitation  Estimates  at  0,140¥  (NotBS  -  0ct87) 


Figure  3.9.  a)  Time  series  of  sateilite-<ierived  monthly  mean  GPCP  prec4>itati<Hi 
estimates  for  the  period  January  1986  Uirough  December  1993  averaged  between  2.S*N 
and  2^"S  and  138.75^  and  141.25*W  (frmn  Dr.  J.  Janowiak).  b)  Comparison  between 
time  series  of  mondily  mean  precipitatioa  estimates  based  on  a  composite  of  GPCP 
climatology  and  observatioas  and  mi  the  SheaA3berhuber  climatology  for  the  period 
November  1983-October  1987.  Units  are  mm  mth'‘. 


48 


12.  Wave  Momentum  Flux 

As  discussed  in  the  introduction,  one  of  the  ctmcems  here  is  the  diurnal 
modulation  of  the  wave  induced  ztxial  momentum  flux.  The  Skyllingstad  and  Denbo 
(1994)  model  results  showed  that  internal  waves  can  provide  an  efficient  mechanism  for 
the  vertical  transport  of  horiz<mtal  mmnentum.  In  their  modeled  surface  mixed  layer,  the 
internal  waves  gained  eastward  momentum  at  the  expense  of  the  background  flow.  In 
stxne  cases,  this  momentum  was  transferred  back  to  the  mean  flow  at  a  critical  level 
resulting  in  a  deceleration  below  the  undo'current  core.  Otherwise,  the  waves  tended  to 
decrease  the  current  velocity  above  the  EUC  core.  In  an  effort  to  compare  their  model 
results  to  observations,  wave  momentum  flux.  u'w'  .  shear  production  of  turbulent 
kinetic  energy,  {u'w^yduldz  ,  and  wave  stress  divergence,  dfdziu'w') ,  were  computed  using 
the  mowed  buoy  data.  As  with  IDV  calculi^wis,  the  temperature  and  velocity  data  were 
effectively  high-pass  filtered  to  remove  the  dominant  tidal  effects  by  segmenting  the  data 
into  6-hour  time  bins  and  fitting  each  to  a  cubic  polynomial.  The  residuals,  t/  and  T, 
were  then  used  to  cwnpute  the  above  variables  as  follows: 


ii 


/ 


R-1,23 


<- 


(C, 


where 


C, 


^az 


au'w'  ^  -  «>Az2) 

*  "  AZ 

As  befwe,  the  Ar  fw  all  of  the  observaticms  was  15  minutes,  however,  the  individual 
observations  of  velocity  rq)resent  15-minute  averaged  quantities,  whereas  some  of  the 
tempwatures  are  15-minute  averages  and  some  are  15-minute  spot  san^les.  The  15- 
minute  averaged  temperatures  and  velocities  have  had  die  higher  frequendes  effectively 
removed  by  averaging.  In  the  calculation  of  w',  the  mean  vertical  temperature  gradient, 
dlT0z,  was  computed  using  centered  differencing.  In  the  calculation  of  the  shear 
production  of  TKE,  the  vertical  velocity  shear  was  computed  as  described  previously. 
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C.  SENSITIVITY  STUDY  OF  COMPUTED  IDV 


Since  the  primary  focus  of  this  investigation  is  an  examination  of  the  diurnal  cycle 
of  high-frequency  temperature  variability,  the  following  sensitivity  studies  were 
performed  to  verify  the  suitability  of  the  method  of  computing  monthly  mean  IDV.  Since 
the  sampling  period  of  IS  minutes  is  insufficient  to  resolve  all  of  the  high-frequency 
waves  of  interest,  it  was  necessary  to  rely  on  statistical  methods  to  examine  the  character 
of  these  high-frequency  waves.  It  was  hypothesized  that  with  sufficient  averaging  of  the 
variances  of  the  IS-minute  spot  sampled  temperatures  it  would  be  possible  to  enhance  the 
diurnal  signal  of  the  high-frequency  waves.  To  examine  this  hypothesis,  two  simple 
studies  comparing  the  computations  of  IDV  using  the  15-minute  temperatures  to  two 
high-frequency  data  sets  were  performed. 

1.  Comparison  With  TROPIC  HEAT  I  Results 

Values  of  IDV  were  cmnputed  for  the  month  of  November  1984  and  for 
the  week  of  November  25  -December  1, 1984  (Hgure  3.10).  This  latter  period  was  chosen 
to  coincide  with  the  TROPIC  HEAT  1  observation  period  when  Wijesekera  and  Dillon 
(1991)  presented  isopycnal  di^lacement  variance  as  a  function  of  time  of  day  (Figure 
3.10a).  Their  calculations  of  isopycnal  displacement  variance  were  conc^uted  from 
microstructure  profiles  obtained  using  the  riq>id  sampling  vertical  jM'ofiler  (RSVP) 
(Caldwell  et  al.,  1985).  Measured  variables  were  dq>th,  temperature,  salinity,  potential 
density,  and  turbulent  kinetic  eneigy  (TKE)  dis84>ation  rate.  The  high  vertical  and 
temptxal  resoluticm  of  thw  data  allowed  isopycnal  di8placeme.its  to  be  computed  directly 
fr(»n  individual  Thoipe-ardered  density  profiles  (Thoq)e,  1977).  Afra-  ranoving  the  0.1 
running  mean  from  these  profiles,  they  computed  IDV.  Comparison  between  their 
isopycnal  displacement  variances  (Figure  3.9a)  and  IDV  computed  hoe  (Figure  3.10b) 
shows  reasonable  agreement,  particularly  with  regard  to  the  vertical  structure  of  the 
diurnal  variations.  Generally,  the  vertical  structure  is  similar,  but  the  magnitudes  of  their 
displacement  variances  are  larger  by  about  a  hmtor  of  two.  The  reason  for  the  difference 
in  magnitudes  remains  unclear,  but  could  be  partially  due  to  the  higher  vertical  resolution 
of  their  data,  the  indtnion  of  salinity  in  dieir  observations,  and  their  shorter  sampling 
interval  of  7.5  minutes,  as  opposed  to  die  IS-minute  sampling  intoval  used  hoe.  The 
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bopycnal  Displaceinent  Variance  (m^ 
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IDV  at  0",  14(nv  (Nov  25  •  Dec  1984) 


Figure  3.10.  Comparison  between  a)  mean  isopycnal  displacement  variance  from 
TROPIC  HEAT  microstmctuie  measurements  by  Wijesekera  and  Dillon  (1991)  during 
November  2S-December  1.  1984  at  QT,  lAffW,  and  b)  mean  isotherm  diylacement 
variance  (IDV)  for  four  six-hour  time  bins  using  moored  buoy  data  from  the  same  time 
period  and  location. 
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magnituctes  show  a  [xt>nounced  diumai  cycle  of  isopycnal  displacement  variance 
between  IQS  m  and  130  m  that  is  not  evident  in  the  IDV  results  ccanputed  from  the 
mooing  data.  Given  the  resolution  differences,  and  the  fact  that  in  this  case  die  mooring 
data  were  averaged  only  over  seven  days,  the  corre^ndence  between  IDV  results  is 
noteworthy. 

2.  Comparison  With  l>Miniite  TR  Data 

Although  the  comparison  with  Wijesekera  and  Dilloi  showed  reasonable 
agreement,  differences  are  apparent  and  the  results  are  not  conclusive.  Therefore,  1- 
minute  spot-sampled  temperature  data  at  40  m  dq>th  for  the  period.  May  12  -  June  1 1, 
1987,  were  used  to  further  validate  our  methodology.  This  is  a  subset  of  the  data 
examined  by  McPhaden  and  Peters  (1992).  who  showed  a  diumai  cycle  of  IDV  at 
frequencies  near  N  (10-30  q)h).  The  identical  procedure  for  calculating  IDV  with  die  15- 
minute  data  was  used  to  calculate  IDV  with  the  1-minute  data,  including  segmenting  the 
data  into  six-hour  bins  and  fitting  the  data  in  each  bin  with  a  cubic  polynomial. 
Temperature  variances  were  then  ctmiputed  from  the  6-hour  residuals,  which  in  this  case 
were  composed  of  360  data  points.  Tliese  variances  were  then  normalized  by  the  vertical 
temperature  gradient  to  compute  IDV.  Thus,  30  daily  IDVs  were  computed  for  each  of 
the  four  6-hour  time  bins.  Using  these  1-minute  IDVs  as  a  standard,  this  (vocedure  was 
then  rqieated  after  subsampling  the  l-minute^iata  in  two  ways.  First,  for  each  15-minute 
period,  the  1-minute  temperatures  were  randomly  subsampled.  Values  of  IDV  were  then 
computed  from  the  resulting  unevenly  spaced  time  series  using  the  above  procedure.  The 
resulting  30  IDV  values  for  each  time  bin  were  then  correlated  with  the  1-minute  IDV 
values  (Table  3.4).  Corelations  were  str<Mig  with  values  ranging  between  0.9173  and 
0.9867.  Next,  the  1-minute  temperatures  were  subsampled  every  15  minutes  and  the  IDVs 
and  correlations  with  the  1 -minute  IDVs  woe  omputed  (Table  3.4).  Again,  correlations 
were  strong  with  values  ranging  between  0.9100  and  0.9858.  Both  sets  of  correlation 
coefficients  indicate  strong  agreement  between  the  1-minute  IDVs  and  foe  15-minute 
IDVs,  indicating  foat  the  methodology  is  suitable  for  characterizing  the  variability  of  high- 
frequency  intonal  waves. 
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Table  3.4.  Correlation  coefficients  between  30  day  time  series  measurements  of  1 -minute 
temperature  data  at  40  m  depth  and  IS-minute  randomly  subsampled  and  periodically 
subsampled  temperatures  for  the  four  local  time  bins  listed. 


LOCAL  TIME  BIN 

15*ininute  RANDOM 

15-Minute  PERIODIC 

2200  -  0400 

0.9173 

0.9100 

0400  -  1000 

0.9472 

0.9535 

1  1000  -  1600 

0.9867 

0.9858  1 

1  1600  -  2200 

0.9843 

0.9168  1 
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IV.  BACKGROUND  CONDITIONS 


A.  INTERANNUAL  VARIABILITY:  ENSO 

Oceanic  and  atmospheric  conditions  in  the  central  equatorial  PaciHc  exhibit 
significant  temporal  variability  over  scales  ranging  from  diurnal  to  interannual.  The  ENSO 
phenomenon  is  the  most  notable  and  pronounced  example  of  interannual  global  climate 
variability.  A  major  ENSO  warm  q)isode,  such  as  that  which  occurred  during  1982/1983, 
leads  to  massive  dislocations  of  tropical  rainfall  regimes,  bringing  drou^t  to  vast  areas 
and  torrential  rains  to  otherwise  arid  regions.  The  related  atmospheric  circulation 
anomalies  extend  deep  into  the  extratrc^ics,  where  they  e  associated  with  imusual  and 
sometimes  severe  climatic  conditions  over  a  variety  of  regions  around  the  globe. 

Although  each  ENSO  warm  episode  has  unique  characteristics,  with  q>isodes 
difrering  in  strength,  phase,  and  behavior,  it  is  useful  to  describe  a  composite  evolution 
over  a  period  of  18-24  months.  The  anomalous  ocean  wanning  in  the  Pacific  often  begins 
near  the  Ecuador-Peru  coast  early  in  the  calendar  year,  tiben  spreads  westward  into  the 
central  equatorial  Pacific.  The  coastal  warming  often  peaks  during  April- June,  but  the 
warming  in  the  central  equatorial  Pacific  nrxmally  continues  fw  several  more  months,  as 
the  high  SSTs  of  the  western  Pacific  spread  eastward.  The  global  atmospheric  climate 
anomalies  are  most  widespread  and  intense  near  the  end  of  the  Erst  year  and  during  the 
early  months  of  the  second  year  of  the  episode,  i.e.,  during  the  boreal  winter.  This  period, 
t^roximately  one  year  following  the  initial  ^)pearance  of  warm  water  in  the  eastern 
Pacific,  is  often  referred  to  as  the  mature  phase  of  the  qiisode.  Following  the  mature 
phase,  the  anomaly  patterns  enter  a  period  of  slow  decay  that  spans  several  months. 

In  November  1983,  the  beginning  of  our  study  period,  the  monthly  Climate 
Diagnostics  Bulletin.  rqxN^  Aat  the  global  1982/1983  ENSO  warm  evmit  had  ended. 
From  November  1983  through  early  1985,  c(mditi<ms  of  the  tropical  Pacific  were 
indicative  of  ENSO  cold  conditions.  During  this  period,  satellite-sensed  mitgoing 
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longwave  radiation  (OLR)  was  anomalously  positive,  indicating  relatively  dry  conditions. 
In  trc^ical  regions,  OLR  provides  an  indication  of  the  intensity  of  atmospheric  convection 
and  the  related  precipitation.  Lower  level  winds  were  anomalously  easterly  in  the  central 
equatorial  Pacific.  The  Southern  Oscillation  Index  (SOI),  computed  from  differences  of 
sea  level  pressure  anomalies  at  Tahiti  and  Darwin,  was  near  normal  over  much  of  this 
period.  Temperature  measurements  at  0®.  140"W  revealed  anomalously  low  SSTs  and 
shallow  thermocline  depths  throughout  this  period.  Progressing  through  1985,  conditions 
across  the  trt^ical  Pacific  suggested  a  trend  toward  the  build-up  phase  of  ENSO  that 
often  precedes  ENSO  warm  episodes  (Climate  Diamostic  Bulletins.  1985).  This  phase 
was  characterized  by  slowly  increasing  SOI,  anomalous  low  level  easterlies,  below 
normal  equatorial  SST,  high  sea  level  in  the  western  Pacific,  low  sea  level  in  the  eastern 
Pacific,  low  convective  activity  in  the  central  Pacific,  and  high  convective  activity  in  the 
western  Pacific.  In  early  1986,  several  of  the  ENSO  indices  suggested  the  possible 
development  of  an  ENSO  warm  episode,  but  by  May  1986  conditions  across  the  txqpical 
Pacific  had  returned  to  near  normal,  where  they  remained  through  the  btxeal  summer. 
During  the  fall  of  1986,  most  of  the  atmospheric  and  oceanic  indices  indicated  the 
develq)ment  of  the  1986/1987  ENSO  warm  episode.  By  October  1986,  positive  SST 
anomalies  were  observed  across  much  of  the  equatcxial  Pacific,  weak  westerly  wind 
anomalies  were  observed  in  the  central  and  eastern  equatorial  Pacific,  and  negative  OLR 
anomalies  were  found  near  the  date  line  at  die  equator.  By  November  1987,  most  of  the 
indices  had  started  a  slow  return  to  normal  conditions  across  the  equaUxial  Pacific. 

B.  INTERANNUAL  VARIABILITY;  LOCAL  CONDITIONS 

In  this  section,  local  surface  forcing  and  oceanic  conditions  observed  during  the 
four-year  period  chronicled  in  the  previous  section  on  large-scale  interannual  ENSO 
variability  are  presented.  Monthly  mean  time  series  of  the  following  observed  or 
computed  variables  for  the  48-month  rectxd  are  presented:  ztmal  and  meridional  winds, 
wind  speed,  shortwave  radiation,  longwave  radiation,  sensible  heat  flux,  latent  heat  flux, 
net  surface  heat  flux,  surface  buoyancy  flux,  subsurface  tempo-ature  (0-250  m),  subsurface 
zonal  and  meridional  velocity  (10-250  m),  buoyancy  frequency  (10-250  m),  velocity  shear 
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(10-250  m),  gradient  Richardson  number  (10-250  m),  mixed-layer  depth,  diurnal  range 
of  mixed-layer  depth,  and  isotherm  displacement  variance  (35  m,  60  m,  1(X)  m.  140  m, 
and  2(X)  m). 

1.  Surface  Fordng 

Monthly  mean  winds  from  the  motxing  at  0°,  140°W  are  shown  in  Figure 
4.1.  The  solid  curves  are  monthly  means  of  the  observed  winds  and  dotted  lines  represent 
the  climatological  mean  annual  cycle  computed  from  mooring  data  during  1983-1991 
(Mcl%aden  and  McCarty.  1992).  Monthly  mean  winds  were  anomalously  weak  (westerly) 
during  the  period  following  the  1982/1983  ENSO  warm  event  through  September  1984 
and  again  starting  in  the  fall  1986/Mdnter  1987  and  lasting  through  the  1986/1987  warm 
episode.  Anomalous  easterlies  prevailed  during  the  later  stages  of  the  1984  cool  episode 
and  through  the  1985/1986  ENSO  build-up  phase.  Near  normal  zonal  winds  during  the 
1986/1987  winter  were  followed  by  a  colli4)8e  of  the  zonal  winds  during  the  mature  phase 
of  the  1986/1987  ENSO  warming.  Since  weak  local  winds  were  observed  during  both 
the  1984  cold  episode  and  the  mature  phase  of  the  1986/1987  warm  episode,  ctxrelation 
between  local  surface  winds  and  the  lar^r  scale  ENSO  conditions  did  not  appear  strong, 
suggesting  non-local  forcing  of  ENSO  events. 

Monthly  mean  estimates  of  Q^,  Q,,  Q^,  and  Q„,  are  presented  in  Figures 
4.2  and  4.3.  In  Figure  4.2,  solid  curves  represent  computed  monthly  mean  estimates  and 
dotted  lines  rq)resent  climatological  valu^.  Values  of  Q,,  and  Q,  (Figure  4.2b,c)  with 
closed  circle  symbols  were  based  on  COADS  mondily  means  or  climatological  values 
of  air  temperature  and  specific  humidity,  respectively.  The  climatological  values  of  Q^, 
Qi,,  and  C2«  are  from  the  Oberhuber  (1988)  climatology,  which  is  based  on  the  2°  x  2° 
COADS  data  set  for  the  period  1950-1979.  The  climatology  of  Q,  was  computed 
following  Reed  (1977)  using  COADS  cloud  data  and  buoy  data  for  the  remaining  terms. 
The  climatology  of  Q„,  was  computed  from  the  climatologies  of  the  four  terms. 

Net  longwave  radiation  was  negative  throughout  the  study  period,  with  a 
mean  net  loss  of  50.2  Wm'^from  the  ocean  to  die  atmosph^-e  (Figure  4.2a).  Comparison 
with  Figure  4.4a  shows  the  dependence  of  estimated  longwave  radiation  on  SST,  and 
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Monthly  Mean  U,  V  and  Wind  Speed  at  0,  140  W 


1983  1964  1985  1986  1967 


Figure  4.1.  Time  series  of  monthly  mean  values  of  a)  meridional  wind,  b)  zonal  wind 
and  c)  wind  speed  at  0",  140*W  during  the  period  November  1983  -  October  1987.  Solid 
lines  indicate  monthly  mean  observations.  Dotted  lines  indicate  clinifttological  means  from 
McPhaden  and  McCarty  (1992). 
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Shortwave,  Longwave.  Latent,  Sensible  &  Net  Heat  Flux 


1983  1984  1985  1988  1987 


Figure  4  Monthly  mean  surface  heat  flux  terms  at  0”,  14(y^V  for  the  period  November 
1983  -October  1987:  a)  longwave  ra(liati<m,  b)  sensible  heat  flux,  ,  c)  latent 
heat  flux,  Q^,  d)  shortwave  radiation,  by  FAST  radiative  scheme  (solid  line)  and 
BULK  formula  (dashed  line),  and  e)  net  surface  heat  flux,  Climatological  values 
are  given  by  dotted  lines.  Downward  fluxes  are  positive. 
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therefore  the  ENSO-related  interannual  variability  of  SST.  During  the  1984  ENSO  cold 
episode  net  outgoing  longwave  radiation  was  reduced  to  about  -40  to  -45  Wm'^.  Low 
values  of  longwave  radiation  were  during  the  build-up  phase  of  the  1986/1987  ENSO 
wanning.  During  the  mature  phase  of  the  warm  episode,  net  outgoing  longwave  radiation 
was  increased  to  about  -55  to  -60  Wm*^ . 

Sensible  heat  flux  estimates  varied  between  +5  and  -5  Wm'^  (Figure  4.2b). 
These  estimates  were  typically  about  2-7  Wm'^  greater  than  the  climatological  estimates 
of  Oberhuber  (1988),  which  were  always  negative  between  about  -2  and  -5  Wm  ^  The 
small  relative  magnitudes  suggest  that  sensible  heat  flux  had  only  minm'  relative 
significance  to  the  generation  of  strong  convection  by  diurnal  changes  in  sensible  heat 
(Figure  4.3).  Generally,  sensible  heat  flux  was  negative  during  the  cold  episode  and 
positive  during  the  warm  episode.  Unfortunately,  air  temperature  measurements  were 
available  for  only  23  of  the  48  months  of  the  study,  so  climatological  values  of  T,  from 
McPhaden  and  McCarty  (1992)  were  necessary  fcr  much  of  the  time  series.  Although  not 
shown,  time  series  of  15-minute  SST-T,  values  show  significant  variations  at  periods  from 
hours  to  weeks. 

Latent  heat  flux  represents  the  largest  of  the  surface  heat  loss  terms  (Figure 
4.3),  with  an  average  monthly  mean  value  of  -79.5  Wm'^.  With  a  standard  deviation  of 
25.0  Wm'^,  the  variability  of  is  the  Itf^t  of  the  four  heat  flux  terms.  Latent  heat  flux 
estimates  range  between  about  -30  and  -135  Wm'^  (Figure  4.2c).  Fran  March  1984 
through  March  1985,  latent  heat  flux  loss  estimates  of  -25  to  -70  Wm'^  were  about  20-50 
Wm'^  less  than  climatological  values.  From  about  May  1985  through  Ajvil  1987,  latent 
heat  estimates  were  near  climatological  values  of  about  -75  to  -120  Wm'^.  Latent  heat  flux 
estimates  varied  widely  between  about  -35  and  -1 15  Wm'^  during  the  mature  phase  of  the 
1986/1987  warm  episode. 

Shortwave  radiation  is  the  largest  term  of  the  surface  heat  balance  with 
average  values  of  283.7  Wm'^  (Figure  4.3).  Shortwave  radiation  computed  from  both  bulk 
formulae  (dashed  line,  BULK)  and  the  radiative  transfer  scheme  (solid  line,  FAST)  of 
Bish(^  and  Rossow  (1991)  are  presented  in  Figure  4.2d.  Both  BULK  and  FAST 
estimates  were  consistently  greater  that  the  Oberhuber  (1988)  climatological  values  (not 
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shown),  by  71.8  sml  49.8  Wm'^,  respectively.  The  BULK  estimates  were  lower  than  the 
FAST  estimates  by  about  20  Wm'^.  The  shortwave  radiation  values  computed  using  the 
bulk  formulae  were  erratic  and  cannot  be  simply  related  to  seasonal  or  ENSO  interannual 
variability.  The  PAST  shortwave  radiation  estimates  show  a  well-defined  seasonal  cycle 
with  a  semi-annual  period  and  interannual  characteristics  apparendy  related  to  ENSO. 
The  seasonal  cycle  appeared  gready  diminished  during  the  1984  cool  q>isode,  and  there 
was  a  significant  reduction  in  downward  shortwave  naliatitm  during  the  1986/1987  warm 
episode.  Since  dus  reductitm  in  incoming  shmtwave  radiatimi  during  the  warm  event  did 
not  i^pear  in  the  simpler  bulk  formula,  which  takes  into  account  cloud  fraction  but  not 
cloud  type  or  vertical  structure  in  the  atmospheric  colunm,  it  appears  that  their  inclusion 
in  the  radiative  transfer  model  provides  important  imin-ovements. 

The  monthly  mean  net  surface  heat  flux.  was  computed  from  the 
monthly  mean  heat  flux  terms  described  above  (Figures  4.2e  and  4.3).  The  average  net 
surface  heat  flux  into  the  ocean  was  1S4.5  Wm'^.  with  monthly  mean  values  ranging 
between  223  Wm'^  and  78  Wm'^.  The  standard  deviation  of  the  monthly  mean  heat  flux 
was  29.7Wm'^.  indicating  large  month  to  month  variability.  Figure  4.3  shows  the  relative 
contribution  of  the  separate  terms  to  the  mt  surfac:«  heat  flux.  As  mentioned  for  the 
individual  terms.  is  the  dominate  surface  heat  source  and  and  are  tiie  dmninate 
heat  loss  terms.  The  high  variability  of  Q„,  is  largely  determined  by  the  variability  of  Q^, 
and  to  a  lesser  extent  Q^.  Interaimually.  Q„,  was  relatively  high  during  tiie  1984  cold 
episode,  with  values  averaging  about  170  Wm  ^  and  relatively  low  during  the  1986/1987 
warm  episode,  with  values  averaging  about  130  Wm'^. 
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Downward  Heat  Flux  (Wm 


Surface  Heat  Balance  at  0,1 40W  (Nov83  —  0ct87) 


Figure  43.  Monthly  mean  values  of  the  terms  of  the  surface  heat  balance  at  0°,  140^ 
during  November  1983-October  1987.  Terms  include  shortwave  radiation,  longwave 
radiation,  sensible  heat  flux,  latent  heat  flux,  and  net  surface  heat  flux, 
Downward  fluxes  are  positive. 
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The  surface  buoyancy  flux  determines  the  turbulent  vertical  motion  and 
amount  of  convective  mixing  in  the  ui^)er  ocean.  Figure  4.4  shows  mcmthly  mean 
estimates  of  the  net  surface  buoyancy  flux,  and  the  heat  and  moisture  components 

of  the  buoyancy  flux,  B,^,  and  B— _ _ respectively.  These  estimates  were  based  on  the 

above  estimates  of  the  net  surface  heat  flux,  and  calculations  of  the  surface  moisture 
flux.  Surface  moishue  flux  was  computed  using  both  the  composite  of  GPCP  observations 
and  climatology  and  the  Shea/Oberhuber  climatology  of  monthly  mean  precipitation,  as 

deso’ibed  in  section  3.2.11.  Although  B^.-. _ values  computed  using  the  different 

precipitation  estimates  were  quantitatively  different,  the  magnitudes  in  each  case  were 
approximately  an  (xxler  of  magnitude  less  than  B,^.  Even  during  the  large  precipitation 
event  during  the  spring  of  1987,  B„^  was  largely  dominated  by  B,^. 

The  average  value  of  B^u  was  -I.ISS  x  10'^  Wkg*'.  These  estimates  are  in 
rough  agreement  with  the  observed  daytime  and  nighttime  values  of  surface  buoyancy 
flux  during  the  TROPIC  HEAT  e]q)eriment,  -6  x  10*’  Wkg  ‘  and  +1.5  x  lO"’  Wkg  ', 
respectively  (Mourn  et  al.,  1989).  The  standard  deviation  of  monthly  mean  B^^  was  1.965 
X  lO'*  Wkg'.  Since  Bg^  was  determined  almost  entirely  by  the  surface  heat  flux,  the 
seasonal  and  interaimual  variability  of  B.^  was  almost  identical  to  that  described  for  Q^,. 
The  downward  buoyancy  flux  was  greater  than  average  during  1984  and  early  1985  and 
less  than  average  during  the  1986/1987  warm  episode. 
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Surface  Buoyancy  Balance  at  0,1 40W  (Nov83  —  0ct87) 


Figure  4.4.  Mcmthly  mean  values  of  the  tenns  of  the  surface  buoyancy  balance  at  (f, 
14(yw  during  November  1983-October  1987.  Terms  include  contributimi  to  buoyancy 

flux  by  moisture  flux,  B,.., _ heat  flux,  and  net  buoyancy  flux,  B,,.,.  Upwaid 

buoyancy  fluxes  are  positive. 
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2.  Subsurface  Temperature 

The  subsurface  thermal  field  from  the  surface  to  250  m  for  the  four-year 
study  period  is  presented  in  Figure  4.5a.  The  data  were  smoothed  with  a  17-day  triangular 
low-pass  filter  (20-  day  half-power  point)  to  rraiove  higher  frequency  oscillations  and  to 
emphasize  seasonal  and  interannual  variability.  The  subsurface  temperatures  show  a 
progression  from  cool  conditions  in  1984,  as  indicated  by  cool  SSTs  and  a  shallow 
thermocline,  followed  by  the  build-up  phase,  which  was  characterized  by  a  steadily 
deepening  thermocline.  and  ending  with  the  moderate  1986/1987  ENSO  warm  q}isode. 
which  was  characterized  by  warm  SSTs  and  a  deq>  thermocline.  The  depth  of  the 
thermocline  is  often  described  by  the  depth  of  the  20  C  isotherm,  which  is  shown  to 
follow  the  center  of  the  Aermocline  in  Figure  4.5a. 

In  addition  to  the  interannual  variations,  the  thermocline  underwent  depth- 
coherent  60-90  day  osdllations  that  were  particularly  apparent  during  the  boreal  winters 
of  1984  and  1986.  These  60-90  day  oscillations  were  rqxxled  by  McPhaden  and  Taft 
(1988)  and  ascribed  to  eastward  propagating  equatorial  Kelvin  waves.  Kessler  et  al. 
(1994),  in  analyzing  data  for  the  period  1983-1993  at  0®,  140“W,  report  that  feese 
intraseasonal  waves  are  typically  more  prevalent  during  the  boreal  winter  months  and  that 
interarmual  variability  of  these  waves  is  not  uncommon. 
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1984  1985  1986  1987 


Figure  4^.  Time  soies  measurements  of  a)  temperature  and  b)  zonal  current  velocity 
at  0",  140W  for  the  poiod  November  1983  -  Octobo- 1987,  These  plots  were  romputed 
from  daily  mean  values  which  were  then  smoothed  using  a  17-day  triangular  low-pass 
filter. 
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3.  Subsurface  Velocity 

The  subsurface  zonal  velocity  fteld  from  10  m  to  250  m  for  the  four-year 
study  period  is  presented  in  Figure  4.Sb.  As  with  the  subsurface  temperature  field,  the 
data  were  smoothed  with  a  17-day  triangular  low-pass  filter.  The  annual  vertical  migration 
of  the  EUC  core,  ranging  from  about  75  m  during  the  boreal  spring  and  summer  months 
to  about  120  m  during  fall  and  winter  months,  was  observed  to  vary  interannually  as  well. 
This  annual  cycle  was  also  accompanied  by  an  intensification  during  the  spring/summer 
shoaling.  During  the  1984  cold  episode,  the  spring/summer  shoaling  of  the  EUC  core  was 
more  pronounced  than  during  the  other  years.  During  the  1986/1987  ENSO  warm  episode, 
the  usual  fall/winter  relaxation  and  deepening  of  the  EUC  core  were  particularly 
pronounced,  with  the  EUC  velocity  reduced  to  as  low  as  25  cms'‘  during 
January/Fbbruary  1987.  As  observed  in  the  subsurface  temperature  data,  the  zonal  velocity 
experienced  pronounced  fluctuations  at  periods  of  60-90  days,  with  peak  amplitudes  of 
about  20  cms'‘. 

4.  Buoyancy  Frequency 

Monthly  mean  fields  of  buoyancy  frequency  (N),  velocity  shear  (S)  and 
gradient  Richardson  Number  (Ri)  during  November  1983  through  October  1987  are 
contoured  in  Figures  4.6a-c.  As  with  subsurface  temperatures,  N  showed  significant 
interannual  variability  associated  with  ENSO.  A  core  of  high  N  {N  >  0.014  s  ')  was 
observed  to  follow  the  depths  of  the  thermocline  shown  in  Figure  4.5a.  The  ctm  of  high 
N,  indicating  strong  stratification,  was  shallow  (0(50  m))  during  the  1984  cold  event. 
During  the  ENSO  build-up  phase  frmn  late  1984  until  early  1986,  the  high-N  core  slowly 
and  steadily  deepened  from  about  90  m  to  about  110  m,  with  monthly  maxima  ranging 
between  0.016  s  '  and  0.018  s  '.  During  the  1986/1987  ENSO  warm  episode,  the  high-N 
core  increased  in  magnitude  to  0.022  s  '  and  deepened  abruptly  in  December  1986/January 
1987  to  about  150  m.  apparently  related  to  the  passage  of  an  intense  equatorial  Kelvin 
wave.  Since  these  data  consist  of  monthly  means,  60-90  day  intraseasonal  fluctuations  are 
not  as  readily  apparent  as  in  the  subsurface  temperature  figure  (Figure  4.5a).  In  addition 
to  this  progression  of  the  high-N  core,  it  is  useful  to  examine  the  temporal  variability  of 
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the  low-N,  weakly  stratified  upper  SO  m.  At  times  when  N  was  relatively  high  in  the 
surface  layer,  indicating  strong  stratification,  vertical  turbulent  mixing  is  e]q>ected  to  have 
been  low.  In  early  1984,  the  upper  ocean  was  well-stratified  up  to  at  least  10  m.  and 
mixing  rates  would  be  expected  to  have  been  low.  When  N  was  low  in  the  upper  layer, 
such  as  during  the  1986/1987  ENSO  warm  qpisode,  increased  vertical  mixing  is  expected 
to  have  occurred. 

5.  Velocity  Shear 

The  monthly  mean  values  of  the  total  vertical  velocity  shear  contoured  in 
Figures  4.4a-c  include  both  zonal  and  meridional  velocity  components.  As  with  zonal 
velocity  (Figure  4.3b).  velocity  shear  underwent  significant  interannual  variability 
associated  with  ENSO.  During  early  1984,  velocity  shears  in  the  upper  10-40  m 
associated  with  the  shoaling  of  the  EUC  during  the  1984  cold  q>isode  exceeded  0.03  s'. 
By  June  of  1984  the  core  of  high  velocity  shear  had  weakened  and  deq)ened  to  about 
0.02  s  '  and  about  60  m.  respectively.  R'om  November  198S  through  January  1986  a 
pronounced  peak  in  5  was  observed  at  35  m,  with  a  maximum  S  exceeding  0.03  s  '.  The 
November  1985  and  January  1986  zonal  velocity  profiles  show  an  unusual  eastward 
acceleration  of  the  flow  at  45  m.  During  the  1986/1987  ENSO  warm  ^isode,  the 
weakening  of  the  EUC  resulted  in  significant  deepening  and  reduction  in  strength  of  the 
high-5  con.  The  differences  between  5  during  the  ENSO  cold  hnd  warm  q}isodes  were 
large. 

6.  Richardson  Number 

As  discussed  in  Chapter  n,  the  gradient  Richardson  Number,  Ri,  is  (tefined 
as  the  ratio  of  to  5^.  Ri  provides  a  measure  of  the  dynamic  stability  of  the  water 
column.  Figure  4.6c  shows  contours  of  the  monthly  mean  values  of  logio(R0  over  the 
four-year  period  of  study.  Since  these  values  are  3()-day  averages  of  RVs  computed  from 
6-hour  mean  temperature  and  velocity  profiles  based  on  limited  vertical  resolution 
observations,  it  is  expected  that  computed  mcmthly  mean  Ri  will  be  higher  than  that 
observed  on  microstructure  scales.  Each  of  the  (^)erations  performed  in  computing  diese 
Ri  values  tends  to  smooth  low  Ri  values.  It  is  not  unreasonable  to  assume  that  values  of 
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Ri  less  than  unity  might  r^resent  conditions  of  marginal  stability  in  a  monthly  mean 
sense  (as  opposed  to  a  critical  =  0.25  in  a  local  sense).  The  heavy  black  contour  line 
indicating  a  logioC^O  =  0.0  corresponds  to  Ri  equal  to  unity.  The  darker  regions  inside  the 
Ri  =1.0  contour  thoefore  represent  regitms  having  Ri  <  1  and  a  possibility  of  marginal 
dynamic  stability.  The  dashed  contours  inside  the  darkly  shaded  region  represent  Ri 
contours  of  0.63.  0.40,  and  0.25,  respectively. 

As  one  would  expect,  Ri  contours  exhibit  a  pattern  of  ENSO  interannual 
variability  similar  to  that  observed  for  N  and  S.  Since  this  study  is  concerned  with 
relationships  between  the  diurnal  variability  of  high-frequency  waves  and  turbulent  mixing 
by  dynamic  instabilities,  increased  attention  is  given  to  low-/?/  r^ons  above  the  EUC. 
The  low-/?/  region  underwent  an  annual  deq>ening  during  the  boreal  winter  months  down 
to  about  1(X)  m  and  a  shoaling  during  the  boreal  summer  months.  The  low-/?/  core 
qjpeared  to  shoal  to  depths  above  40  m  during  the  1984  cold  q>isode,  or  about  35  m 
shallower  than  usual.  During  the  1986/1987  warm  q>isode.  the  low-/?/  region  dis^)peared 
as  the  weakening  of  the  EUC  decreased  the  S  to  small  values. 

Each  of  die  Figures  4.6a-c  has  the  symbols  THl  and  TH2  along  the  upper 
time  axis  to  indicate  the  times  of  dm  TROPIC  HEAT  1  and  TROPIC  HEAT  2 
microstructure  measurements.  It  is  noteworthy  that  the  two  8a^^>ling  periods  were 
characterized  by  largely  different  oceanic  ctxiditiais,  as  rqireaenied  by  N,  5,  and  Ri. 
FurthermOTe,  these  figures  indicate  that  considerable  ncaanir  variability  occurs  at  (P. 
140°W,  and  it  is  not  unreasonable  to  expect  corresponding  variability  of  the  diurnal  cycle 
of  high-frequency  temperature  variability. 
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Figures  4.6a-c.  Time-depth  contour  plots  of  monthly  mean  values  of  a)  buoyancy 
frequency,  N  b)  velocity  shear,  S  and  c)  gradient  Richardson  number,  Ri,  from  10  m  to 
250  m  over  the  period  November  1983  -  October  1987.  THl  and  TH2  at  the  top  of  each 
plot  indicate  the  approximate  time  of  the  TROPIC  HEAT  experiment  microstructure 
measurements. 
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7. 


Mixed  Layer 

Mixed-layer  depfH  is  related,  directly  or  indirectly,  to  each  of  the  observed 
forcing  or  background  conditions  discussed  above.  Monthly  mean  estimates  of  mixed- 
layer  depth  for  each  of  the  four  six-hour  time  bins  are  presented  in  Figure  4.7a.  These 
mixed-layer  depths  were  computed  as  the  depth  where  the  vertical  temperature  gradient 
fu^t  reached  O.OlCm  '.  Although  the  10-15  m  vertical  instrument  spacing  limits  the 
accuracy  of  these  estimates,  the  diurnal  variations  were  consistent  throughout  the  record 
and  agree  with  diurnal  variations  of  mixed-layer  depth  computed  using  higher  resolution 
data.  Monthly  means  of  the  maximum  diurnal  ranges  in  mixed-layer  d^th,  i.e.,  difference 
between  daily  maximum  and  minimum  mixed-layer  depths  among  the  four  time  bins,  are 
presented  in  Figure  4.7b.  These  figures  show  considerable  interannual  variability  in  both 
mixed-layer  depth  and  the  diurnal  range  of  mixed-layer  depth  superimposed  over  the 
annual  cycle.  The  nighttime  mixed-layer  depth  had  a  mean  of  31.8  m  and  a  standard 
delation  of  15.1  m.  The  daytime  mixed-layer  dq)th  had  a  mean  of  23.1  m  and  a  standard 
deviation  of  13.7  m. 

During  the  1984  ENSO  cold  episode,  mixed-layer  depth  was  shallow, 
averaging  about  20  m  for  the  period  November  1983  through  May  1984.  During  this 
period,  the  diurnal  range  in  mixed-layer  depth  was  similarly  very  small  (0(5  m)). 
Following  this  cold  episode,  mixed-layer  depths  and  diurnal  ranges  of  mixed-layer  depths 
progressively  increased  through  the  ENSO  build-up  phase.  After  a  poiod  of  shoaling 
between  February  1986  and  ^ril  1986,  mixed-layer  depths  deepened  significantly  during 
the  1986/1987  ENSO  warm  episode.  The  diurnal  range  of  mixed-layer  deptiis  followed 
a  similar  progression. 

Since  surface  layer  mixing  processes  are  possibly  related  to  high-frequency 
temperature  variability  in  the  thermocline.  statistical  relationships  between  the  monthly 
mean  surface  forcing  and  backgroimd  oceanic  conditions  and  mixed-layer  dq>ths  are 
examined.  Of  the  fourteen  surface  forcing  variables  examined,  monthly  mean  mixed-layer 
depths  were  significantly  correlated  at  the  95%  confidence  level  only  with  wind  speed, 
zonal  wind,  meridional  wind  (for  day  and  evening  time  bins  only),  latent  heat  flux 
(evi^KMration),  net  surface  heat  flux,  and  net  surface  buoyancy  flux  (Table  4.1).  Ftx  these 
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Figure  4.7,  Time  series  of  monthly  mean  values  of  a)  mixed-layer  d^th  and  b)  diurnal 
range  of  mixed-layer  depth  over  the  poiod  November  1983  -  October  1987.  Mixed-layer 
depths  are  given  for  each  of  the  four  6-hour  time  bins:  2200-0400, 0400-1000, 1000-1600, 
and  1600-2200. 


calculations,  mixed-layer  depths  were  negative. 

The  strong  correlations  with  wind  speed  and  zonal  wind,  about  -0.51  and 
0.51,  respectively,  indicate  that  mechanical  energy  input  into  the  mixed  layer  by  wind 
forcing  plays  a  significant  role  in  mixing  the  surface  layer  at  0°,  140^.  The  sign  of  the 
correlations  indicate  that  increased  local  wind  stress  deepened  the  mixed  layer.  It  is 
interesting  to  point  out  that  over  larger  spatial  scales,  increased  easterlies  along  the 
equator  drive  Ekman  pumping,  which  would  raise  the  thermocline  and  result  in  shallower 
mixed-layer  depths.  Meridional  winds,  which  were  typically  weaker  than  the  zonal  winds 
(Figure  4.1),  were  significantly  correlated  with  the  day  and  evening  mixed-layer  depths, 
but  were  not  significantly  correlated  witl.<  the  night  and  morning  mixed-layer  depths.  This 
diurnal  difference  of  the  correlations  could  be  explained  by  the  diurnal  cycle  of 
meridional  winds  reported  by  Deser  (1994).  Her  analysis  of  the  diurnal  cycle  of  winds 
over  the  equatorial  Pacific  indicated  diurnal  meridional  wind  variations  of  about  0.6-0.8 
ms  '  at  most  locations,  with  southward  anomalies  of  the  flow  at  night.  Since  the 
meridional  winds  at  0°,  140°W  were  usually  northward  during  this  four-year  period  of 
study,  southward  anomalies  of  the  flow  at  night  indicate  a  deceleration,  which  could 
translate  into  decreased  nighttime  input  of  mechanical  energy  by  the  meridional  wind.  An 
alternative,  or  perhaps  complementary,  explanation  for  the  decreased  correlation  between 
the  meridional  winds  and  the  nighttime  mixed-layer  depth  could  be  that  the  weaker 
nighttime  meridional  winds  do  not  penetrate  to  the  deeper  nighttime  nuxed-layer  d^th. 
By  this  hypothesis,  the  daytime  mixed  layer  was  sufficiently  shallow  to  be  affected  by  the 
stronger  daytime  meridional  winds. 

Of  the  terms  of  the  surface  heat  flux,  only  latent  heat  flux,  and  the  net 
surface  heat  flux,  were  significantly  correlated  with  the  monthly  mean  mixed-layer 
depths.  The  positive  correlations  with  Q„,  and  G,  indicate  that  the  mixed-layer 
deepened  during  periods  with  increased  net  and  latent  heat  loss  from  the  ocean  surface. 
In  both  cases,  surface  heat  loss  drives  convective  mixing  and  mixed  layer  deqtening.  On 
these  monthly  time  scales,  Q,,  Q^,  and  Qi,  were  not  well  cmrelated  with  mixed-layer 
depths.  It  is  noteworthy  that  the  largest  term  of  the  surface  heat  flux.  G,.  was  not  well 
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correlated  with  mixed-layer  depth  on  monthly  time  scales.  In  section  4.3.2,  on  the  other 
hand,  good  correspondence  between  mixed-layer  depth  and  Q,  are  observed  on  seasonal 
time  scales. 


Table  4.1.  Correlation  coefficients,  r,  between  monthly  mean  mixed-layer  depths  for  each 
of  the  four  6-hour  time  bins  and  surface  forcing  variables  during  November  1983  - 
October  1987.  Statistical  significance  at  the  95%  confidence  level  was  calculated  using 
computed  integral  time  scales  between  time  series  to  determine  indq>endence  of  data 
points  (Davis,  1976).  Statistically  significant  ctxrelations  are  indicated  by  Bold  type. 


VARUBLE 

2200>0400 
95%c.i.  r 

0400^1000 
95%c.i.  r 

1000-1600 
95%c.i.  r 

1600-2200  1 
95%c.i.  r  1 

Wind  Speed 

.411 

•311 

.424 

•324 

399 

-.511 

.416 

-310  1 

Zonal  Wind 

.436 

325 

396 

.487 

.425 

mm 

Merid.  Wind 

.486 

-.377 

.501 

-.360 

Q3i 

AOA 

•305 

Wind  Dir. 

.467 

-.237 

.484 

-.226 

.438 

-.365 

.461 

-.352 

a 

.485 

-.100 

.491 

-.177 

.448 

.482 

-.057 

QJE 

369 

.438 

.454 

373 

.513 

Qh 

.419 

.156 

.429 

.184 

.412 

.078 

.420 

.098 

Qt 

.398 

.202 

.405 

.158 

.405 

.274 

.406 

.251 

.389 

.388 

.386 

.349 

J62 

.420 

.389 

.486 

P 

.316 

.066 

.308 

-.078 

.316 

.027 

£  -  P 

.330 

-.244 

.339 

-.262 

.322 

-.224 

.333 

-.318 

.320 

.134 

.330 

.158 

.315 

.005 

.318 

.106 

Bm. 

.319 

.275 

.317 

.272 

.314 

.286 

322 

.331 
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Comparison  between  mmithly  mean  mixed-layer  depth  (Figure  4.7a)  and 
monthly  mean  buoyancy  frequency  (Figure  4.6a)  shows  considerable  correspondence. 
When  the  high-A^  core  was  shallow,  such  as  during  the  1984  cold  episode,  the  mixed  layer 
was  likewise  shallow.  Similarly,  when  the  high-N  core  was  deep,  the  mixed  layer  was 
deep.  In  other  words,  the  correspondence  between  monthly  mean  mixed-layer  depths  and 
N  largely  reflects  the  fact  that  large-scale  vertical  displacements  of  the  thermocline  lead 
to  changes  in  mixed-layer  depth,  as  deftned  here,  through  mass  conservation,  an  adiabatic 
process  that  requires  no  mixing. 

The  correlation  between  monthly  mean  mixed-layer  depths  and  N  in  the 
upper  100  m  for  each  of  the  four  six-hour  time  bins  are  given  in  Table  4.2.  Statistically 
significant  positive  correlation  was  observed  between  mixed-layer  depth  and  N  at  the 
depths  between  35  m  and  80  m  for  all  four  time  bins,  and  at  25  m  and  100  m  during  the 
night  time  bin.  As  described  above,  the  positive  correlatimis  indicate  that  on  monthly 
time  scales  the  mixed-layer  depth  shallowed  when  the  stratification  increased.  It  is 
noteworthy  that  ctxrelations  were  consistently  stronger  during  the  night  and  mtxtiing  time 
bins,  perhaps  suggesting  a  more  direct  relationship  between  mixing  and  stratification  at 
times  when  active  convective  mixing  was  occurring  than  during  times  of  restratification. 

Significant  positive  correlations  were  observed  between  monthly  mean 
mixed-layer  depth  and  velocity  shear  at  25  m  and  45  m  during  all  four  time  bins,  and 
velocity  shear  at  35  m  during  the  day  time  bin  (Table  4.3  and  Figure  4.8).  This  depth 
range  of  significant  correlations  coincides  with  the  typical  depth  of  the  mixed  layer.  The 
positive  sign  of  the  ctNrelations  indicates  that  increased  shear  was  observed  during  times 
when  the  mixed-layer  was  relatively  shallow  and  vice  versa.  Furthermore,  ctxrelations 
in  this  dq}th  range  increased  during  the  day  time  bin,  when  the  mixed  layer  was 
diumally-shallowest  and  generally  above  this  depth  range.  These  results  imply  increased 
shear  at  the  base  of  the  mixed  layer,  when  the  mixed  layer  was  deep,  momentum  in  the 
mixed-layo-  was  well-mixed  and  shear  was  reduced.  When  the  mixed  layo-  was  shallow, 
momentum  was  not  well-mixed,  and  shear  was  high.  On  these  montttly  and  diurnal  time 
scales,  mixing  apparently  worked  to  deq>en  the  mixed  layer  and  decrease  the  shear. 
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Table  4J2.  Correlations  coefficients,  r.  between  monthly  mean  mixed'layer  depths  for 
the  four  6-hour  time  bins  and  buoyancy  frequency  in  the  upper  100  m  during  November 
1983-October  1987.  Statistical  signifrcance  at  the  95%  confidence  level  was  determined 
based  on  integral  time  scales  between  time  series  (Davis,  1976).  Statistically  significant 
correlations  are  indicated  by  Bold  type. 


VARIABLE 

2200-0400 

0400-1000 

1000-1600 

1600-2200  1 

95%c.i. 

r 

95%c.i. 

r 

95%c.i. 

r 

95%c.i. 

r 

N  at  010m 

.425 

.339 

.430 

.312 

.414 

.348 

.449 

.344 

N  at  025m 

.412 

.417 

.415 

.410 

.365 

.320 

.400 

.352 

N  at  035m 

J88 

.526 

J97 

MS 

J55 

.439 

J80 

.432  1 

1  N  at  04Sm 

.429 

J05 

.463 

J59 

.415 

.429 

.431 

.427  1 

.445 

J62 

.480 

.620 

.429 

.466 

.446 

.480  1 

1  fV  at  080m 

.417 

.641 

.464 

.688 

.420 

J60 

.619  1 

1  Af  at  100m 

.379 

.422 

.392 

.390 

.348 

.316 

J63 

.401  1 

Table  4J.  CcKrelations  coefficients,  r,  between  monthly  mean  mixed-layer  depths  for 
the  four  6-hour  time  bins  and  velocity  shear  in  the  uppier  1(X)  m  during  November  1983- 
October  1987.  Statistical  significance  at  the  95%  confidence  level  was  determined  based 
on  integral  time  scales  between  time  series  (Davis,  1976).  Statistically  signifrcant 
correlations  are  indicated  by  Bold  type. 


1  VARIABLE 

2200-0400 

0400-1000 

1000-1600 

1600-2200  1 

95%ci. 

r 

95%c.l. 

r 

95%cJ. 

r 

95%c.i. 

r 

S  at  010m 

.356 

.209 

.344 

.181 

.337 

.105 

.347 

.167 

S  at  025m 

J98 

.549 

.410 

J46 

mm 

.618 

397 

.478 

S  at  035m 

.357 

.269 

.346 

.255 

372 

.439 

.371 

.256 

S  at  045m 

.663 

.418 

.669 

.402 

.686 

396 

.634 

S  at  060m 

.400 

.210 

.423 

.237 

.407 

.272 

.410 

.263  1 

S  at  080m 

.450 

.059 

.457 

.064 

.419 

.089 

.432 

.122  1 

S  at  ICOm 

.480 

-.334 

.500 

-.411 

.450 

-.387 

.480 

-.330  1 
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Correlations  between  monthly  mean  mixed-layer  d^ths  and  Richardson 
number  in  the  upper  100  m  were  generally  weak  (Table  4.4  and  Figwe  4.8).  Statistically 
significant  correlations  were  found  only  with  Ri  in  the  upper  45  m,  and  only  for  qjecific 
time  periods.  At  10  m.  signiftcant  correlation  was  observed  only  during  the  night  time 
bin.  when  the  mixed  layer  was  deep.  At  25  m.  significant  correlations  were  observed 
during  the  night,  morning,  and  evening  time  bins,  with  the  largest  correlation.  -0.595. 
occurring  during  the  evening  time.  At  45  m,  significant  correlation  was  observed  only  for 
the  evening  time  bin.  The  sign  of  each  of  the  significant  cwrelations  was  negative 
indicating  that  the  mixed  layer  was  relatively  deep  during  times  when  Ri  in  the  upper  45 
m  was  high,  and  vice  versa.  This  reflects  the  fact  that  when  the  mixed  layer  was  deep, 
mixing  of  momentum,  and  to  a  lesser  extent  heat,  resulted  in  increased  Ri. 


Table  4.4.  Correlations  coefficients,  r.  between  monthly  mean  mixed-layer  dq>ths  for 
the  four  6-hour  time  bins  and  gradient  Richardson  Number  in  the  upper  100  m  during 
November  1983-C)ctober  1987.  Statistical  signiricance  at  the  95%  cmifidence  level  was 
determined  based  on  integral  time  scales  between  time  series  (Davis.  1976).  Statistically 
significant  correlations  are  indicated  by  Bold  type. 


1  VARUBLE 

2200-0400 

0400-1000 

1000-1600 

1600-2200  1 

95%c.i. 

r 

95%c.l. 

r 

95%c.l. 

r 

95%c.i. 

r 

Ri  at  lOm 

334 

.320 

-.317 

.317 

-.133 

.298 

-.181 

Ri  at  25in 

J91 

-.447 

-.499 

J99 

-395 

.396 

-.343 

Ri  at  35m 

.347 

-.216 

.362 

-.243 

377 

-.461 

.353 

-.243 

Ri  at  45m 

.308 

-.254 

.346 

-.336 

318 

-379 

.287 

-.164 

Ri  at  60m 

.381 

.070 

.398 

-.062 

.372 

-.047 

.398 

-.037 

Ri  at  80m 

.489 

.269 

.496 

.270 

.461 

.199 

.468 

.199 

Ri  at  100m 

.528 

.490 

542 

.523 

.499 

.445 

326 

.413 
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The  above  relationships  between  monthly  mean  mixed-layer  d^ths  and  N, 
S,  and  Ri  in  the  upper  100  m  are  summarized  in  Figure  4.8.  In  these  figures,  the  heavy 
green  dashed  lines  represent  the  four-year  mean  nighttime  and  daytime  mixed-layer 
depths.  The  green  dotted  lines  represent  the  uppermost  and  lowermost  standard  deviations 
about  the  respective  mean  mixed-layer  depths.  Figure  4.8a  graphically  shows  that  mixed- 
layer  depths  are  significantly  correlated  with  N  at  depths  weU  below  the  mixed  layer. 
These  correlations  below  the  base  of  the  mixed  layer  and  the  comparison  between 
monthly  mean  N  (Figure  4.6a)  and  mixed-layer  depths  (Figure  4.7a)  suggest  that  these 
correlations  are  probably  more  related  to  the  large-scale  dynamics  responsible  for  the 
vertical  migrations  of  the  thermal  structure  than  by  local  mixing  processes.  The 
relationships  with  S  and  Ri,  on  the  other  hand,  are  strong  only  near  the  base  of  the  mixed- 
layer,  suggesting  the  possibility  of  local  forcing. 


Correlations  Between  Mixed-Layer  Depth  and  N,  S,  k  Ri 


1.00-0.S0  0.00  0.80  1.00  -1.00-0.80  0.00  0.80  1.00  -1.00-0.80  0.00  0.80  1.00 
Correlation  Coetttclent  Correlation  Coefficient  Correlation  Coefficient 


Figure  4.8.  Qxrelation  coeffidents,  r,  between  monthly  mean  mixed-layer  depths  and 
buoyancy  frequency,  N,  velodty  shear,  5,  and  gradient  Richardson  number,  Ri,  in  the 
upper  100  m  during  November  i983-C)ctoba'  1987.  Dotted  black  vertical  lines  represmt 
95%  confidence  intervals  determined  using  an  integral  time  scale  between  time  series 
(Davis,  1976).  Dashed  heavy  green  lines  indicate  four-year  mean  daytime  (iq>per)  and 
nighttime  (lower)  mixed-layo*  dq>ths.  Dotted  green  lines  indicate  iq)per  and  lower  bounds 
of  the  standard  deviations  about  the  respective  mean  mixed-layer  depths. 
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C.  SEASONAL  CYCLES 

In  the  previous  section,  monthly  mean  time  series  of  several  surface  forcing  and 
oceanic  background  condition  variables  were  presented  and  discussed  in  terms  of  their 
interannual  variability.  In  several  of  those  discussions,  annual  cycles  were  observed  and 
in  some  instances  briefly  mentioned.  In  this  section,  seasonal  cycles  of  many  of  these 
variables  are  examined.  Generally,  seasonal  cycles  were  computed  as  means  of  monthly 
means  for  each  calendar  month  from  all  of  the  months  with  available  data  during  the 
period  November  1983  through  October  1987. 

1.  Surface  Forcing 

Climatological  seasonal  cycles  of  wind  speed,  zonal  wind,  and  meridional 
wind  are  shown  as  dotted  lines  in  Figure  4.1.  These  were  computed  from  all  of  the 
available  buoy  data  during  1983-1991  (Mcl^iaden  and  McCarty,  1992).  The  seasonal 
cycle  of  the  relatively  weak  meridional  winds  varied  from  a  minimum  near  zero  during 
February  and  March  to  a  maximum  of  about  1.5  ms'‘  during  August  and  September.  The 
seasonal  cycle  of  zonal  winds  ranges  from  a  minimum  easterly  of  about  -4.8  ms  '  in  May 
to  a  maximum  easterly  of  about  -6.8  ms  '  in  December  and  January.  The  seasonal  cycle 
of  wind  speed  is  dominated  by  the  zonal  component  with  minima  during  the  boreal  spring 
between  March  and  June  and  maxima  during  the  b(»eal  winter  between  October  and 
February.  Each  of  the  wind  components  and  wind  speed  show  a  weak  semi-annual  cycle 
indicated  by  a  weak  maximum  or  leveling-ofr  of  wind  speed  during  the  boreal  summer 
months  August  through  October. 

Climatological  seasonal  cycles  of  the  surface  heat  flux  terms  are  presented 
as  dotted  lines  in  Figures  4.2a-e.  The  Oberhuber  (1988)  climatology  was  used  for  the 
seasonal  cycles  of  (2a>  and  Q^.  Longwave  radiation  has  a  annual  maximum  magnitude 
of  about  -58  Wm'^  during  the  boreal  fall  between  October  through  December  and  a 
minimum  magnitude  of  about  -48  Wm'^  during  the  boreal  spring  between  March  and 
May.  The  climatological  seasonal  cycle  of  the  sensible  heat  flux  varies  by  only  about  3-4 
Wm'^.  The  estimates  of  Q,,  show  an  annual  minimum  of  about  -5  Wm'^  during  the  winter 
and  a  maximum  of  about  -i4  Wm'^  during  the  spring  and  summer.  These  seasonal  cycles 
of  (2a  are  probably  too  small  and  variable  to  play  a  significant  role  in  modulating  the 
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seasonal  cycles  of  convective  mixing  or  high-frequency  temperature  variability  in  the 
thermocline.  The  seasonal  cycle  of  latent  heat  flux,  rni  the  other  hand,  shows  a 
pronounced  annual  cycle  with  a  broad  maxima  of  about  -115  Wm  ^  during  the  boreal 
summer  and  fall  and  a  minimum  of  about  -75  Wm'^  from  February  through  April. 

The  seasonal  cycle  of  shortwave  radiation  is  presented  as  the  dotted  line 
in  Figure  4.2d.  This  was  computed  using  the  formula  of  Reed  (1977)  and  COADS  cloud 
data  and  buoy  data  for  the  period  1983-1992.  A  pronounced  semi-annual  cycle  was 
observed  in  response  to  the  semi-annual  cycle  of  the  ecliptic.  Asymmetric  shortwave 
radiation  minima  occur  about  the  time  of  the  solstices  in  Jime/July  and 
November/December  and  maxima  occur  during  the  equinoxes  in  March  and  Sq)tember. 
The  minimum  in  June/July  is  about  232  Wm'^  compared  with  the  November/December 
minimum  of  about  248  Wm'^.  The  September  maximum  reaches  about  273  Wm'^ 
compared  to  the  March  maximum  of  about  272  Wm'^  .  The  seasonal  cycle  of  the  data 
during  the  four-year  study  period  had  a  similar  pattern,  but  the  magnitudes  were  about  30 
Wm'^  greater. 

The  seasonal  cycle  of  the  net  surface  heat  flux  was  computed  from  the 
seasonal  cycles  of  the  four  component  terms  and  is  presented  as  the  dotted  line  in  Figure 
4.2e.  The  seasonal  cycle  of  Q^,  is  characterized  by  an  asymmetric  semi-annual  cycle  with 
minimum  values  of  54.1  and  83.3  Wm'^  occurring  in  July  and  December,  respectively, 
and  maximum  values  of  143.7  and  1(X).8  Wm'^  occurring  in  March  and  September, 
respectively.  The  annual  climatological  mean  is  99.1  Wm'^  and  the  standard  deviation  is 
25.6  Wm'^.  The  estimated  Q^,  had  a  significantly  greater  mean  of  154.5  Wm'^  and  a 
standard  deviation  of  29.7  Wm'^. 

The  seasonal  cycle  of  the  moisture  flux  (Figure  4.4)  is  characterized  by  a 
pronoimced  annual  cycle  with  a  maximum  during  February  and  March  and  a  minimum 
during  July  through  September.  This  annual  cycle  is  closely  related  to  the  seasonal  cycle 
of  precipitation  presented  in  Figure  3.8.  Since  the  net  surface  buoyancy  flux  was  strongly 
dominated  by  the  net  surface  heat  flux,  the  seastmal  cycle  of  is  almost  identical  to 
that  described  for  Q,^^. 
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2. 


Mixed  Layer 

The  seasonal  cycles  of  mixed-layer  dq}th  for  each  of  the  four  six-hour  time 
bins  and  the  diurnal  range  of  mixed-layer  depth  are  shown  in  Figure  4.9.  As  with  surface 
winds  and  shortwave  radiation,  both  mixed-layer  depth  and  the  diurnal  range  of  mixed- 
layer  depth  exhibited  a  weak  semi-annual  cycle  and  a  pronounced  annual  cycle.  Mixed- 
layer  depths  were  deepest  during  the  boreal  winter  months  December  and  January,  when 
the  nighttime  mixed  layer  deepened  to  about  50  m.  Mixed-layer  depths  were  shallowest 
during  the  boreal  spring  months  March  through  May.  when  the  nighttime  mixed  layer 
penetrated  to  a  depth  of  only  about  20  m.  For  the  period  between  June  and  November, 
nighttime  mixed-layer  depths  remained  relatively  stable  at  about  35  m.  The  seasonal 
cycles  of  mixed-layer  depth  and  surface  wind  speed  were  in  close  agreement,  with  deep 
mixed  layers  occurring  during  the  boreal  winter  months  when  wind  speeds  reached  their 
aimual  peak.  This  strong  relationship  between  increased  mixing  during  periods  of 
increased  wind  speed  are  not  unique  to  this  region  and  have  been  observed  and  modelled 
frequently.  The  semi-aimual  period  of  the  mixed  layer  corresponds  with  the  setni-armual 
period  of  Q„  with  a  shallow  mixed  layer  during  the  boreal  spring  and  fall  when  Q, 
reached  semi-annual  maximum  values.  This  relationship  implies  that  increased  shortwave 
radiation  increases  the  stratification  and  reduces  convective  mixing.  Since  Q„,  and 
include  this  semi-armual  cycle  of  Q^,  their  seasonal  cycles  also  correspond  with  die 
seasonal  cycle  of  the  mixed  layer. 

The  seasonal  cycle  of  the  diurnal  range  of  mixed-layer  depth  was  also 
characterized  by  a  semi-armual  cycle.  Maximum  diurnal  ranges  of  14-16  ra  occurred 
during  the  bweal  winter  months  October  through  February  and  again  during  May. 
Minimum  diurnal  ranges  of  7-11  m  occurred  during  March  and  April  and  again  during 
July  through  September.  The  large  diurnal  ranges  observed  during  the  winter  corresponds 
with  higher  winter  winds.  The  semi-annual  cycle  was  not  clearly  related  to  the  semi- 
armual  cycle  of  Q,. 
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SEASONAL  CYCLE  OF  MIXED-LAYER  DEPTH  ft  DIURNAL  RANGE 


Figure  4.9.  Seasonal  cycles  of  mixed-layer  depth  fw  four  six-hour  time  bins  and  diurnal 
range  of  mixed-layer  depth.  Computed  from  monthly  means  during  the  paiod  November 
1983-October  1987.  Black  votic^  lines  represent  ±one  mean  standard  error  of  the  mean 
computed  from  48  monthly  standard  errors  fw  the  mixed-layer  dq>th  and  12  standard 
OTors  fw  the  diurnal  range  of  mixed-layo'  depth. 
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3.  Buoyancy  Frequency 

The  seasonal  cycles  of  buoyancy  frequency.  N,  zonal  velocity,  U,  velocity 
shear.  5,  and  Richardson  number,  Ri,  over  the  depth  range  10-250  m  are  contoured  in 
Figure  4.10.  The  seasonal  cycle  of  N  in  the  upper  SO  m  (Figure  4.10a)  reflected  a  relation 
to  the  seasonal  cycle  of  the  mixed-layer  depth  (Figure  4.9).  During  the  boreal  spring, 
when  the  mixed  layer  was  shallow,  N  increased  to  above  0.01  s  '  in  the  upper  25  m  and 
above  0.012  s  '  between  25  m  and  50  m.  For  the  remainder  of  the  year,  when  the 
nighttime  mixed  layer  extended  down  to  at  least  35  m,  N  in  the  upper  35  m  or  so 
remained  less  than  0.01  s  '  and  between  May  and  September  N  was  less  than  0.008  s  '  in 
the  upper  25  m. 

The  high-N  core  associated  with  the  depth  of  the  thermocline  migrated 
vertically  during  the  year  from  a  maximum  depth  during  the  btu-eal  winter  of  about  125 
m  to  a  minimum  dq>th  of  about  100  m  during  the  boreal  spring.  The  magnitude  of  the 
high-N  c(ve  varied  during  the  year  from  about  0.015  s  '  in  January  to  over  0.018  s*' 
from  June  through  September.  The  vertical  gradient  of  A/  underwent  noteworthy  seasonal 
variations  with  strong  gradients  observed  between  the  mixed  layer  and  thermocline  core 
(25-75  m)  during  the  period  May  through  September  and  weak  votical  gradients  during 
the  boreal  winter.  These  weaker  vertical  gradients  of  Af  during  the  winter  ctxrespond  with 
the  deq)ening  of  the  mixed  layer. 

4.  Zonal  Velocity 

The  seasonal  cycle  of  zonal  velocity,  I/,  over  the  depth  range  10-250  m  is 
contoured  in  Hgure  4.10b.  The  most  striking  feature  of  the  seasonal  cycle  of  £/  is  the 
annual  vertical  migration  and  intensifreation  of  the  EUC.  During  bweal  spring  and  early 
summer,  the  eastward  EUC  core  intensified  in  strength  to  over  1.2  ms''  and  shoaled  to 
typical  depths  of  about  85  m.  This  shoaling  of  the  EUC  often  extended  to  the  surface, 
where  weak  eastward  cunents  existed  during  May  and  June.  By  July,  tlw  EUC  began  to 
deepen  and  reduce  in  strength.  During  the  boreal  winter,  the  EUC  core  was  generally 
centered  at  a  dq>th  of  about  125  m  and  had  a  reduced  strength  of  less  than  1.0  ms  '.  From 
July  through  March,  the  surface  velocities  wo-e  westward  at  mean  speeds  between  0.2  and 
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0.4  ms  '.  This  westward  surface  flow  was  part  of  the  large-scale  SEC  driven  by  the 
easterly  winds  at  the  surface. 

5.  Velocity  Shear 

The  seasonal  cycle  of  velocity  shear  over  the  depth  range  10-250  m  is 
contoured  in  Figure  4.10c.  This  figure  shows  ctxrespondence  with  the  zonal  velocity 
contour  plot,  eventhough  shear  is  based  on  both  zonal  and  meridional  velocity.  As 
observed  fw  zonal  velocity,  shear  underwent  noteworthy  seasonal  variations.  The  high- 
shear  core  migrated  vertically  from  about  80-90  m  during  the  boreal  fall  and  winter 
months  to  about  35-45  m  during  the  Ixxeal  spring.  While  shear  between  the  surface  and 
the  EUC  core  was  always  relatively  strcmg  with  values  exceeding  -0.012  s  '  throughout 
the  year,  maximum  shears  exceeding  -0.020  s  '  occurred  during  the  boreal  summer 
months  when  the  EUC  core  was  shallow,  but  beginning  to  deqpen,  and  the  strength  was 
intensified.  In  April  and  May,  when  the  EUC  was  strongest  and  extended  to  the  surface, 
maximum  shears  of  -0.020  s  '  occurred  at  about  35  m.  During  June  through  August,  when 
the  EUC  had  begun  to  deepen  and  the  SEC  surface  westward  flow  had  returned, 
maximum  shears  of  -0.024  s''  occurred  at  about  60  m.  During  the  boreal  winter  months, 
when  easterly  winds  were  strcmgest,  the  mixed  layer  was  deqjest.  and  the  EUC  was 
weaker  and  deeper,  the  shear  was  about  -0.012  s''  from  about  30  m  to  about  100  m. 

Below  the  EUC  core,  maximum  shears  of  over  0.012  occurred  from  May 
through  October  at  a  depth  of  about  140  m.  During  the  boreal  winter  months,  the  shear 
below  the  EUC  was  reduced  to  values  near  0.004  s''. 

The  relatively  weak  shear  observed  below  the  EUC  core  might  have  been 
partially  reduced  by  the  poOT  vertical  lesolutirm  of  current  meters  at  these  depths.  The 
comparison  between  shears  cmnputed  from  ADCP  measurements  with  8  m  vertical 
resolution  and  those  used  here  which  was  di«;ussed  in  section  3.2.3  showed  considerable 
underestimates  of  shear  at  depths  between  140  m  and  180  m. 

6.  Richardson  Number 

The  seasonal  cycle  of  Ri  over  the  dq)th  riqge  10-250  m  is  contoured  in 
Figure  4.10d.  As  discussed  in  section  4.1.2,  the  averaging  operators  used  in  confuting 
monthly  mean  Ri  have  the  effect  of  increasing  low-Ri  values.  In  computing  seasonal 
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cycles  of  Ri,  this  effective  smoothing  was  repeated  and  mean  seasonal  values  in  low-/?/ 
regions  were  increased.  The  seasonal  cycle  of  Ri  shows  the  expected  similarity  with  the 
seasonal  cycle  of  shear.  Minimum  values  of  Ri  coincided  with  maximum  values  of  shear. 

Druing  the  period  March  through  June,  when  the  mixed-layer  was  shallow,  mean  Ri  at 
35  m  was  less  than  unity,  with  minimum  values  near  0.5.  This  low-/?/  region  deepened 
to  about  60  m  from  July  through  October  with  mean  values  of  0.5  for  much  of  this  time. 
This  deepening  coincided  with  a  deepening  of  the  seasonal  mixed-layer  of  about  20  m. 
E>uring  the  periods  when  Ri  was  near  0.5,  the  vertical  extent  of  the  low-/?/  core  was 
noticeably  reduced,  typically  to  depth  ranges  of  about  25  m.  During  the  bweal  winter 
months,  when  mean  Ri  in  the  low-/?/  core  was  typically  slightly  greater  than  unity,  the 
vertical  extent  of  the  low-/?/  ccxe  appeared  greater.  The  high-/?/  core  appeared  to  follow 
the  low-shear  contour  of  Figure  4.10c.  The  low-/?/  core  below  the  EUC  core  had  a 
minimum  of  about  2.0  centered  at  about  140  m  during  the  months  ^>ril  through  July. 
These  relatively  high  values  of  Ri  in  the  low-/?/  core  were  probably  overestimates  due  to 
the  likely  underestimates  of  shear  at  these  depths. 
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Figures  4.10.  Tiine-dq)th  contour  plots  of  seasonal  cycles  of  a)  buoyancy  firequency, 

N,  b)  zonal  velocity,  U,  c)  velocity  shear,  S,  and  d)  Richardson  numbw,  Ri,  at  (f , 
14(FW  over  the  dq>th  range  10-250  m.  Computed  from  monthly  mean  values  during 
November  1983-October  1987.  Contour  intairals  are  a)  0.2  x  10^  s‘\  b)  20  ons  c) 

O. 25  X  10*  s*‘,  and  d)  0.2. 
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V.  ISOTHERM  DISPLACEMENT  VARIANCE 


A.  ISOTHERM  DISPLACEMENT  VARIANCE 

The  primary  objective  of  this  study  is  to  characterize  the  diurnal  variability  of 
high-frequency  wavelike  motions.  The  indicattM*  high-frequency  motions  used  in  this 
investigation  is  isotherm  displacement  variance  (IDV),  which  is  interpreted  as  some 
combination  of  high-frequency  internal  waves,  local  Kelvin-Helmholtz  shear  instabilities, 
and  local  turbulence.  In  the  following  five  sections,  IDV  will  be  discussed  in  terms  of  its 
interannual  and  seasonal  variability,  and  the  relatitmships  with  the  observed  surface 
forcing  and  dynamic  stability. 

B.  INTERANNUAL  VARIABILITY  OF  IDV 

This  section  focuses  on  the  interannual  variability  of  IDV  during  the  period 
November  1983  through  October  1987.  Mtmtiily  mean  values  of  IDV  at  the  dq>ths  35  m. 
60  m,  100  m,  140  m  and  200  m  for  each  of  the  four  six-hour  time  bins  are  presented  in 
Figures  5.1a-e.  IDV  at  35  m  varied  substantially  over  interannual  time  scales  (Figure 
5.1a).  During  the  1984  ENSO  cold  ^isode,  IDV  at  35  m  was  low,  with  values  avoaging 
about  2-3  m^,  cwresponding  to  vertical  isotherm  displacements  of  about  1.5  m.  B^inning 
in  the  fall  of  1984,  a  trend  of  increasing  IDV  at  35  m  was  observed,  which  was 
particularly  ^^arent  during  the  night  and  morning  time  bins.  By  the  winter  of  1985/1986, 
typical  values  of  IDV  at  35  m  averaged  about  7-8  m^  corresptmding  to  vertical  isotherm 
displacements  of  about  2.5  m.  This  trend  was  superimposed  over  an  annual  cycle  witii 
minima  typically  occurring  during  the  boreal  sfxmg.  Following  the  annual  minimum  in 
the  spring  of  1986,  IDV  at  35  m  increased  sharply  to  values  as  high  as  17  m^  during  die 
boreal  winter  months,  correspondiiig  to  vertical  isotherm  di^lacements  exceeding  4  m. 
The  nighttime  increases  of  IDV  at  35  m  were  particularly  pronounced  during  periods  of 
relatively  high  IDV  and  diminished  during  periods  of  low  IDV.  It  is  noteworthy  tiiat 
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periods  of  low  IDV  at  35  m  occurred  when  stratification  at  35  m  was  strong,  and  vice 
versa. 

Due  to  three  large  data  gi^s  of  5-6  months,  the  time  series  of  IDV  at  60  m  is 
not  ideally  suited  to  investigating  interannual  (m*  seasonal  variability  (Figure  5.1b). 
Nevertheless,  significant  variations  on  these  time  scales  were  q^arent.  Typical  values  of 
IDV  at  60  m  were  about  4  m^  with  vertical  isotherm  displacements  of  about  2  m.  There 
was  a  noteworthy  reduction  in  IDV  at  60  m  during  the  later  half  of  the  1986/1987 
ENSO  warm  episode.  Between  February  and  October  1987,  IDV  at  60  m  averaged  less 
than  2  m^.  During  this  period  and  during  the  winter  of  1984  diurnal  variations  of  IDV  at 
60  m  were  not  observed. 

At  1(X)  m,  which  was  typically  near  the  high-V  core  of  the  thermocline,  IDV 
amplitudes  were  greater  than  at  35  m,  60  m,  140  m,  and  2(X)  m  (Figure  5.1c).  Only  at  300 
m  was  IDV  greater  than  at  1(X)  m,  possibly  due  to  the  required  use  of  first  differences, 
rather  than  centered  differences,  in  computing  the  vertical  ten^rature  gradient  used  to 
normalize  the  temperature  variances.  Interannual  variations  of  IDV  at  100  m  are  iq>parent. 
During  the  1984  cold  episode,  IDV  at  100  m  was  high,  with  an  average  of  about  10  m^. 
During  this  period,  IDV  at  100  m  exhibited  strong  oscillatioins  having  a  period  of  about 
4  months  and  an  amplitude  of  about  5  m^.  Throughout  1985  and  into  the  spring  of  1986, 
the  average  IDV  decreased  to  about  5-6  m^  with  oscillations  having  a  period  of  about  5-7 
months  and  weaker  amplitudes  of  about  2  m^.  During  the  boreal  winter  of  1986/1987, 
IDV  at  1(X)  m  increased  to  about  10  ml  Generally,  diurnal  variations  of  IDV  at  100  m 
were  less  fq)parent. 

At  140  m,  which  was  typically  below  the  high-velocity  core  of  the  BUG  and 
below  the  hi^-N  core  of  the  thermocline,  magnitudes  of  IDV  were  less  than  at  100  m, 
with  values  averaging  about  6  m^  (Figure  5.1d).  An  exception  occurred  during  the  1984 
ENSO  cold  episode,  when  IDV  at  140  m  averaged  only  1.5  m^.  Aside  from  this  period 
of  low  IDV,  interannual  variations  of  IDV  at  140  m  were  not  jx’onounced.  No  consistent 
pattern  of  diurnal  variatitxis  of  IDV  at  140  m  were  observed. 
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Isotherm  Displacement  Variance  at  0°,140^ 
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Figure  S.la-e.  Monthly  mean  isotherm  displacement  variance  (IDV)  at  5  dq)ths  for  the 
period  November  1983-October  1987  at  0",  140^:  a)  35  m,  b)  60  m,  c)  100  m,  d)  140 
m,  and  e)  200  m.  Gq>s  in  each  record  are  periods  of  no  data  or  data  of  questionable 
quality.  Erro'  bars  indicate  the  mean  standard  enror  of  the  mean  for  all  time  bins  over  the 
entire  record  length.  Units  are  m^. 
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At  200  m,  where  both  the  strength  of  S  and  N  were  typically  weak  (Figure  5.1e), 
IDV  exhibited  subtle  variability  on  interannual  time  scales  not  easily  identified  as  ENSO- 
related.  For  the  first  18  months  of  the  reoxd,  IDV  at  2(X)  m  displayed  a  weak  downward 
trend  with  an  average  IDV  of  about  8  m^.  Following  a  seven-month  data  gap,  the  Enal 
two  years  of  data  reveal  a  noteworthy  oscillation  of  IDV  at  2(X)  m  with  an  annual  period. 
The  average  values  of  IDV  during  this  period  varied  between  about  4  m^  during  the 
bOTeal  winters  and  about  6  m^  during  the  boreal  springs.  No  consistent  pattern  of  diurnal 
variability  of  IDV  at  2(X)  m  was  observed. 

Although  IDV  was  computed  at  3(X)  m,  lack  of  confidence  in  the  vertical 
temperature  gradient  calculation  and  the  strong  sensitivity  of  IDV  to  this  calculation 
rendered  it  questionable.  Weak  temperature  gradients  produced  large  monthly  mean  IDV 
values  in  the  range  15-40  m^.  These  large  IDV  values  were  accompanied  by  large 
standard  errors.  No  consistent  pattern  of  diurnal  variability  of  IDV  at  3(X)  m  was 
observed. 

To  summarize  the  vertical  structure  of  IDV  over  the  four-year  study  period.  Figure 
5.2  shows  the  mean  vertical  {XoEIes  of  IDV  for  each  of  the  four  six-hour  time  bins 
computed  from  the  monthly  mean  IDV  values  during  November  1983  through  October 
1987.  The  mean  IDV  at  35  m  and  60  m  is  highest  during  the  mtxning  time  bin.  The 
mean  IDV  at  35  m  is  lowest  during  the  evening  time  bin  and  the  mean  IDV  at  60  m  is 
lowest  during  the  day  time  bin.  Generally,  IDV  in  Ae  upper  60  m  is  higher  during  the 
nighttime  and  lower  during  the  daytime.  A  striking  feature  of  the  mean  votical  structure 
of  IDV  is  the  high  variability  of  the  mean  at  35  m  and  60  m  compared  with  the  low 
variability  of  the  mean  IDV  at  100  m,  140  m,  and  2(X)  m. 
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Figure  5.2.  Mean  vertical  profiles  of  isothenn  displacment  variance  (IDV)  for  four 
six-hour  time  bins  over  the  d^th  range  35-200  m.  Mean  profites  computed  firom  monthly 
means  of  IDV  during  November  1983  -  October  1987.  Small  arrows  indicate  tone 
standard  error  of  the  mean. 
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C.  SEASONAL  CYCLE  OF  IDV 


The  seasonal  cycles  of  IDV  at  35  m.  60  m,  100  m,  140  m.  and  200  m  for  each 
of  the  four  six-hour  time  bins  are  presented  in  Figures  5.3a-e.  These  seasonal  cycles  were 
computed  from  the  monthly  mean  values  of  IDV  presented  in  Figures  5.1a-e.  The 
seasonal  cycle  of  IDV  at  35  m  shows  a  pronounced  annual  cycle  with  nighttime  maxima 
of  about  10  m^  occurring  during  the  boreal  winuo*  months  November  through  January  and 
nighttime  minima  of  about  2  m^  occurring  during  the  boreal  spring  months  March  and 
April  (Figure  5.3a).  A  pronounced  seasonal  cycle  of  the  diurnal  cycle  of  IDV,  AIDV,  at 
35  m  was  well  correlated  with  the  seasonal  cycle  of  IDV.  Maxima  of  AIDV  at  35  m  of 
about  6  m^  occurred  during  the  months  November  through  January  and  minima  near  zero 
occurred  during  the  months  March  and  A^jril. 

The  seasonal  cycle  of  IDV  at  60  m  had  a  pronounced  annual  cycle  of  the 
nighttime  IDV  with  maxima  of  about  5.2  m^  occurring  during  the  boreal  winter  months 
November  through  January  and  minima  of  less  than  2.0  m^  during  May  and  June  (Figure 
5.3b).  The  pattern  was  similar  to  that  described  for  IDV  at  35  m  excq)t  the  magnitudes 
were  reduced  and  there  appears  to  have  been  a  higher  frequency  seasonal  osdllation  with 
a  period  of  3-4  months.  The  seasonal  cycle  of  MDV  at  60  m  was  also  similar  to  that  at 
35  m  with  maxima  of  about  2.5  m^  occurring  from  November  through  January  and 
minima  near  zero  in  May  and  June. 

The  seasonal  cycle  of  IDV  at  100  m  underwent  iq>parent  annual  and  bi-monthly 
variations  (Figure  5.3c).  Annual  maxima  exceeding  10  m^  occurred  in  October  and 
January  and  an  annual  minimum  of  about  4.5  m^  occurred  in  August  There  does  not 
appear  to  have  been  a  significant  seasonal  cycle  of  AIDV  at  100  m.  There  was,  however, 
a  significant  seasonal  peak  in  AIDV  during  the  month  of  October.  The  seasonal  cycles 
of  IDV  and  AIDV  at  140  m  and  200  m  were  not  significant  (Figures  5.3d,e). 
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Figure  5.3a>e.  Seasonal  cycles  of  isotherm  displacement  variance  (IDV)  for  the  four 
six-hour  time  bins  at  five  depths:  a)  35  m,  b)  60  m,  c)  100  m,  d)  140  m,  and  e)  200 
m.  Seasonal  cycles  are  computed  from  monthly  mean  values  of  IDV  during  the  period 
November  1983-October  1987.  Error  bars  indicate  the  mean  standard  erro'  of  the  mean 
fn*  all  time  bins  over  the  twelve  calendar  months.  Units  are  m^. 
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D.  RELATIONSHIPS  BETWEEN  IDV  AND  SURFACE  FORCING 


In  this  section,  statistical  relationships  between  monthly  means  of  surfia^e  fm'cing 
conditions  and  monthly  means  of  IDV  at  3S  m,  60  m,  100  m,  140  m,  200  m,  and  300  m 
are  examined.  Table  5.1  shows  correlation  coefficients  computed  between  IDV  at  these 
six  depths  and  surface  winds.  Wind  speed  and  zonal  wind  were  significantly  correlated, 
at  the  95%  conAdence  level,  with  IDV  at  35  m  during  the  evening  and  night  time  bins 
and  with  IDV  at  60  m  during  the  night  and  morning  time  bins.  These  time  differences 
of  the  significant  correlations  suggest  a  possible  phase  lag  with  depth  of  about  6  hours 
that  is  coa  tent  with  phase  lags  of  turbulent  dissipation  observed  during  TROPIC  HEAT 
(Peter  •it  al.,  1988;  Mourn  et  al.,  1989)  and  modelled  numerically  (Garwood  et  al.,  1989; 
Schudlich  and  Price,  1992;  Skyllings^  and  Denbo,  1994).  The  lower  correlations 
occurred  after  sunrise,  when  the  surface  layer  underwent  restratification.  The  significant 
correlations  at  35  m  and  60  m  were  positive  for  wind  speed  and  negative  for  the  easto'ly 
(negative)  zonal  wind.  In  each  case,  IDV  increased  with  increased  mechanical  energy 
input  by  the  wind.  Surface  winds  were  not  generally  well  correlated  with  IDV  at  100  m, 
140  m,  200  m.  and  300  m.  with  two  exceptions.  Significant  correlations  were  observed 
between  both  wind  speed  and  zonal  wind  and  IDV  at  200  m  during  the  evening  time  bin 
and  with  IDV  at  3(X3  m  during  the  morning  time  bin.  The  lack  of  significant  correlation 
between  surface  winds  and  IDV  at  intermediate,  depths  argues  against  local  wind 
forcing  and  vertical  wave  propagation  to  explain  the  deqier  correlations.  Mwe  likely, 
they  are  related  to  larger-scale  processes.  Meridional  winds  and  wind  direction  were  not 
well  correlated  with  IDV  at  any  of  the  six  depths  examined  for  any  of  the  four  six-hour 
time  bins. 

Statistical  relationships  were  examined  between  IDV  at  the  six  TR  depths  for  each 
of  the  four  six-hour  time  bins  and  each  of  the  terms  of  the  surface  heat  and  buoyancy 
balances.  Correlations  were  computed  ftx*  the  following  terms:  net  surface  heat  flux  (Q^X 
shortwave  radiation  (QJ,  latent  heat  flux  (QX  sensible  heat  flux  (QX  longwave  radiation 
(QX  surface  buoyancy  flux  (B„J,  moisture  flux  (»iec4>itation  (F),  and 

evaporation  minus  precipitation  (E-F).  Genonlly,  the  surface  fluxes  wo’e  not  well 
cmrelated  with  IDV  at  most  d^ths.  The  following  terms  were  found  to  be  not 
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Table  5.1.  Correlations  coefficients,  r,  between  monthly  mean  isotherm  displacement 
variances  (IDV)  at  six  depths  for  four  6-hour  time  bins  and  surface  winds  (WSPD;  wind 
speed;  U:  zonal  wind;  V:  meridional  wind;  WDIR:  wind  direction)  during  November  1983 
-  October  1987.  Statistical  significance  at  the  95%  confidence  level  was  calculated  using 
computed  integral  time  scales  between  time  series  (Davis,  1976).  Statistically  significant 
correlations  are  indicated  by  Bold  type. 

VARIABLES  2200-0400  0400-1000  1000-1600  1600-2200 

95%c.i.  r  95%c.i.  r  95%c.i.  r  95%c.i.  r 


WSPD/IDV  35 

.418 

.477 

.339 

.241 

.309 

.261 

■Em 

Jll 

J35 

.339 

.346 

.240 

.176 

.263 

.145 

/IDVIOO 

.310 

.009 

.287 

-.163 

.300 

-.063 

.309 

-.108 

/IDV140 

.311 

.011 

.418 

.244 

.360 

.113 

.388 

.043 

/IDV200 

.399 

.263 

.271 

-233 

IQ9 

-.087 

J47 

-291 

/IDV300 

.282 

-.270 

246 

267 

.070 

.283 

.042 

U  /IDV  35 

.430 

-.492 

.341 

-.256 

-263 

J79 

-.403 

J15 

•J46 

.341 

-J59 

.246 

-.178 

.263 

-.144 

/IDViOO 

.311 

.011 

.287 

.189 

.299 

.094 

.311 

.136 

/IDV140 

.306 

.015 

.423 

-279 

.361 

-.121 

/IDV200 

.401 

-294 

.270 

.226 

129 

.074 

252 

289 

/IDV300 

.281 

.263 

236 

-253 

.250 

-.077 

.278 

-.051 

.486 

.202 

252 

.145 

.258 

.121 

.432 

.304 

/IDV  60 

.073 

.264  .096 

.278 

.114 

/IDVIOO 

.302 

.155 

.307 

.152 

.311 

.181 

.204 

/IDV140 

.312 

.153 

.341 

-.176 

.247 

.019 

.307 

-263 

/IDV200 

.359 

-.193 

.239 

-.064 

.293 

.022 

.391 

1 

00 

/IDV300 

.344 

1 

00 

.281 

.162 

.302 

-.042 

WDIR/IDV35 

.486 

.183 

IQ 

.174 

.253 

.079 

^99BS 

/IDV  60 

.267 

-.027 

IQH 

.032 

.272 

.076 

/IDVIOO 

.302 

.061 

.291 

.100 

.295 

.074 

.310 

.128 

/IDV140 

.076 

.341 

-.144 

.266 

.118 

.294 

-.180 

/IDV200 

.339 

-.246 

.229 

.293 

.152 

.393 

-203 

/IDV300 

.337 

-.219 

.264 

.093 

.267 

-.158 

.293 

.038 

95 


significantly  correlated  with  IDV  at  any  of  the  six  TR  depths  for  any  of  the  time  bins  and 

will  not  be  specifically  discussed:  Qj,,  Bn-. _ P,  and  E-P.  Correlations  between  monthly 

mean  IDV  and  remaining  terms  of  the  surface  heat  balance  are  given  in  Table  5.2.  and 
the  correlations  between  IDV  and  the  remaining  terms  of  the  surface  buoyancy  balance 
are  given  in  Table  5.3. 

As  shown  in  Table  5.2,  correlations  with  the  terms  of  the  surface  heat  flux  were 
generally  small.  The  only  term  having  statistically  significant  correlation  at  the  95% 
confidence  level  with  IDV  in  the  upper  100  m  was  Q,.  which  had  a  ccwrelation  of  -0.422 
with  IDV  at  35  m  during  the  evening  time  bin.  The  negative  sign  of  this  correlation 
indicates  that  IDV  was  high  when  heat  loss  from  the  ocean  by  (2,  increased.  The  latent 
heat  flux  was  also  significantly  correlated  with  IDV  at  200  m  (r  =  0.310).  Correlations 
between  IDV  at  140  m  and  Q,,  and  Q„,  were  significant  during  the  morning  and  day  time 
bins  and  day  and  evening  time  bins,  respectively.  It  does  not  seem  likely  that  the  small 
changes  in  Q,,  could  be  strong  enough  to  locally  force  temperature  variability  at  140  m, 
particularly  since  (2*  was  not  well  ccMrelated  with  IDV  in  the  upper  100  m.  Although  Q„, 
could  possibly  generate  vertically-propagating  waves,  the  lack  of  correlation  in  the  upper 
100  m  suggests  non-local  fcHTcing. 

The  cmrelations  between  IDV  at  the  six  TR  dq>ths  for  the  four  six-hour  time  bins 
and  B..,  and  E  were  also  generally  weak  (Table  53).  Both  Bg^  and  E  were  significantly 
correlated  with  IDV  at  35  m  during  the  evoiing  time  bin.  The  correlation  with  B.^  was  - 
0.374,  indicating  that  IDV  was  higher  when  buoyancy  decreased.  Similarly,  the  correlation 
with  E  of  0.421  indicates  that  IDV  was  high«'  when  evaporation  increased. 
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Table  5.2.  Correlation  coefficients,  r,  between  monthly  mean  isotherm  displacement 
variances  (IDV)  at  six  depths  fcx*  four  time  bins  and  surface  heat  flux  terms;  net  surface 
heat  flux,  shortwave  radiation.  Q,^  latent  heat  flux.  Q^,  and  sensible  heat  flux.  Q„. 
Statistical  signiflcance  at  the  95%  confldence  level  was  calculated  based  on  the  integral 
time  scale  between  time  series  (Davis.  1976).  Statistically  significant  correlations  are 
indicated  by  Bold  type. 


VARIABLES 

2200-0400 
95%c.i.  r 

0400-1000 
95%c.i.  r 

1000-1600 
95%c.i.  r 

1600-2200 
95%c.i.  r 

QJTD\  35oi 

.372 

-.195 

.368 

-.070 

.316 

-.108 

.357 

-.339 

.300 

-.003 

19 

.289 

-.024 

/IDVlOOm 

.317 

.027 

BSB 

.120 

.303 

.038 

/IDV140in 

.317 

.027 

.340 

.292 

Jll 

J45 

206 

259 

/mVlQOin 

.321 

.100 

.273 

.238 

.279 

.041 

-.045 

/IDV300I1I 

.255 

-.038 

.235 

-.193 

.273 

.121 

Q.  /IDV  35ni 

.255 

.391 

.058 

.325 

.017 

.438 

.173 

/IDV  60m 

.299 

.116 

.318 

.174 

.281 

-.007 

.278 

-.001 

/IDVlOOm 

.327 

.019 

.311 

.112 

.315 

.019 

.326 

.044 

/IDV140m 

.325 

.021 

.378 

.330 

.313 

.169 

.303 

.242 

/IDVlOOm 

-.051 

iOl 

fm 

268 

-499 

/IDV300m 

J61 

-JOl 

211 

.059 

.221 

.013 

.285 

.237 

Q,nD\  35m 

.393 

-.307 

wsm 

294 

-.093 

257 

-422 

/IDV  60m 

-.073 

.337 

-.071 

244 

-.014 

269 

-.043 

/IDVlOOm 

.326 

.017 

280 

.284 

-.037 

.317 

.016 

/IDV140m 

.317 

.013 

.406 

.146 

.370 

297 

.249 

/IDVlOOm 

.395 

286 

JIO 

.107 

.356 

.270 

/IDV300m 

.288 

.190 

-245 

.244 

-.118 

ES 

QJtD\  35m 

.426 

-.317 

.328 

-.186 

.317 

-.033 

.391 

-.105  1 

/IDV  60m 

.276 

-206 

.304 

-224 

.230 

-.033 

.271 

/IDVlOOm 

.304 

.149 

.284 

-.033 

.288 

.159 

.311 

/IDV140m 

.306 

.157 

J72 

-A32 

216 

•263 

.342 

-262 

/IDVlOOm 

.347 

-253 

.263 

-.196 

.292 

-.158 

.367 

.045 

/IDV300m 

.293 

-.063 

248 

-.133 

.254 

-.121 

WM_ 

-273 
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Table  5.3.  Ccnrelation  coefficients,  r,  between  monthly  mean  isotherm  displacement 
variances  (IDV)  at  six  dq>ths  for  four  time  bins  and  the  net  surface  buoyancy  flux,  B.,. 
Statistical  significance  at  the  95%  confidence  level  was  calculated  based  rni  the  integral 
time  scale  between  the  time  series  (Davis.  1976).  Statistically  significant  correlations  are 
indicated  by  Bold  type. 


VARIABLES 

2200-0400 
95%c.i.  r 

0400-1000 
95%c.i.  r 

1000-1600 
95%c.i.  r 

IW  * 

95%c.i. 

B,^  /IDV35in 

.280 

-.123 

11^1 

BBS 

J87 

-.313  1 

.295 

-.057 

.309 

-.131 

.275 

.108 

.286 

.163 

/IDVlOOm 

.329 

.139 

.326 

.131 

.321 

.054 

.301 

.050 

/IDVldOm 

.329 

.146 

.303 

.079 

.317 

.149 

.320 

.220 

/IDV200m 

.303 

-.121 

.302 

.086 

.285 

.218 

.266 

.023 

/IDVJOOm 

.272 

-.001 

.276 

-.191 

.291 

-.191 

.283 

.107 
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E.  RELATIONSHIP  BETWEEN  IDV  AND  DYNAMIC  STABILITY 


In  this  section,  statistical  relationships  between  monthly  means  of  the  dynamic 
stability  parameters,  N,  S,  and  Ri,  and  IDV  at  35  m,  60  m,  and  100  m  are  discussed. 
Correlation  coefficients  for  each  of  the  above  parameters  and  IDV  at  these  three  dq>ths 
for  the  four  six-hour  time  bins  are  presented  in  Table  5.4.  Generally,  correlations  between 
IDV  and  these  stability  parameters  were  weak.  Monthly  mean  IDV  was  candy 
correlated  at  the  95%  confidence  level  with  N  at  35  m  during  the  day,  evening  and  ni^t 
time  bins.  These  correlations  were  negative  indicating  that  IDV  at  35  m  increased  with 
decreasing  stratification.  This  result  is  consistent  with  linear  theory  in  that  increased 
isotherm  displacements  are  expected  with  reduced  stratification.  Correlations  between 
IDV  and  N  at  60  m  and  100  m  were  not  significant  at  the  95%  confidence  level. 
Correlations  between  both  IDV  and  S  and  IDV  and  Ri  were  not  statistically  significant 
at  35  m,  60  m,  or  100  m. 


Table  5.4.  CtMrelation  coefficients,  r,  between  mtmthly  mean  isotherm  displacement 
variances  (IDV)  and  buoyancy  frequency  (AO,  velocity  shear  (5),  and  Richardson  number 
(Ri)  at  35  m,  60  m,  and  100  m  fix  four  six-hour  time  bins  during  November  1983  - 
October  1987.  Statistical  significance  at  the  95%  confidence  level  was  calculated  based 
on  the  integral  time  scale  between  the  time  series  (Davis,  1976).  Statistically  significant 
ctxrelations  are  indicated  by  Bold  type. 


1  Variables 

22004)400 
95%ci.  r 

0400-1000 
!)5%c.i.  r 

1000-1600 
95%c.i.  r 

1600-2200  1 
95%ci.  r  1 

1  IDV/)V_35in 

J82 

•M4 

.358 

-.356 

321 

-351 

339 

-.394 

1  60m 

.323 

-.194 

.356 

-.323 

.219 

3i9 

.252 

-.162 

1  100m 

.291 

-.161 

.262 

-.031 

.281 

-.032 

.306 

-.046 

IDV/S_35m 

.326 

-.076 

.380 

-.159 

.326 

-.155 

.326 

-.147 

60m 

.287 

-.044 

.287 

-.086 

.287 

-.223 

lOQm 

.287 

-.137 

.287 

-.102 

IDV//?L35m 

.326 

-.088 

.326 

-.018 

333 

.326 

.022 

60m 

.323 

-.075 

.323 

-.215 

.323 

.074 

.323 

.141 

l(X)m 

.290 

-.032 

.290 

-.107 

El 

.015 

99 


Correlaticms  between  IDV  and  the  dynamic  stability  parameters  at  140  m,  200  m. 
and  300  m  were  weak  and  could  not  be  simply  related  to  the  diurnal  cycle  of  turbulence 
or  high-frequency  internal  waves  under  investigation  here. 

F.  RELATIONSHIP  BETWEEN  IDV  AND  MIXED-LAYER  DEPTH 

Mixed-layer  d^th  and  the  diurnal  range  of  mixed-layer  dq>th  provide  an 
indication  of  the  cumulative  effects  of  both  surface  forcing  and  the  static  and  dynamic 
stability  of  the  upper  ocean.  As  such,  it  represents  an  integral  of  the  processes  diought  to 
be  associated  with  upper  ocean  mixing  and  the  generation  of  high-frequency  temperature 
variability  in  the  mixed  layer  and  in  the  thermocline  below  the  mixed  layer.  Gxrelation 
coefficients  between  monthly  mean  mixed-layer  ctepths  ahd  the  diurnal  range  of  mixed- 
layer  depth  and  IDV  at  35  m,  60  m.  100  m,  140  m,  200  m,  and  300  m  for  the  four  time 
bins  are  shown  in  Table  5.5. 

The  correlations  between  mixed-layer  dq>th  and  IDV  at  35  m  were  significant  at 
the  95%  confidence  level  during  the  evening,  night,  and  morning  time  bins.  These 
correlations  were  particularly  strong,  -0.583  ami  -0.627,  during  the  evening  and  night  time 
bins,  respectively.  The  negative  sign  of  the  conelations  indicates  that  IDV  at  35  m 
increased  as  the  mixed  layer  deqxnwd,  in  agreement  with  the  microstructure  observations 
of  mixing  and  internal  waves  in  the  iq^)er  ocean  discussed  in  Chapter  2.  ^xed-layer 
dq[)th  and  IDV  at  35  m  were  not  well  conelated  during  the  day  time  bin,  when  tiie  mixed 
layer  was  shallow  due  to  daytime  lestratificatitm.  With  the  mixed  layer  being  shallow 
during  this  day  time  bin,  it  is  not  unreasonable  to  expect  a  decoiq>ling  with  IDV  at  35  m. 
The  strong  correlation  between  IDV  at  35  m  and  mixed-layer  dq>tii  during  the  evmiing 
time  bin  might  have  been  less  e}q)ected  since  the  mixed  layer  was  typically  shallow 
(Figures  4.9  and  4.11).  Close  examination  of  these  figures  shows  that,  exoq>t  for  during 
the  boreal  qning,  the  mixed  layer  was  de^wx  during  tiie  evening  time  bin  than  during  tiie 
day  time  bin,  indicating  that  the  active  nighttime  ctmvection  had  bqgun.  Significant 
negative  correlatimis  were  likewise  observed  between  mixed-layer  dq>tii  and  IDV  at  60 
m  during  the  morning,  day,  and  evening  time  bins  (Table  5.5).  The  corrdation  between 
mixed-layer  dq)tii  and  IDV  at  60  m  during  the  ni^t  time  bin  was  just  below  tiie  95% 
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confidence  level.  The  correlations  were  largest  during  the  morning  and  evening  time  bins, 
-0.497  and  -0.471,  respectively.  These  stronger  correlations  during  the  morning  and 
evening  time  bins  might  signify  that  during  those  hours  active  changes  were  occurring  in 
both  mixed-layer  depth  and  IDV  at  60  m.  IDV  at  100  m  and  140  m  was  not  significantly 
correlated  with  the  mixed-layer  depth  on  monthly  time  scales.  Positive  correlation 
between  IDV  at  200  m  and  mixed-layer  depth  was  observed  during  the  evening  and 
morning  time  bins.  Also,  positive  correlation  was  observed  between  mixed-layer  depth 
and  IDV  at  300  m  during  the  night  time  bin.  Although  it  is  difficult  to  determine  with 
certainty  whether  ex'  not  these  correlations  with  IDV  at  200  m  and  300  m  reflect  a  direct 
physical  relationship,  the  lack  of  significant  corelation  at  100  m  and  140  m  suggests  non¬ 
local,  large-scale  fexoing. 

The  diurnal  range  of  mixed-layer  depth  can  be  thought  of  as  a  measure  of  the 
intensity  of  mixing  on  diurnal  time  scales.  The  diurnal  range  of  mixed-layer  depth  was 
significantly  ctxrelated  with  IDV  at  35  m  during  the  night  and  morning  time  bins,  when 
IDV  was  at  its  diurnal  high  (Figures  S.1-S.3}.  These  correlations  indicate  that  when  the 
monthly  mean  diurnal  range  of  mixed-layer  depth  was  large,  indicating  high  diumally- 
modulated  mixing  intensity,  the  high-frequency  temperature  variability  at  35  m  was 
increased.  During  the  daytime  hours,  when  IDV  at  35  m  was  at  its  diurnal  low,  the 
correlation  with  diurnal  range  of  mixed-layer  dq>th  was  weak.  The  correlations  between 
IDV  at  60  m  and  100  m  and  diurnal  range  of  mixed-layer  deptii  were  not  significant. 
Significant  correlations  between  the  diurnal  range  of  mixed-layer  dq>th  and  IDV  at  140 
m  during  the  mtHtung  and  evening  time  bins  and  IDV  at  2(X)  m  during  the  night  and  day 
time  bins  cannot  be  simply  related  to  local  forcing.  As  with  the  significant  cmrelatims 
between  mixed-layer  depth  and  IDV  at  2(X)  m  and  3(X)  m,  these  relaticxishqK  between  the 
diurnal  range  of  mixed-layer  dq>th  and  IDV  at  140  m  and  200  m  are  probably  related  to 
large-scale  processes. 
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TaMe  5J.  CorrelatiiMi  coefficients,  r.  between  monthly  mean  isotherm  displacement 
variances  (IDV)  at  six  depdis  for  each  of  the  four  6-hour  time  bins  and  mixed-layer  dq)th 
(MLD)  and  diurnal  range  of  mixed-layer  depth  (AMLD)  during  the  period  November 
1983  -  October  1987.  Statistical  significance  at  the  95%  ctmfidence  level  was  calculated 
based  on  the  integral  time  scale  between  the  time  series  (Davis,  1976).  Statistically 
significant  correlations  are  indicated  by  Bold  type. 
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VI.  DIURNAL  CYCLE  OF  ISOTHERM  DISPLACEMENT 

VARIANCE 


In  the  preceding  two  chapters,  the  interannual  and  seasonal  variability  of  the 
surface  ftxoing,  background  dynamic  stability,  mixed-layer  dq>th.  and  high-frequency  IDV 
have  been  described,  and  statistical  relationships  among  these  fields  were  examined.  In 
this  chapter,  the  focus  shifts  to  the  diurnal  cycle  of  the  high-frequency  variations  of 
temperature,  as  measured  by  AIDV. 

A.  VERTICAL  EXTENT  OF  DIURNAL  CYCLE  OF  IDV 

The  time  soles  of  monthly  mean  IDV  at  five  depths  for  the  four  six-hour  time 
bins  presented  in  Figure  S.la-e  showed  pronounced  increases  of  IDV  during  the  night  and 
morning  time  bins,  collectively  referred  to  as  nighttime,  at  bodi  35  m  and  60  m  for  much 
of  the  48  month  record.  At  3S  m,  the  nigh^me  mtmdily-mean  IDV  increase  frequently 
exceeded  5  m^  ccmipaied  with  daytime  values.  During  a  few  months  this  nighttime 
increase  of  IDV  exceeded  10  m’.  At  60  m,  die  ni^tdme  monthly  mean  IDV  values 
frequendy  exceeded  daytime  IDV  values  by  about  5  m^.  At  100  m  and  below,  these 
diurnal  differences  in  monthly  mean  IDV  were  rarely  observed.  Similarly,  examination 
of  individual  mmithly  mean  vertical  profiles  of  IDV  (not  shown)  and  the  long-term  mean 
vertical  profile  of  IDV  (Figure  5.2)  showed  that  the  diurnal  cycle  of  mtmdily  mean  IDV 
only  rarely  penetrated  to  100  m.  Figure  5.2  also  indicated  hi^  variability  of  monthly 
mean  IDV  at  35  m  and  60  m  and  relatively  low  variability  of  mtmdily  mean  IDV  at  100 
m,  140  m,  and  200  m. 

Further  examinatim  of  the  vertical  extent  of  the  diurnal  cycle  of  IDV  was 
accomplished  by  computing  monthly  means  of  die  daily  differences  between  IDV  during 
different  six-hour  time  bins.  Since  diurnal  differences  of  IDV  i^ipeared  dqiendent  on 
the  magnitude  of  IDV,  i.e.  larger  differences  were  observed  when  IDV  was  high  (Rgure 
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Diurnal  Differences  of  IDV  (0400-1000  -  1600-2200) 
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Figure  6.1.  ^toithly  mean  normalized  diurnal  differences  of  isotherm 
variance  (AIDV^  computed  from  daily  values  of 0400*1000  IDV  minus  1600-2200  IDV 
at  a)  35  m,  b)  60  m,  c)  100  m,  d)  140  m,  and  e)  200  m.  Daily  diurnal  diffieienoes 
uwre  normalized  by  the  mean  of  the  moniing  and  evening  IDV  values  for  each  day. 
Dashed  lines  indicate  ±ooe  standard  error  ci  die  mean. 
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5.1).  the  diunud  differences  have  been  normalized  by  the  mean  of  the  daily  IDV  values 
for  the  two  time  bins.  Normalization  provides  a  relative  measure  of  the  diurnal  differences 
of  IDV,  allowing  mwe  direct  examination  ovo^  seasonal  and  interannual  time  scales.  The 
monthly  mean  ntxmalized  diurnal  differences  between  IDV  during  the  mtxning  time  bin 
and  IDV  during  the  evening  time  bin,  AIDV„  at  35  m,  60  m.  100  m,  140  m,  and  200 
m  are  shown  in  Figures  6.1a-e.  Hie  morning  time  bin  corresponds  with  the  times  when 
the  mixed  layer  was  deepest  (Figures  4.9  and  4.11).  Both  microstructure  measurements 
and  numerical  model  results  have  shown  the  diurnal  peak  in  turbulent  dissipation  in  the 
thermocline  (Mourn  et  al.,  1989;  Garwood  et  al.,  1989;  Schudlich  and  Price,  1992).  The 
evening  time  bin  corresponds  to  die  time  when  the  mixed  layer  was  still  shallow,  and 
microstructure  measurements  have  shown  weak  turbulent  dissipatimi  in  the  thermocline. 
Fot  comparison,  the  monthly  mean  diurnal  differences  between  IDV  during  the  night  time 
bin  and  IDV  during  the  day  time  bin.  AIDV,,,  are  shown  in  Figures  6.2a-e.  Examination 
of  the  IDV  time  series  suggests  that  differences  between  ATDV^  and  AIDV^  are 
generally  small. 

The  monthly  means  of  AIDV^  at  35  m  show  that  the  morning  IDV  is  almost 
always  greater  than  the  evening  IDV,  and  diat  for  81%  of  the  record  this  difference  is 
significant  within  a  standard  eirw  of  the  mean  (Figure  6.2a).  The  AIDV„  at  60  m  shows 
that  for  65%  of  the  record  length  the  morning  IDV  is  greater  than  the  evening  IDV  by 
at  least  one  standard  error  of  the  mean  (Hgure  6.2b).  By  die  100  m  dqith,  AIDV.^  is 
greato*  than  one  standard  error  of  the  mean  only  29%  of  the  record,  a  considerable 
decrease  from  diat  at  60  m  (Figure  6.2c).  At  140  m  and  200  m,  diese  percentages  were 
reduced  to  17%  and  12%,  respectively  (Figures  6.2d,e). 

Examination  of  AID V^  shows  similar,  tfaou^  weaker,  diurnal  differences  of  IDV 
(Figure  6.2a).  AIDVg^  at  35  m  shows  that  the  ni^t  IDV  values  were  greater  than  the  day 
IDV  valu^  by  at  least  one  standard  error  for  72%  of  the  record  (Figure  6.2a).  At  60  m, 
this  percentage  was  reduced  to  48%.  At  100  m,  140  m,  and  700  m,  Uiese  percentages 
were  29%,  8  %,  and  12  %,  respectively  ^ven  though  each  of  diese  percentages  were 
smaller  than  those  cominited  for  AIDV,^  the  same  general  pattern  of  decreased  diurnal 
differences  of  IDV  with  dqith  was  observed. 
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Diurnal  Differences  of  IDV  (2200-0400  -  1000-1600) 
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Figure  62,  Monthly  mean  diurnal  differences  of  isothenn  diq>laoenient  variance 
(AIDVn^  computed  from  daily  valiws  of  2200-0400  IDV  minus  1000-1600  IDV 
normalized  by  die  mean  of  the  two  IDV  values  at  a)  35  m,  b)  60  m,  c)  100  m.  d)  140 
m,  and  e)  200  m.  Dariied  lines  indicate  ±oat  standard  error  of  die  mean. 
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Figure  63.  Long-tain  mean  vertical  protiles  of  a)  Isotherm  diylacement  variance  for 
the  0400-1000  (solid  line)  and  1600-2200  (dashed  line)  time  bins,  b)  zonal  velocity,  and 
c)  temperature  confuted  from  all  available  mmithly  means  during  the  paiod  November 
1983  -  October  IS^.  Standard  error  of  the  mean  is  indicated  by  arrows  in  a),  and 
stamlard  deviations  are  indicated  by  dashed  lines  in  b)  and  c). 
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To  further  evaluate  the  vertical  structure  of  the  diurnal  cycle  of  high-frequency 
IDV  and  possible  long-term  consequences  to  the  equatorial  mcMnentum  balance,  statistics 
of  the  vertical  profile  of  IDV  were  computed  from  the  48  monthly  mean  values  for  each 
of  the  four  six-hour  time  bins  (Figure  5.2).  These  long-term  mean  IDV  profiles  were 
consistent  with  the  above  results  that  showed  a  pronounced  diunud  difference  of  IDV  at 
35  m  and  60  m,  particularly  between  the  morning  and  evening  time  bins.  The  long-term 
mean  IDV  profiles  for  these  two  time  bins  and  the  long-term  mean  vertical  profiles 
of  temperature  and  zonal  velocity  are  shown  in  Figures  6.3a-c.  Figure  6.3a  shows 
pronounced  diurnal  differences  in  IDV  at  35  m  and  60  m,  with  the  mean  morning  IDV 
greater  than  the  mean  evening  IDV  by  about  4  m^  and  2  m^.  respectively.  This  figure  also 
shows  a  sharp  decay  with  depth  of  the  long-tenn  mean  diurnal  cycle  of  high-frequency 
IDV.  Although  the  time  series  of  monthly  mean  diurnal  differences  showed  occasional 
penetration  of  AIDVg^  to  100  m  and  deqxr  (Figure  6.1),  the  long-term  mean  clearly 
shows  that  on  average  the  diurnal  cycle  was  not  evident  at  100  m.  140  m,  or  200  m. 
Figure  6.3  also  shows  the  depths  of  the  observed  AIDV^e  relative  to  the  mean  EUC  and 
thermocline.  In  the  long-term  mean,  the  AIDV^  did  not  penetrate  through  the  EUC  core 
or  high-V  region  of  the  thermocline.  As  was  shown  in  the  time  series  (Figure  6.2),  the 
long-term  mean  diurnal  difference  of  IDV  between  night  and  day  time  bins  is  ai^arent, 
but  weaker  than  the  morning/evening  difierence. 

Another  method  of  examining  the  vertical  structure  of  the  high-fiequency  AIDV 
is  to  compute  vertical  correlations  between  the  AIDV  at  different  dq>tbs.  The  advantage 
of  this  method  is  that  it  provides  an  indication  of  the  relationship  between  the  diurnal 
cycle  of  the  high-frequency  variability  at  diffoent  depths.  Strong  cmrelations  over  a  large 
vertical  extent  would  be  consistent  with  hypotheses  suggesting  vertical  internal  wave 
propagation.  Correlation  coefficients  between  the  monthly  mean  AIDV^  at  35  m  and 
AIDV^  at  60  m,  100  m,  140  m,  200  m.  and  300  m  are  shown  in  Figure  6.4.  Based  <m 
the  result  that  the  most  pronounced  diurnal  difference  in  IDV  occurred  betwem  the 
mmning  and  evening  time  bins,  this  is  the  diurnal  difference  used  in  the  computation  of 
the  correlation  coefficients  shown  in  Figure  6.4.  Mean  profiles  of  zcmal  velocity  and 
temperature  illustrate  the  posititms  of  the  EUC  and  thermocline  relative  to  the  vortical 
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structure  of  AIDV. 

The  correlation  between  the  monthly  mean  AIDV,,,  at  35  m  and  AIDV^  at  60  m 
was  strong,  +0.61,  indicating  that  when  the  AIDV„  at  35  m  increased  there  was  likely 
to  be  a  corresponding  increase  of  AIDV„^^  at  60  m,  and  vice  versa.  Hie  correlations 
between  AIDV^  at  35  m  and  AIDV„  at  100  m,  140  m,  200  m,  and  300  m  were  not 
significant  at  the  95%  confidence  level. 

On  shcMter  times  scales,  the  15-minute  temperature  measurements  were  used  to 
examine  vertical  propagation  by  computing  correlations  as  a  function  of  time  lag  (every 
15  minutes)  and  depth.  In  addition  to  sharp  decay  with  depth,  in  agreement  with  monthly 
mean  results,  those  correlations  indicated  highest  vertical  correlations  at  zero  lag, 
suggesting  little  downward  phase  propagation. 

These  combined  results  indicate  that  vertical  propagation  of  high-frequency 
internal  waves  from  the  surface  layer,  if  occurring,  did  not  appear  to  propagate  through 
the  EUC  core  using  monthly-averaged  statistics.  These  results  do  not  rule  out  the 
possibility  that  a  diurnal  cycle  of  IDV  may  be  occasionally  observed  at  depths  greater 
than  60  m.  For  example,  AIDV^  at  100  m  was  greata*  than  one  standard  error  of  die 
mean  for  29%  of  the  months  examined  (Figure  6.1).  However,  on  average  the  results 
indicate  that  the  diurnal  cycle  of  IDV  decayed  sharply  below  60  m. 
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Correlation  of  Diurnal  Cycle  of  IDV,  U  &  T 
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Figure  6.4.  Vertical  profiles  of  a)  correlatitm  coefficients  between  monthly  mean  diurnal 
differences  of  isotherm  displacement  variance  (AIDV^  at  35  m  and  those  at  60  m,  100 
m,  140  m.  200  m,  and  300  m,  b)  Itmg-term  mean  zonal  velocity,  and  c)  long-term  mean 
temperature.  Dnhed  lines  indicate  a)  95%  confidrace  intervals  based  on  intq^ral  time 
scales  between  time  series  (Davis,  1976),  and  b)  and  c)  standard  deviations. 
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B.  INTERANNUAL  VARIABILITY  OF  THE  DIURNAL  CYCLE  OF  IDV 


In  discussing  the  interannual  variability  of  IDV  for  the  four  six-hour  time  bins 
shown  in  Figures  S.la-e,  mention  was  made  of  the  interannual  variability  of  the  diurnal 
cycle  of  IDV  at  35  m  and  60  m.  Further  examination  of  Figures  5.1a,b  shows  that  during 
the  1984  ENSO  cold  episode,  a  diurnal  cycle  of  IDV  was  not  observed  at  either  35  m  or 
60  m.  This  disiq>pearance  of  the  diurnal  cycle  during  the  cold  episode  was  also  q^arent 
in  the  time  series  of  monthly  means  of  AIDV„  (Figures  6.1a.b).  Beginning  in  mid-1984, 
AIDV  at  35  m  reappeared  and  remained  present  during  most  of  the  next  two  years  (Figure 
6.1a).  Comparison  with  the  subsurface  temperature  record  (Figure  4.7)  shows  that  the 
reappearance  coincided  with  the  beginning  of  a  two-year  deepening  of  the  thermocline. 
Although  the  daytime  IDV  values  increased  during  this  ENSO  build-up  phase,  it  appears 
that  AIDVn,  was  dominated  by  larger  increases  of  IDV  during  the  nighttime  hours 
(Figure  5.1a).  During  the  1986/1987  ENSO  warm  q>isode,  this  did  not  q>pear  to  be  the 
case.  Both  daytime  and  nighttime  IDV  values  increased  above  the  annual  boreal  winter 
maximum  (Figure  5.1a)  and  AIDV„  decreased  (Figure  6.1a).  Although  the  long-term 
mean  vertical  profiles  of  IDV  did  not  show  a  diurnal  cycle  at  100  m  (Figures  5.2  and 
6.3),  diurnal  differences  of  IDV  at  100  m  wo'e  not  uncommon  during  the  ENSO  build-up 
phase  between  November  1984  and  March  1986  (Figure  6.1c). 

C.  SEASONAL  CYCLE  OF  DIURNAL  CYCLE  OF  IDV 

As  discussed  in  section  5.3,  there  was  a  {x-onounced  seasonal  cycle  of  IDV  at  35 
m  and  60  m  and  a  weaker  seasonal  cycle  of  IDV  at  1(X)  m.  Examination  of  Figure  5.2 
shows  a  similar  seasonal  cycle  of  the  diunud  IDV  signal,  as  indicated  by  the  relative 
divergence  of  the  nighttime  and  daytime  IDV  curves  at  35  m  and  60  m.  As  with  IDV, 
AIDV  was  oteerved  to  be  strongest  during  the  boreal  winter  months  and  weakest  during 
the  boreal  spring.  Seasonal  cycles  of  AIDV  from  the  monthly  mean  values  of  AIDV„e 
(Figure  6.1)  are  shown  in  Figure  6.5.  Hgures  6.5a.b  indicate  a  pronounced  annual  cycle 
of  AIDV  at  35  m  and  60  m,  with  maxima  occurring  during  the  boreal  winter  and  minima 
occurring  during  the  boreal  spring.  Indications  of  an  apparent  semi-annual  cycle  of  AIDV 
at  35  m  were  also  obsoved.  An  (q>parent  weak  annual  cycle  of  AIDV  at  100  m  was  observed 
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Seasonal  Cycles  of  Diurnal  Differences  of  IDV 


Figure  6^.  Seasonal  cycles  of  the  diurnal  difference  ci  IDV  at  a)  35  m,  b)  60  m.  c) 
100  m,  d)  140  m,  and  e)  200  m.  Computed  from  normalized  mondily  mean  differences 
of  IDV  between  0400>1000  and  1600*2000.  firor  bars  represent  aven^  of  all  the 
individual  standard  errors  for  eadi  of  die  twelve  ob>nrf«r  months. 


with  relative  maxima  occurring  during  die  boreal  summer  mcmths  and  near  zero  diurnal 
differences  occurring  during  the  boreal  winter  months.  This  weak  aimual  cycle  at  100  m. 
opposite  in  phase  to  the  annual  cycles  of  AIDV  at  3S  m  and  60  m.  was  not  significant 
based  on  ±one  standard  error  of  the  mean.  The  seasonal  variability  of  AIDV  was  not 
significant  at  both  140  m  and  200  m  since  diurnal  variations  of  IDV  were  usually  not 
evident. 

D.  SURFACE  FORCING  AND  THE  DIURNAL  CYCLE  OF  IDV 

In  this  section,  statistical  relationships  between  the  diurnal  cycle  of  IDV  and 
several  indices  of  surface  forcing  are  (vesented.  Based  on  the  earlier  finding  that  the 
strongest  diurnal  difference  of  IDV  occurred  between  the  morning  and  evening  time  bins, 
monthly  means  of  ATDV,^  at  35  m,  60  m,  and  100  m  were  used  as  the  reference  of  the 
diurnal  cycle  of  IDV  to  examine  these  relatimishqis  (Figure  6.1).  With  the  jxior  finding 
that  there  was  no  qiparent  diurnal  cycle  of  IDV  at  140  m,  200  m,  and  300  m,  the 
statistical  relationsh4>s  between  AIDV^  at  these  dqiths  and  the  surface  forcing  indices 
are  not  of  interest.  The  monthly  mean  time  series  of  each  of  the  surface  forcing  conditions 
described  in  section  4.2.1  were  correlated  widi  AIDV„e- 
1.  Surface  Winds 

In  section  5.4,  surface  wind  speed  and  zonal  wind  speed  were  shown  to  be 
strongly  ctxrelated  with  IDV  at  35  m  and  60  m  mi  mondily  time  scales.  In  this  sectimi, 
cmrelations  between  monthly  mean  wind  speed  and  directitm,  and  the  zmial  and 
meridional  wind  compmients,  and  mtmdily  mean  AIDV^  at  35  m,  60  m,  and  1(X)  m  are 
examined.  Table  6. 1  shows  these  correlations  and  dieir  significance  at  the  95%  confidence 
level  based  on  the  integral  time  scale  between  time  smies  (Davis,  1976).  Mondily  mean 
AIDV^  at  35  m  had  positive  correlation  with  wind  speed,  0.345,  and  significant 
negative  correlaticm,  -0.361,  with  z<mal  wind.  AIDV^  at  60  m  had  a  strtmg  positive 
correlation  with  wind  speed,  •^0.557,  and  strcmg  n^adve  correlatimi  with  zonal  winds,  - 
0.580.  These  strong  cmrelations  and  an  examinatimi  of  the  monthly  mean  time  series  in 
Figure  6.6  show  that  AIDV^  increased  with  increased  wind  speed  and  easterly  zonal 
winds.  It  is  important  to  note  diat  since  diurnal  variatimis  of  the  local  winds  were  small. 
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Table  6.1.  Correlatioii  coefficients  between  normalized  mtmthly  mean  diurnal  differences 
of  isotherm  di^lacement  variance  (AIDV)  at  3Sm.  60m,  and  100m  and  surface  wind 
speed,  zonal  and  meridional  winds,  and  wind  direction  during  the  period  November  1984 
-  October  1987.  AIDV  computed  as  0400-1000  IDV  minus  1600-2200  IDV.  Statistical 
significance  at  the  95%  confidence  level  was  calculated  based  on  the  integral  time  scale 
between  the  time  series  (Davis,  1976).  Statistically  significant  correlations  are  indicated 
by  Bold  type. 


Variable  1 

Variable  2 

95%  C.L 

r 

Wind  Speed 

AIDV_035m 

±0.353 

0.345 

±0.392 

0.557 

alDV.lOOm 

±0.324 

0.216 

Zonal  Wind 

AlDV.OaSm 

±0.354 

•0361 

±0.397 

-0380 

alDV.lOOm 

±0.321 

-0.230 

Merid.  Wind 

AlDV_035m 

±0.351 

0.023 

AlDV_060m 

±0.428 

0.133 

alDV.lOOm 

±0.238 

-0.022 

Wind  Dir. 

AlDV_035m 

±0.353 

0.102 

AlDV_060m 

±0.396 

0.110 

alDV.lOOm 

±0.334 

-0.061 

about  0.2  ms'‘,  it  is  unlikely  that  these  correlations  between  the  wind  q>eed  and  zonal 
wind  and  AIDV„  at  35  m  and  60  m  were  related  to  diurnal  fluctuations  of  the  local 
winds.  The  stronger  correlatitms  at  60  m  could  have  been  due  to  the  shorter  record  length 
and  large  data  gaps  at  60  m  during  the  boreal  summers.  The  mcmthly  mean  values  of 
AIDV„  at  100  m  were  not  well  correlated  with  wind  speed  or  zonal  wind.  The  monthly 
mean  meriditmal  wind  and  wind  direction  were  not  well  correlated  witii  AIDV^  at  any 
of  the  TR  depths. 
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2.  Surface  Heat  and  Buoyancy  Fluxes 

In  Cluster  V.  statistical  relationships  between  monthly  mean  IDV  and 
monthly  mean  estimates  of  the  surface  heat  and  buoyancy  fluxes  were  examined. 
Generally,  correlations  were  weak.  The  <mly  significant  correlatimis  were  between  IDV 
at  35  m  and  Q,,  and  during  the  evening  time  bin.  In  this  section,  relatitmships 
between  the  surface  heat  and  buoyancy  fluxes  and  the  diurnal  cycle  of  IDV  in  the  upper 
100  m  are  examined.  Monthly  means  of  tl»  confuted  estimates  of  the  terms  of  the 
surface  heat  flux  and  AIDV^  at  35  m  and  60  m  are  shown  in  Figure  6.7.  Correlatim 
coefficients  between  die  surface  heat  and  buoyancy  fluxes  and  AIDVg^  at  35  m,  60  m,  and 
100  m  are  given  in  Table  6.2. 

Earlier  results  showed  that  IDV  in  the  upper  100  m  was  not  well  correlated 
with  Q,.  The  correlations  between  mondily  mean  AIDV^  at  35  m  and  60  m  and  were 
0.434  and  0.453,  respectively.  The  conelation  widi  AIDV^  at  35  m  is  statistically 
significant  at  the  95%  ctxifidence  level  based  on  the  integral  time  scale  of  the  two  time 
series  (Davis,  1976).  Although  the  ccrrelatiaii  widi  AIDV^  at  60  m  was  sli^tly  greater 
than  that  at  35  m,  it  was  not  statistically  significant  due  to  greater  confidence  intervals 
associated  with  the  shorter  record  length.  The  positive  sign  of  these  correlations  indicates 
that  AIDV^  increased  during  periods  of  increased  shortwave  radiadoo  into  the  ocean. 
When  monthly  mean  Q,  was  anomalmisly  low,  such  as  during  the  1984  ENSO  cold 
qnsode  and  during  the  wanning  in  1987.  the  diurnal  differences  of  IDV  at  35  m  and  60 
m  were  low  (Figure  6.2).  Althou^  iimeased  shortwave  radiation  would  increase  die 
daytime  stratification  of  the  upper  ocean,  and  hence  reduce  the  daytime  generation  (rf 
ctmvective  mixing  and  possibly  inleraal  wave  generation,  the  fact  that  Q,  is  n^gible 
during  the  nighttime  hours  means  that  an  increase  in  daytime  shortwave  radiation 
corresponds  to  an  increase  in  the  diurnal  difference  in  shortwave  radiation.  It  is  not 
unreasonable  to  eiqiect  an  increase  in  die  diurnal  response  of  IDV  when  there  was  an 
increase  in  the  diurnal  solar  forcing,  and  vice  versa.  Therefore,  these  results  could  be 
intenmted  to  indicate  diat  increased  results  in  an  increase  in  the  diurnal  difference 
of  IDV  because  of  reduced  daytime  ^neration  of  IDV  rather  thm  due  to  increased 


116 


Diurnal  Cycle  of  IDV  and  Surface  Heat  Fluxes 


1983  1984  1985  1986  1967 


1983  1984  1985  1986  1987 


1983  1984  1985  1986  1967 

Figure  6.7.  Time  series  of  monthly  mean  a)  diurnal  differences  of  isotherm 
displacement  variance  (AIDV^  at  35  m  and  60  m,  b)  net  surface  net  flux,  c) 
shortwave  radiation.  d)  latent  heat  flux,  and  e)  longwave  radiation,  (2».  and 
sensible  heat  flux.  during  Novembv  1983-October  1987.  AIDV.^  computed  as  0400- 
1000  IDV  minus  1600-2200  IDV  normalized  by  die  mean  IDV. 
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Table  6JI.  Correlation  coefficients  between  normalized  mmithly  mean  diurnal  differences 
of  isotherm  displacement  variance  (AIDV)  at  35  m,  60  m.  and  100  m  and  surface  heat 
and  buoyancy  fluxes,  including  net  surface  heat  flux.  shortwave  radiation,  latent 
heat  flux,  Q^,  sensible  Iwat  flux,  (2*>  longwave  radiation.  net  surface  buoyancy  flux. 
Bg,(.  and  ev^Ktration,  E,  during  the  poiod  November  1983-October  1^7.  AIDV 
computed  as  0400-1000  IDV  minus  1600-2200  IDV.  Statistical  significance  at  die  95% 
confidence  level  was  calculated  based  on  tl^  integral  time  scale  between  the  time  series 
(Davis,  1976).  Statistically  significant  cmrelations  indicated  by  Bold  type. 
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nighttime  generatitm  of  IDV.  Although  ctxreUtions  were  weak,  crmiparistMi  between 
daytime  values  of  IDV  at  35  m  and  60  m  (Figures  S.la.b)  and  Q,  supports  this 
interpretation. 

E)aytime  IDV  values  were  lowest  during  early  1984  and  during  most  of 
1987,  both  periods  of  anomalously  low  Q^. 

The  only  other  term  of  the  surface  heat  balance  having  significant 
ctMTelation  with  AIDV^in  the  upper  100  m  was  Q,,.  The  correlations  between  and  the 
diurnal  differences  of  IDV  at  35  m  and  60  m  were  -0.359  and  -0.415,  req)ectively.  The 
negative  sign  of  the  correlations  indicates  that  AIDV^  increased  during  periods  when  the 
sensible  heat  flux  was  from  the  ocean  to  the  atmosphere.  Although  the  sign  agrees  with 
intuition,  the  small  relative  magnitude  of  Q„  would  not  f4>pear  sufficient  to  generate 
significant  convective  mixing.  Since  (2*  is  a  fimction  of  wind  speed,  this  significant 
correlation  is  most  likely  due  to  the  strong  ccarelation  between  wind  q)eed  and  AIDV^. 

The  remaining  terms  of  the  surface  heat  flux  and  the  moisture  flux  terms 
of  the  surface  buoyancy  flux  were  not  well  correlated  with  AIDV^. 

E.  DYNAMIC  STABILITY  AND  THE  DIURNAL  CYCLE  OF  IDV 

In  this  section,  relationships  between  the  dynamic  stability  during  each  of  the  four 
six-hour  time  bins  and  AIDV^  at  35  m,  60  m,  and400  m  are  examined.  In  Figures  4.8a- 
c,  time-dq>th  crmtour  plots  of  monthly  mean  N,  5,  and  Ri  were  i»esented.  In  examining 
the  relationships  between  the  diurnal  cycle  (tf  IDV  and  the  dynamic  stability,  diurnal 
variations  of  diese  dynamic  stability  parametm  are  discussed. 

1.  Buoyancy  FrequoiQr 

Earlier  results  showed  diat  N  varied  seasonally  and  interannually,  and  that 
these  variatirms  of  V  at  several  dq>ths  in  die  iqiper  100  m  were  negatively  correlated  with 
IDV  at  35  m  and  60  m,  but  not  with  IDV  at  100  m.  The  correlatitms  between  the  mtmdily 
mean  AIDVg^  at  35  m,  60  m,  and  100  m  and  monthly  means  of  local  values  of  N  for  each 
of  the  four  six-hour  time  bins  are  {xesented  in  Figure  6.8a.  Significant  negative 
correlation  was  observed  at  35  m  and  60  m  for  each  of  the  four  time  bins.  The  negative 
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Correlations:  Diurnal  Difference  of  IDV  &  N,  S,  &  Ri 


a)  Ba«juM7  HraqiMiMgr  b)  Valoelt^  Shwur  c) 


Figure  61^.  Correlation  coefficients  confuted  between  diurnal  differences  of  isotbenn 
displacement  variance  (AIDV)  and  a)  buoyancy  frequency,  N,  b)  vdiocity  shear,  S,  and 
c)  Richardson  number,  logioC^O  dutiog  November  1983-October  1987.  AIDV  CQnq>uted 
as  0400-1000  IDV  minus  1600-2200  IDV  normalized  by  mean  IDV. 
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sign  of  the  correlations  indicates  that  alDV^  at  3S  m  and  60  m  increased  during  periods 
of  low  N,  and  decreased  during  periods  of  high  N.  This  suggests  that  strong  stratification 
may  have  aided  in  suppressing  the  diurnal  cycle  of  hirbulence  and  high-frequency  internal 
waves.  Although  daily  values  of  AEDVn,  were  generally  eliminated  when  the  TRs  were 
in  the  mixed  layer  before  computing  monthly  statistics,  the  (xxrespondence  between  high 
AIDVg^  and  low  N  suggests  that  the  instruments  were  in  the  entrainment  zone  and  that 
active  mixing  may  have  been  occurring.  If  so.  it  would  suggest  that  AIDV^  is  largely  a 
reflection  of  turbulence  rather  than  internal  waves. 

Diurnal  differences  were  computed  from  the  fields  of  monthly  mean  N 
values  between  each  of  the  four  time  bins.  The  resulting  difference  fields  showed  that  the 
maximum  diurnal  differences  occurred  between  the  mcMuing  and  evening  time  bins,  as 
was  observed  for  diurnal  differences  of  IDV.  Figure  6.9  is  a  time-depth  contour  plot  of 
the  monthly  mean  evening  minus  morning  differences  of  N,  Large  were 
observed  in  the  upper  40  m  or  so.  reflecting  the  strong  diurnal  mixing  cycle  associated 
with  daytime  stratification  by  solar  heating  and  nighttime  convective  mixing.  In  the  dq}th 
range  of  about  50-100  m.  the  sign  of  appeared  to  alternate  with  an  annual  period. 
During  the  boreal  fall  and  winter  months,  values  of  in  this  dq}th  range  were  smaU 
and  positive,  indicating  that  stratification  increased  slightly  during  the  evening  time  bin. 
This  penetration  of  the  diurnal  increase  of  N  into  the  thermocline  coincided  with  the 
pronounced  increase  of  the  diurnal  cycle  of  IDV  at  35  m  and  60  m  during  the  winter 
mrniths  (Figure  5.3a,b).  Thus,  this  diurnal  cycle  of  V  in  die  thermocline  t^pears  to  reflect 
the  deep  diurnal  cycle  of  turbulence  and  internal  waves  that  has  been  reported  (Peters  et 
al.,  1988;  Mourn  et  al.,  1989;  McPhaden  and  Peters,  1992).  During  the  boreal  spring  and 
summer  months,  values  of  in  this  50-100  m  depth  range  were  strongly  negative, 
indicating  weaker  stratification  during  the  evening  time  bins,  i.e.  t^posite  the  evening 
increase  in  stratification  observed  in  the  surface  layer.  The  opposite  sign  of  AN^  across 
the  base  of  the  mixed  layer  during  these  months  indicates  conventitmal  mixed-layer 
dynamics  with  a  diurnal  buoyancy  flux  at  the  base  of  the  mixed  layer  balanced  by 
diumally-steady  upwelling  by  Ekman  pumping.  These  seasonal  variations  of  AN^  were 
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Figure  6.9.  Time-depth  contour  plot  of  diurnal  differences  of  monthly  mean  a) 
buoyancy  frequency,  during  November  1983-October  1987.  A/V„  computed  as 

1600-2200  N  minus  0400-1000  N.  Units  are  10^*  x  s'*. 
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observed  to  vary  interannually  as  well.  During  both  the  1984  ENSO  cold  episode  and 
the  1986/1987  ENSO  warm  episode,  the  daytime  stratification  of  the  surface  layer  was 
noticeably  weaker. 

Although  pronounced  diurnal  differences  of  ID V  and  N  were  observed  and 
AIDVok  at  35  m  and  60  m  was  well  correlated  with  N,  cmrelations  between  and 
AIDVqn  St  35  m,  60  m,  and  100  m  were  not  significant  at  the  s5%  conEdence  level.  The 
fact  that  the  diurnal  cycle  of  IDV  was  correlated  with  N,  but  not  AA/„,  indicates  that 
stratification  in  the  surface  layer  may  have  provided  a  necessary,  but  not  sufficient, 
condition  for  the  diurnal  cycle  of  the  high-frequency  variability. 

2.  Velocity  Shear 

In  this  section,  the  aurelations  between  the  diurnal  cycle  of  IDV  and  shear, 
5,  and  its  diurnal  variations  are  examined.  Figure  6.8b  shows  local  correlations  between 
monthly  mean  AIDV^  and  S  at  35  m,  60  m,  and  100  m  for  each  of  the  four  six-hour  time 
bins.  At  each  of  these  three  depths,  cmrelations  between  S  and  AIDV^g,  were  weak. 
Statistically  significant  corelation  at  the  95%  confidence  level  was  found  only  at  100  m 
and  only  during  the  day  time  bin. 

Diurnal  differences  of  monthly  mean  5  between  each  of  the  four  six-hour 
time  bins  showed  that  the  maximum  diurnal  range  occurred  between  the  morning  and 
evening  time  bins,  as  was  the  case  fcx  IDV  and  N.  Figure  6.10  is  a  time-d^th 
contour  plot  of  AS^.  Large  AS„  was  observed  in  the  surface  laya*,  which  by  this  measure 
penetrated  to  between  about  25  m  and  50  m.  Large  positive  AS^  in  the  surface  layer 
indicates  that  shear  increased  during  the  evening  time  bin,  a  result  consistent  with 
observations  and  modeling  studies  showing  a  late  afternoon  diurnal  jet  caused  by  the 
increased  stratification  (Montgomery  and  Strmip,  1962;  Schudlich,  1991).  The  vertical 
extent  and  magnitude  of  AS„  varied  seasonally  and  interarmually.  On  a  few  occasions, 
surface-related  AS^  q>peared  to  penetrate  as  deep  as  100  m,  most  notably  during  the 
1986/1987  ENSO  warm  q)isode  with  the  passage  of  a  pronounced  Kelvin  wave  in 
November  1986.  The  seasonal  cycle  of  AS^  in  the  surface  layer  was  characterized  by 
minima  during  the  Ixxeal  spring  and  maxima  during  the  bweal  summer  and  fall. 
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Diurnal  Difference  of  Velocity  Shear 


Figure  6.10.  'nme-dq>th  contour  plot  of  diurnal  differences  of  monthly  mean  velocity 
shear,  during  November  1983-October  1987.  computed  as  1600-2200  S  mitii.a 
0400-1000  S.  Units  are  10^*  x  s‘‘. 


124 


As  with  the  diurnal  cycle  of  N,  correlations  between  the  diurnal  cycle  of  S  and  the 
diurnal  cycle  of  IDV  at  35  m.  60  m,  and  100  m  were  not  statistically  significant  at  the 
95%  confidence  level. 

3.  Richardson  Number 

Local  correlations  between  IDV  and  Ri  in  the  upper  100  m  were  shown  to 
be  weak.  Significant  non-local  correlations  were  found  between  IDV  at  60  m  and  Ri  at 
80-120  m  during  the  morning  time  bin  and  between  IDV  at  100  m  and  Ri  at  10  m  during 
the  morning  and  evening  time  bins.  In  this  section,  the  relationships  between  AIDVg,^  and 
Ri  and  the  diurnal  variations  of  Ri  are  examined.  Figtire  6.8c  shows  the  local  correlations 
between  monthly  mean  AIDV,^  and  Ri  at  35  m.  60  m,  and  100  m  for  each  of  the  four 
six-hour  time  bins.  Correlations  were  weak  at  all  three  depths  during  all  four  time  bins 
indicating  that  AIDV^  and  Ri  were  not  strongly  related  on  these  time  scales. 

Difierences  between  monthly  mean  Ri  for  each  of  the  four  six-hour  bins 
showed  that  the  maximum  diurnal  range  of  Ri,  ^Ri,  occurred  between  the  morning  and 
evening  time  bins.  A  time-d^tfa  contour  plot  of  Figure  6.11,  shows  large 
occurred  over  seasonal  and  interannual  dme  scales.  llne]q)ected  positive  ARi^ 
indicates  that  surface  layer  Ri  values  w«e  usually  larger  during  the  morning  time  bin  than 
during  the  evening  time  bin.  Generally,  this  positive  penetrated  to  about  25-50  m. 
There  were  occasions,  however,  when  the  votical  extent  of  these  positive  ^i^  extended 
to  about  100  m,  such  as  during  the  boreal  winter  of  the  1986^1987  ENSO  warm  q>isode. 
I>uring  the  1984  ENSO  cold  q)isode,  small  were  observed  in  the  surface  layer 
while  large  positive  were  observed  between  about  60  m  and  130  m.  The  ARi^  in 
this  latter  depth  range  alternated  in  sign,  with  occasitMud  periods  oil  lar^  negative 
occurring  during  the  1986/1987  ^SO  warm  q>isode. 

Correlations  between  AIDV^,  and  were  not  significant  at  35  m,  60 
m,  and  100  m.  The  values  of  AIDV,.  at  35  m  were  significantly  crnrelated  with  the  night 
minus  day  Ri  diffo-ence,  ^i,^,  at  35  m.  This  correlation  was  -0.36,  indicating  that  the 
moming/evening  AIDV  at  35  m  increased  when  the  night/day  diurnal  difference  of  Ri 
decreased.  Physically,  it  is  unclear  how  a  decrease  in  the  diurnal  difference  in  Ri  could 
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Figure  6.11.  Time-depth  contour  plot  of  diurnal  differences  of  mmthly  mean  Richardscm 
number,  ARi^  during  November  1983-October  1987.  computed  as  0400-1000  Ri 
minus  1600-2200 /{i. 


be  related  to  an  increase  in  the  diurnal  cycle  of  turbulence  or  internal  waves. 

4.  Longterm  Mean  Dynamic  Stability 

Although  the  correlations  between  mondily  mean  AIDV^  and  the  mmthly 
mean  diurnal  differences  of  dynamic  stability.  AA/^,  and  were  weak,  the 
longterm  mean  diurnal  differences  suggest  possible  relationships.  Figure  6.12  shows  the 
four-year  mean  vertical  profiles  of  ATDV^,  AS„.  and  AJ?(„  and  the  standard  errtxs 
of  the  means.  The  longterm  mean  vertical  structure  of  the  diurnal  cycle  of  IDV  was 
discussed  in  section  6.1.  Unlike  Figures  5.2  and  6.3.  which  showed  the  Icmgterm  mean 
IDVs  for  each  of  the  four  time  bins.  Figure  6.12a  shows  the  Itmgterm  mean  AIDV^. 
Recall  that  AID was  computed  as  the  monthly  mean  of  the  daily  mcxoing  IDV  minus 
the  evening  IDV  normalized  by  the  daily  mean  of  the  two  IDV  values.  It  is  noteworthy 
that  the  standard  err<xs  of  the  mean  AID Vg.  are  much  smaller  than  the  standard  errms  of 
the  mean  IDV.  Although  by  this  measure  of  significance  AIDV^  at  100  m  was 
significant,  the  sharp  decay  with  dqith  of  AIDV.^  below  60  m  agrees  with  the  sharp 
decay  with  dqith  of  IDV  described  previously.  This  difference  in  significance  is  a  result 
of  the  normalization  procedure,  which  acts  to  reduce  die  mondi-to-mondi  variability  by 
providing  a  relative  measure  of  the  diurnal  difference  of  IDV. 

The  Icmgterm  mean  A/V^  was  signiffcantly  positive  in  the  upper  40  m  aixl  slightly 
negative  between  about  45  m  and  100  m  to  me  standard  error  of  the  mean  (Figure  6. 12b). 
In  discussing  Figure  6.9,  a  pronounced  seasonal  cycle  of  this  change  of  sign  of  A/V  across 
the  base  of  the  mixed-layer  was  described.  During  the  boreal  winter,  the  diurnal  cycle  of 
N  penetrated  as  deep  as  1(X)  m,  reflecting  the  deep  diurnal  cycle  ci  turbulent  mixing. 
During  the  Ixxeal  summer,  AA/^  was  negative  in  the  thermocline  region,  indicating 
conventional  mixed-layer  physics  widi  a  diurnal  buoyancy  flux  at  the  base  of  the  mixed 
layer  balanced  by  Ekman  pumping.  In  the  Imgterm  mean  (Figure  6.12b),  a  weak  diurnal 
buoyancy  flux  from  the  surface  layer  to  the  thermocline  is  observed.  The  deep  diurnal 
mixing  cycle  of  the  winter  mrmths  iq^arently  acts  to  reduce  the  buoyancy  flux  at  die  base 
of  the  mean  mixed-layo-. 

A  similar  feature  is  ^iparent  in  the  longterm  mean  {vofrle.  The  longtom 
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mean  was  significantly  positive  in  the  upper  50  m  and  negative  around  60  m. 
possibly  indicating  a  diurnal  momentum  flux  from  the  surface  layer  to  the  thermocline. 
It  is  impmtant  to  note  that  the  surface  lay^  increase  of  S  was  nearly  double  the  surface 
layer  increase  of  N  during  die  evening  time  bin.  This  is  also  reflected  in  the  longterm 
mean  proAle  of  A/?i^  which  shows  significantly  decreased  Ri  in  the  ui^r  SO  m  during 
the  evening  time  bin.  Investigators  using  TROPIC  HEAT  microstructure  data  indicated 
comparable  or  off-setting  diurnal  difierences  of  N  and  S  with  no  af^arent  diurnal 
differences  of  Ri  (Chereskin  et  al.,  1986;  Mourn  et  al..  1989).  The  Skyllingstad  and  Denbo 
(1994)  model  results  showed  increased  instantaneous  Ri  in  the  mixed  layer  during  the 
evening  and  decreased  instantaneous  Ri  in  the  mixed  layer  during  the  morning.  Similarly, 
preliminary  results  from  TIWE  data  (Lien  et  al..  1994)  have  indicated  increased  Ri  in  the 
mixed  layer  during  the  evening  and  decreased  Ri  in  the  mixed  layer  during  the  morning. 

The  reasons  for  the  i^arent  differences  between  the  longterm  mean  diurnal 
differences  of  Ri  computed  here  and  those  found  during  TIWE  and  in  the  Skyllingstad 
and  Denbo  model  results  remain  uncertain.  Examinatitm  of  Figure  6.1 1  shows  times  v^en 
ARi^  was  negative  and  agreed  with  the  TIWE  findings,  i.e.  higher  Ri  in  the  mixed  layer 
during  the  afternoon.  It  is  possible  that  siiKe  TIWE  was  ccmducted  during  die  mature 
phase  of  the  1991/1992  ENSO  warm  event,  the  observed  diurnal  cycle  ci  Ri  may  have 
indicated  anomalous  conditions.  Figure  6.11  shows  negative  ARi^  during  March  through 
May  1987,  which  coincided  with  die  mature  phase  of  the  1986/1987  ENSO  warm  event 
Strnie  differences  between  results  would  also  be  eiqiected  due  to  different  vertical  and 
tempOTal  resolution  of  the  data.  Whereas  diurnal  difference  in  Ri  shown  in  Figure  6.12 
include  both  time  when  the  instruments  (TRs  and  VACMs)  were  in  the  mixed  layer  and 
time  when  the  instruments  wee  below  the  mixed  laye,  the  highe  vertical  reolution 
velocity  data  during  HWE  enabled  Lien  e  al.  to  examine  these  time  sqMrately. 
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Longterm  Mean  Diurnal  Differences  of  IDV,  N.  S,  &  Ri 

AIDV^  AM«.  ASiM  ABI^ 

■I-/-  Std  liTor  ■«■/-  Std  Irror  -f/-  8td  bror  -•>/-  8td  trrer 


Figure  6.12.  Umgtenn  mean  vertical  prcrfiles  of  diurnal  diffnences  of  a)  isothenn 
diaplacement  variance,  AIDV^  b)  buoyancy  fiequency,  c)  velocity  ahear, 
and  D)  Richardson  number,  AR4.,  computed  over  the  period  November  1983-October 
1987.  Dashed  lines  indicate  ±one  standvd  error  of  die  mean. 
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F.  MIXED-LAYER  DEPTH  AND  THE  DIURNAL  CYCLE  OF  IDV 

In  section  S.6,  strong  correlation  was  observed  between  mixed-layer  dq>di  and 
IDV  at  both  35  m  and  60  m.  In  this  section,  relationships  between  mtmthly  mean  AIDV^ 
at  35  m,  60  m,  and  100  m  and  mmithly  mean  mixed-layer  dq>ths  for  each  of  the  four  six- 
hour  time  bins  and  monthly  mean  diurnal  ranges  of  mixed-layer  dq>th  are  examined. 
Table  6.3  lists  the  correlatitxi  coefficients  between  each  of  these  variables.  Statistically 
significant  correlations  were  observed  between  AIDV„  at  60  m  and  monthly  mean  mixed- 
layer  depths  for  each  of  the  evening,  night,  and  morning  time  bins.  These  correlations 
were  particularly  strong  during  the  night  and  morning  time  bins,  -0.532  and  -0.581. 
respectively.  These  negative  correlations  indicate  that  AIDV^  at  60  m  increased  during 
periods  when  the  mixed  layer  was  deep  and  decreased  during  periods  when  the  mixed 
layer  was  shallow.  In  other  words,  when  the  mixed  layer  was  deep,  indicating  either 
increased  surface  forcing  or  weakened  background  stability  conditions  or  bodi,  the  diurnal 
difference  of  the  high-frequency  temponture  variability  increased  at  60  m.  Correlatioos 
between  mixed-layer  dq>th  and  AIDV^  at  35  m  and  100  m  were  not  significant  at  the 
95%  ctmfidence  level. 

Significant  positive  correlatitms  woe  obsoved  between  AIDV^  at  35  m  and  60 
m  and  the  diurnal  range  of  mixed-layer  depth.  The  diurnal  range  of  mixed-layer  depth  can 
be  c<»sidered  a  measure  of  die  ni^tdme  mixing  intensity.  The  positive  correlations 
indicate  that  AIDV^  increased  when  die  mixing  intensity  increased.  This  suggests  a 
reladtmship  between  mixed  layer  dynamics  and  the  diurnal  cycle  of  the  higfa-frequmicy 
temperature  variability. 
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Table  6.3.  Coirelation  coefficients,  r.  between  monthly  mean  diurnal  differences  of  IDV 
(ATDV^  at  35  m,  60  m,  and  100  m  and  monthly  mean  mixed-layer  depths  (MLD)  for 
four  six-hour  time  bins  and  diurnal  ran^  of  mixed-layer  depth  (AMLD)  during 
November  1983-October  1987.  AIDV^  computed  as  0400-1000  IDV  minus  1600-2200 
IDV  normalized  by  mean  IDV.  AMLD  computed  as  monthly  mean  of  daily  maximum 
MLD  minus  daily  minimum  MLD.  Statistical  significance  at  the  95%  confidence  level 
was  calculated  based  on  the  integral  time  scale  between  time  series  (Davis.  1976). 
Statistically  signiticant  correlations  are  indicated  by  Bold  type. 


VARIABLES 

0400-1000 

1000-1600 

1600-2200 

95%c.i. 

r 

95%c.i.  r 

95%c.i. 

r 

95%c.i.  r 

MLD/AIDV_35 

.359 

-.197 

.363  -.242 

.354 

-.054 

.365  -.070  1 

.406 

-.532 

.417  -.581 

.389 

-.353 

.409  -.428  1 

/AEDV.IOO 

.268 

-.079 

.269  -.092 

.290 

-.007 

.280  -.015 

aMLD/AIDV_35 

.327 

.443 

H||[|Q[Q2E9 

.361 

.455 

/AIDV_100 

.295 

.164 

G.  DIURNAL  CYCLE  OF  WAVE  MOMENTUM  FLUX 

In  section  6.1,  the  diurnal  cycle  of  high-frequency  temperature  variability  was 
shown  to  be  confined  above  the  ewe  of  the  EUC,  and  in  most  cases  above  100  m.  To 
frulher  examine  the  diurnal  cycle  of  high-frequency  variability  in  the  upper  equatorial 
ocean,  15-minute  averaged  velocity  data  wore  used  to  estimate  die  high-frequency  wave 
momentum  flux  during  November  1984.  This  period  was  chosen  because  it  coincided  with 
the  microstructure  measurements  of  the  TROPIC  HEAT  1  experiment  and  because  the 
monthly  mean  profrle  of  IDV  shows  a  pronounced  diurnal  cycle  penetrating  well  into  the 
thermocline  (Figure  6.13).  During  November  1984,  the  velocity  shear  was  strong  and 
stratification  was  relatively  weak  resulting  in  low  Richardson  numbers,  Ri  less  than  unity, 
down  to  100  m.  Surface  winds  were  anomalously  stremg,  with  a  mean  wind  speed  of  - 
7.01  ms'*.  The  easterly  z<mal  wind  pseudo-stress  was  -51.89  mV^  (hfcPhaden  and 


McCarty,  1992).  The  monthly  mean  mixed-layer  depth  ranged  between  1 1.8  m  during  the 
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IDV  and  Ri  at  0®,  140®W  During  November  1984 


+/-  Sld  Deviation 
Kicliiirtlson  Ntiinlurr,  l.iii  Ui ) 


Richardson  Number,  LnCRO  Isotherm  Displacement  Variance  (m^ 


Figure  6.13.  Monthly  mean  vertical  prdiles  of  isotherm  di^lacement  variance  (IDV) 
during  November  1984  for  four  six-hour  time  bins  and  mcmthly  mean  Richardson  number 
(x-ofile.  Error  bars  indicate  ±1  standard  error  of  the  mean. 
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day  time  bin  and  22.S  m  during  the  morning  time  bin. 

Figure  6.14  shows  profiles  of  u'w'  for  each  of  the  four  six-hour  time  bins  for  the 
month  of  November  1984.  Positive  values  of  u'w'  were  observed  in  the  depth  range  45- 
120  m,  indicating  a  downward  flux  of  wave-induced  westward  momentum  from  the  wind- 
driven  surface  flow  to  the  EUC,  decelerating  die  eastward  flow.  It  appears  that  this 
downward  flux  of  westward  m<Mnentum  did  not  penetrate  the  EUC  core.  The  wave 

momentum  flux  divergence,  was  divergent  down  to  about  75  m  and  convergent 

below  75  m. 

Diurnal  difierences  of  u'w'  between  the  four  time  bins  were  not  significant  to  (me 
standard  error  of  the  mean,  excqit  at  25  m.  where  u'w'  during  the  day  time  bin  was 
greater  than  u'w'  during  the  mcxning  and  evening  time  bins.  Although  these 
indicate  a  vertical  transport  of  energy  fr(»n  the  surface  layo*  to  the  EUC  by  hi^- 
frequency  wave  momentum  flux,  it  appears  that  this  flux  was  not  related  to  diurnal 
forcing.  It  is  likely  that  the  observed  wave  momentum  flux  during  November  1984  was 
driven  by  mechanical  ftMoing  by  the  anomalously  large  surface  wind  stress.  It  is  not 
unreasonable  to  expect  seasonal  and  interaimual  variability  of  the  wave  momentum  flux, 
but  since  it  does  not  appear  related  to  the  diurnal  cycle  of  internal  wave  energy,  diese 
aspects  are  beyond  the  scope  of  diis  investigaticm. 
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Figure  4.14.  Montfily  mean  profiles  of  wave  momentum  flux  fw  four  six-hour  Hine  King 
during  November  1984.  Method  of  computation  described  in  section  3.2.8.  Error  bars 
indicate  ±one  standard  error  of  the  mean.  Units  are  cmV^. 
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H.  OTHER  POSSIBLE  SOURCES  OF  VARIANCE 


1.  Mooring  Motion 

Since  the  observations  of  high-frequency  temperature  variability  used  in 
the  foregoing  analyses  were  derived  frcun  moored  buoy  measurements  that  are  known 
to  be  subject  to  vertical  aiul  horizontal  motion,  it  is  not  unreasonable  to  assume  that  some 
of  the  temperature  and  velocity  variability  in  the  data  were  due  to  mooring  motion.  No 
indqyendent  measurements  of  mooring  motion  are  available  to  accurately  determine  the 
contribution  of  mooring  motion  to  such  errcM'.  Unlike  Boyd  (1989)  and  Levine  (1987), 
who  were  concerned  about  the  effect  of  mooring  motion  on  the  vertical  coherences  of  the 
internal  wave  spectrum,  the  ccmcem  here  is  whether  <x  not  mooring  motion  could  degrade 
the  observations  of  a  diurnal  cycle  of  high-frequency  temperature  variability.  The  most 
likely  causes  of  increased  mooring  motitm  are  wind  and  sea  conditions.  It  is  possible  that 
part  of  the  seasonal  and  interannual  fluctuations  of  IDV  that  have  been  described  are 
related  to  mooing  motion  caused  by  seasonal  and  interannual  fluctuations  of  wind  and 
sea  conditions.  However,  since  the  diurnal  variability  of  the  surface  winds  is  small  at  0”, 
140nv,  it  is  not  likely  that  the  observed  diurnal  variatioDs  of  IDV  were  related  to 
mooring  motion  caused  by  these  forces.  Another  possible  cause  of  diurnal  modulation  of 
the  temperature  variability,  either  directly  or  indirectly  by  mooring  mmion,  are  the 
gravitational  tides. 

2.  Gravitational  Tides 

Spectra  of  the  unfiltered  temperature  and  velocity  data  from  the  moorings 
have  fM-onounced  peaks  at  diurnal  and  semi-diurnal  tidal  frequencies.  Since  this  study 
focuses  on  the  diiunal  variability  of  high-frequency  wavelike  modmis,  eliminatitm  or 
signifrcant  reduction  of  the  diurnal  and  smni-diurnal  tide  from  the  data  without 
signifrcantly  reducing  the  diurnal  modulation  of  hi^-frequency  waves  is  critical.  Udal 
analyses  of  both  the  temperature  and  velocity  records  showed  10-50%  of  the  diurnal 
variance  was  removed  by  32  harmonic  crastituents.  These  analyses  were  dq)endait  on 
the  (^th  and  on  the  recm'd  length  examined,  which  varied  between  29  and  180  days. 
Accofding  to  diese  harmonic  analyses,  the  remaining  50-90%  of  the  diurnal  variance  was 
determined  to  be  unrelated  to  the  barotropic  ticfe.  Although  useful  as  an  indicator  of  the 
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amount  of  tidal  variance,  the  harmonic  analyses  are  subject  to  errors  associated  with  non- 
stationarity  of  the  data.  Fcm*  the  most  part,  the  voiical  temperature  structure  is  affected  by 
the  baroclinic  tide,  but  not  the  barotropic  tide.  Velocity,  on  the  other  hand,  is  affected  by 
both  baroclinic  and  barotropic  tides. 

The  primary  advantage  of  the  procedure  used  in  this  study  to  filter  the 
diurnal  and  semi-diurnal  tide,  described  in  section  3.1,  was  that  it  was  applied  locally  and 
was  therefore  less  subject  to  errors  associated  with  the  non-stationarity  of  the  internal 
tides.  A  disadvantage  of  the  procedure  is  die  inability  to  easily  quantify  the  role  of  the 
filter  in  reducing  the  variance.  Although  it  is  possible  that  some  of  the  diurnal  modulation 
of  variance  from  the  residual  fields  was  tidally  forced,  the  fact  that  the  resulting  residual 
variance  had  a  well-detined  vertical  structure  suggests  that  the  barotropic  tide  was  mostly 
removed.  In  their  tidal  analysis  of  velocity  records  in  the  eastern  equatorial  Pacific, 
Weisburg  et  al.,  (1987)  found  that  the  barotropic  tide  accounted  for  78.8%  of  the  tidal 
energy.  The  remaining  baroclinic  tide  was  random,  isotropic,  and  without  mode 
preference.  Their  partition  between  barotropic  and  baroclinic  tidal  energy  is  midway  in 
the  range  of  50-90%  given  by  Wunsch  (1975)  for  the  global  ocean.  Thus,  it  appears  that 
most  the  tidal  energy,  that  due  to  the  barotropic  tide,  was  reasonably  well  accounted  for. 

The  remaining  concern  then  lies  with  determining  the  impact  of  the 
baroclinic,  (m*  internal  tide,  on  calculations  using  velocity  data,  such  as  velocity  shear, 
Richardson  number,  and  wave  momentum  flux.  The  internal  tides  are  simply  a  fmm  of 
internal  wave  with  quasi-tidal  periodicity  (Wunsch,  1975).  Theory  and  observations 
suggest  that  internal  tides  are  generated  by  the  passage  of  the  barotropic  tide  over 
bathymetric  features,  diough  it  is  not  possible  to  say  with  certainty  that  any  of  several 
mechanisms  dominates  (Wunsch,  1975).  Variatimis  in  the  density  and  velocity  field  near 
the  generation  region  produce  variations  in  the  response  of  the  internal  tide.  Away  fnxn 
the  source,  the  internal  tide  is  reflected  and  refracted  by  the  density  and  velocity  field 
though  which  it  propagates,  with  the  resulting  internal  tide  being  highly  variable  in  time 
and  space  (Levine  and  Richman,  1989). 

In  evaluating  the  possible  importance  of  the  intonal  tide  on  the 
interpretation  of  results  presented  here,  die  surrounding  bathymetry  was  compared  with 


the  surrounding  bathymetry  for  the  analysis  of  Weisberg  et  al.  (1987)  from  mooring  data 
at  0“,  110"  W.  By  comparison,  the  bathymetry  near  0",  140"W  is  relatively  featureless. 
A  few  small  seamounts  are  located  about  250  km  to  the  west.  The  shallowest  of  these  is 
at  a  depth  of  about  2800  m,  compared  with  the  surrounding  4(X)0  m  abyssal  plain.  The 
nearest  islands  are  the  Marquesas  located  about  900  km  to  the  south.  Aside  from  the 
small  seamounts,  the  nearest  large  badiymetric  features  in  a  zonal  direction  are  about 
1700  km  to  the  west.  The  bathymetry  surrounding  the  moored  array  used  in  Weisberg 
et  al.  (1987)  is  complex,  with  the  east  Pacific  Rise  and  numerous  fracture  zones  ir  close 
proximity.  If  internal  tides  are  primarily  generated  by  barotrq^ic  tides  interacting  with 
bathymetric  features,  as  evidence  suggests,  it  is  not  umvasonable  to  assume  that  the 
internal  tide  at  0",  140"W  would  be  no  greater  than,  and  possibly  much  less  than,  that 
described  by  Weisberg  et  al.  (1987). 

Examination  of  diurnal  differences  between  four-year  means  of  IDV  and 
the  dynamic  stability  parameters  provide  further  indication  that  diurnal  gravitational  tides 
are  not  likely  to  have  significantly  modified  the  interpretation  of  results.  The  dominant 
diurnal  tidal  constituents,  Kl,  PI,  and  01,  have  pmods  of  23.93  hours,  24.07  hours,  and 
25.82  hours,  respectively.  These  periods  cause  the  constituents  to  shift  relative  to  the  solar 
day,  with  the  Kl  constituent  retarding  3.93  minutes  per  day,  the  PI  constituent  advancing 
3.95  minutes  per  day,  and  the  01  constituent-advancing  almost  two  hours  per  day.  Each 
year  the  Kl  and  PI  constituents  cycle  completely  through  a  solar  day,  effectively 
averaging  out  shorter  period  contributions  to  the  solar  diurnal  cycle.  Calculations  of 
means  over  the  four-year  record  average  over  four  complete  Kl  and  PI  tidal  cycles  and 
over  a  hundred  complete  01  tidal  cycles.  Hie  fact  that  significant  diurnal  variations  wo-e 
observed  in  these  longterm  mean  calculations  and  that  the  diurnal  variations  display 
significant  vertical  structure  indicates  that  these  diurnal  cycles  wo'e  ftxced  by  non-tidal 
processes. 
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VII.  CONCLUSIONS 


A.  SUMMARY  AND  DISCUSSION 

In  diis  research,  the  diurnal  cycle  of  high-frequency  temperature  fluctuations 
attributable  to  turbulence  and  internal  waves  in  the  ui^)er  central  equatorial  Pacific  Ocean 
was  examined  over  seasonal  and  interannual  time  scales.  Using  moored  temperature  and 
velocity  data  in  the  upper  300  m  from  OP,  14QPW  during  the  period  November  1983- 
October  1987,  IDV,  AIDV,  and  the  vertical  structure  of  AIDV  were  examined  under 
surface  forcing  and  dynamic  stability  conditions  that  varied  over  diurnal,  seasonal,  and 
interannual  time  scales.  Using  monthly  mean  values  of  AIDV  as  a  proxy  for  turbulence 
and  internal  waves,  a  pronounced  diurnal  cycle  of  overturning  and  vertical  displacement 
was  observed  at  35  m  and  60  m,  and  thk  diurnal  cycle  was  found  to  vary  seasonally  and 
interannually  in  response  to  the  varying  surface  forcing  and  background  cmiditions. 

The  diurnal  cycle  of  turbulence  and  high-frequency  internal  waves  was 
characterized  by  nighttime  increases  in  IDV  with  maximum  valties  occurring  during  die 
morning  (0400-1000)  time  bin  and  minimum  values  occurring  during  die  evening  (16(X)- 
2200)  time  bin.  This  diurnal  cycle  of  IDV  is  consistent  with  microstructure  observadmis 
over  a  few  short  periods,  showing  that  IDV  is  a  useful  {voxy  for  turbulence  and  internal 
wave  energy  during  periods  when  microstructure  measurements  are  not  available. 

Since  diurnal  variations  in  the  surface  wind  stress  were  small,  IDV  is  presumed 
to  have  been  diumally-modulated  by  diurnal  variadons  of  the  surface  buoyancy  flux, 
which  were  predominandy  govoned  by  the  strtmg  diurnal  cycle  of  the  shortwave 
radiation.  In  addition  to  regulating  the  diurnal  cycle  of  IDV,  the  surface  buoyancy  flux 
induced  pronounced  diurnal  variations  of  the  mixed-layer  depth  and  associated  dynamic 
stability  in  the  upper  50-100  m.  As  with  AIDV,  the  diurnal  differences  of  mixed-layer 
dq>th,  N,  S,  and  Ri  wore  largest  between  the  morning  and  evening  time  bins.  As  eiqiected, 
N  and  5  in  the  ufqjer  50  m  increased  during  the  afternoon  lestratificatimi.  An  interesting 
finding  is  that  S  increased  by  a  factor  of  two  more  than  N,  resulting  in  decreased  Ri  in 
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the  upper  50-60  m  during  the  afternoon.  This  four-year  mean  diurnal  cycle  of  Ri  was 
contrary  to  sh<»1-term  microstructure  observations  showing  no  apparent  diurnal  variation 
of  Ri  (Chereskin  et  al.,  1986;  Mourn  et  al..  1989)  or  an  of^XMite  diurnal  cycle  of  Ri  in  the 
mixed-layer  (Lien  et  al.,  1994).  There  were  a  few  times  during  the  four-year  study  period, 
most  notably  during  the  mature  phase  of  the  1986/1987  ENSO  warm  event,  when  Ri's 
computed  from  the  mooring  data  were  larger  during  the  evening,  in  agreement  with  the 
Lien  et  al.  (1994)  results,  which  were  based  on  observations  taken  during  the  mature 
phase  of  the  1991/1992  ENSO  warm  event.  Differences  in  Ri  results  could  also  be  related 
to  large  differences  in  the  vertical  and  tempcxal  resolution  of  the  data.  The  longterm  mean 
diurnal  cycle  of  Ri  rq)orted  here  includes  observations  when  the  instruments  were  in  the 
mixed  layer  and  observations  when  the  instruments  were  below  the  mixed  layer. 

In  the  four-year  mean,  the  diurnal  cycle  of  turbulence  and  internal  waves,  as 
proxied  by  AIDV,  decayed  sharply  with  dq)tfa  below  60  m  and  was  not  significant  at  100 
m  tx  deq)er.  Likewise,  vertical  correlations  between  AIDV^  at  35  m  and  AIDV^  at  60 
m,  100  m,  140  m,  200  m,  and  300  m  indicated  sharp  decay  with  depth  below  60  m.  This 
is  not  to  say  that  AIDV  never  penetrated  below  60  m.  In  fact,  on  monthly  time  scales  a 
weak  diurnal  cycle  of  IDV  was  observed  at  1(X)  m  during  about  29%  of  the  rntmths 
examined.  Nevertheless,  the  combined  longterm  means  and  vertical  correlations  indicate 
that  the  diurnal  cycle  of  turbulence  and  internal  waves  generally  did  not  extend  below  the 
EUC.  Thereftxe,  it  appears  unlikely  that  vertical  propagatitm  of  diumally-forced  waves 
from  the  surface  layer  to  the  deq>er  waters  below  the  EUC  ctmtributed  significantly  to 
the  equatorial  momentum  balance,  as  has  been  hypothesized  by  Dillon  et  al.  (1989)  and 
Wijesekera  and  Dillon  (1991). 

Although  it  cannot  be  determined  with  certainty  whether  these  observations  of  IDV 
rq)resent  high-firequency  internal  waves,  local  shear  instabilities,  turbulent  mixing,  or 
some  cmnbination  of  these  processes,  the  vertical  structure  of  AIDV  and  the  relati(»ish4>s 
between  IDV  and  surface  forcing,  dynamic  stability  conditions,  and  mixed-layer  dqpth 
provide  imptxtant  clues.  The  vertical  structure  of  AIDV  indicates  that  diumally-fnoed 
v«lically-propagating  internal  waves  q^jear  not  to  be  transporting  significmit  eneigy 
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from  the  surface  layer  through  the  EUC  to  the  deep  ocean.  If  vertically-propagating 
internal  waves  were  generated  in  or  near  the  base  of  the  mixed  layer,  it  appears  that  they 
did  not  prt^agate  far  vertically  before  breaking  or  triggering  local  shear  instabilities  and 
turbulent  mixing. 

Seasonally,  both  IDV  and  AIDV  reached  annual  nuxima  during  the  boreal  winter 
and  minima  during  the  boreal  spring.  These  seasonal  extrema  coincided  well  with 
maximum  wind  speed  and  easterly  zonal  wind  during  the  boreal  winter  and  minimum 
wind  speed  and  easterly  zonal  wind  during  the  boreal  spring.  Similarly,  mixed-layer  dq>th 
and  the  dimnal  range  of  mixed-layer  dq)th  had  maxima  during  the  winter  and  minima 
during  the  spring.  AddititMially,  the  seasonal  cycle  of  the  dynamic  stability  of  the  tqjpo' 
ocean  was  characterized  by  a  narrow,  high-5,  low-/?/  rpgion  at  abmit  35  m  during  the 
txxeal  spring  and  a  broad,  weaker-5,  low-/?/  region  at  about  the  60-m  dq>th  dunng  the 
boreal  winter.  In  summary,  the  boreal  spring  was  characterized  by  weak  wind'i,  a  shallow 
mixed  layer,  decreased  diurnal  range  of  mixed-layer  dq>th.  a  shallow  hi^-5,  low-/?/ 
region,  and  low  IDV  and  AIDV.  The  boreal  winter,  on  the  other  hand,  was  characterized 
by  strong  easterly  winds,  a  deep  mixed  layer,  increased  diurnal  range  of  mixed-layer 
depth,  a  relatively  deep,  broad,  low-/?/  regicm,  and  high  IDV  and  AIDV. 

The  fact  that  IDV  and  AIDV  at  35  m  were  at  seasonal  minima  i^e  the  dynamic 
stability  at  35  m  was  weakest  sug^ts  that  IDV  did  not  rq)r68ent  only  local  shear 
instabilities.  In  other  words,  the  high-5,  low-/?/  region  at  35  m  during  the  boreal  spring 
is  expected  to  have  been  most  likely  to  produce  local  shear  instabilities.  If  local  shear 
instabilities  were  generated,  they  were  not  reflected  in  the  monthly  mean  IDV  and  AIDV. 
As  the  seasonal  mixed  layer  and  low-/?/  region  deq)ened  to  about  35  m  and  60  m, 
respectively,  during  the  boreal  summer,  IDV  and  AIDV  both  increased  at  35  m,  but 
remained  low  at  60  m.  The  low  dynamic  stability  at  60  m  did  not  correspmid  widi 
increased  IDV  or  AIDV.  It  is  notewothy  that  IDV  and  AIDV  increased  at  35  m,  ^diich 
was  at  die  base  of  the  mixed  layer.  During  the  boreal  winter,  whmi  the  mixed  layer 
deepened  to  the  seastmal  maximum  of  about  50  m  and  the  dynamic  stability  increased, 
IDV  and  alDV  increased  to  seasonal  maxima  at  35  m  and  60  m,  respectively.  In  diis 
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case,  temperature  measurements  at  35  m  were  well  within  the  mixed  layo-,  as 
defined  in  this  study.  Although  efforts  were  taken  to  remove  IDV  values  at  times  when 
the  instrument  was  within  the  mixed  layer,  it  is  possible  that  the  remaining  measurements 
of  IDV  at  35  m  during  the  winter  months  were  biased  by  turbulent  mixing  in  the  mixed 
layer,  rather  than  internal  waves  or  local  shear  instabilities.  The  wdnter  IDV  measurements 
at  60  m,  on  the  other  hand,  were  generally  below  the  base  of  the  mixed  layer,  so  it  is  not 
clear  whether  AIDV  represented  increased  IDV  caused  by  internal  waves,  turbulent 
mixing,  or  a  ctmibination  of  these  processes. 

Over  interannual  time  scales.  IDV  and  AIDV  showed  some  axreqxmdence  to 
ENSO  variability,  but  this  relationship  was  not  always  clear.  The  value  of  IDV  at  35  m 
was  observed  to  be  low  during  the  1984  ENSO  cold  q>isode,  steadily  inovasing  during 
the  ENSO  build-up  phase,  and  high  during  the  1986/1987  ENSO  warm  q>isode.  The 
values  of  AIDV^  at  35  m  showed  a  similar  pattern.  The  values  AIDV^e  at  35  m  woe 
similar,  but  values  were  higher  during  the  ENSO  build-up  phase  than  during  the  mature 
phase  of  the  1986/1987  warm  q[>isode.  Data  gi^  in  the  60  m  record  made  examination 
of  the  interannual  variability  more  difficult.  The  values  of  IDV  and  AIDV  at  60  m  were 
low  during  the  1984  cold  q>isode.  high  during  the  build-up  phase  and  beginning  of  the 
1986/1987  warm  event,  and  low  during  the  later  half  of  the  warm  q>isode. 

Monthly  means  of  IDV  at  35  m  and  60  m  over  the  four-year  record  were  well 
correlated  with  wind  speed,  zonal  wind,  mixed-layer  dq>th,  diurnal  range  of  mixed  layer 
dq>th,  latent  heat  flux,  net  surface  buoyancy  flux,  and  evaporation.  Generally,  IDV  in  the 
upper  60  m  was  significantly  correlated  with  the  stratification  in  tee  tq)per  100  m. 
Cwrelations  between  IDV  and  velocity  shear  and  Ri  were  generally  weak. 

Monthly  means  of  AIDV  at  35  m  and  60  m  were  well  correlated  with  wind  q)eed, 
zonal  wind,  shortwave  radiation  flux,  sensible  heat  flux,  diurnal  range  of  mixed-layer 
depth,  and  buoyancy  frequency  in  tee  u|;q;}er  60  m.  In  addition,  AIDV  at  60  m  was  well 
cmrelated  wite  mixed-layer  dq>te.  Each  of  tee  correlatitms  was  greater  at  60  m  tean  at 
35  m,  possibly  a  reflection  of  the  shorter  record  lengte  and  data  gq>s  during  teree  of  the 
boreal  summers.  The  values  of  AIDV  at  100  m  were  not  well  correlated  wite  mixed-layer 
depth  or  any  of  the  surface  fncing  <x  dynamic  stability  conditimis. 
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What  do  these  correlatioiis  between  monthly  mean  values  of  IDV  and  dlDV  and 
the  surface  forcing  and  background  stability  conditions  reveal?  Since  both  IDV  and  AIDV 
at  35  m  and  60  m  had  the  strongest  ctmrelatitms  with  the  surface  wind  fncing,  it  seems 
reasonable  that  diurnal  variations  of  IDV  and  AIDV  were  directly  related  to  the  amount 
of  mechanical  energy  input  into  the  ocean  at  the  surface  by  the  wind  stress.  Similarly, 
mixed-layer  depth  was  strongly  related  to  tl»  surface  wind  stress,  in  agreement  with  most 
mixed  layer  theories.  The  complete  explanation  of  the  relationship  between  IDV  and 
AIDV  and  mixed-layer  depth  is  less  certain.  It  is  not  unreasonable  to  suggest  that 
processes  controlling  mixed-layer  depth  might  have  also  contributed  to  the  generation  of 
the  observed  high-frequency  temperature  variability.  In  addition  to  mechanical  fcxcing  by 
the  wind  stress  and  thermodynamic  forcing  by  the  surface  buoyancy  flux,  mixed-layer 
d^th  was  strongly  correlated  with  the  dynamic  stability  parameters,  N,  S,  and  Ri.  The 
relationships  between  dynamic  stability  and  IDV  and  AIDV  were  not  consistent  on 
monthly  time  scales.  Although  shortwave  radiation  was  significantly  correlated  with  AIDV 
at  35  m,  the  lack  of  correlation  widi  mean  IDV  is  not  readily  e]q[>lained. 

High-frequency  wave  momentum  flux  during  November  1984  indicated  a 
downward  flux  of  westward  momentum  frmn  the  wind-driven  surface  flow  to  the  EUC, 
but  q>paiently  not  through  the  core  of  tibe  EUC.  During  November  1984,  wave 
momentum  flux  did  not  q>pear  to  vary  significantly  on  diurnal  time  scales.  However, 
these  wave  momentum  flux  estimates  were  based  on  15-minute  averaged  VACM  data 
which  effectively  removed  some  of  the  hi^-frequency  variability  of  interest. 

B.  CONCLUSIONS 

The  diurnal  cycle  of  high-frequency  temperature  variability  in  the  upper  layers  cS 
the  central  equatorial  Padfic  Ocean  was  examined  using  moored  buoy  temperature  and 
velocity  data  at  0",  140^,  in  the  ui^)er  300  m  during  the  period  November  1983- 
October  1987.  The  most  noteworthy  findings  from  this  analysis  are  listed  below. 

»>  A  prtmounoed  diurnal  cycle  of  IDV  was  observed  at  35  m  and  60  m,  with  a 
maximum  during  the  morning  hours  and  a  minimum  during  the  evening  hours. 
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*■  The  diurnal  cycle  of  IDV  decayed  nq>idly  with  dq>th  below  60  m  aiui  was  not 
generally  significant  at  100  m  or  deqier,  indicating  that  it  is  unlikely  that  diumally-forced 
waves  propagate  energy  frcmi  die  surface  layer  through  the  EUC  core. 

>  The  above  two  findings  agree  widi  microstructure  observations  and  suggest  that 
AIDV  is  a  useful  proxy  for  turbulence  and  internal  waves  in  the  equatorial  Pacific  Ocean. 

*■  The  diurnal  cycle  of  turbulence  and  internal  waves  at  both  35  m  and  60  m,  as 
proxied  by  AIDV,  varied  seasonally,  widi  a  maximum  magnitude  during  the  boreal  winter 
and  a  minimum  magnitude  during  the  boreal  spring. 

*■  The  diurnal  cycle  of  turbulence  and  internal  waves  at  both  35  m  and  60  m,  as 
proxied  by  AIDV,  was  strongly  correlated  with  wind  speed,  zonal  wind,  mixed>layer 
depth,  diurnal  range  of  mixed-layer  dqith,  and  stratification  in  the  upper  60  m. 

»>  Pronounced  diurnal  cycles  of  the  dynamic  stability  parameters  were  obsoved  in 
the  surface  layer,  with  smaller  Ri  during  the  evening  dian  during  the  morning. 
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